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PREFACE 


In this book authors bring out the innovative applications of 
matrices defined, described and developed by them. Here they 
do not include the natural product on matrices newly described 
and defined by them in the book on ‘natural product x, on 
matrices’. 

This book is organized into seven chapters. The first one is 
introductory in nature. In the second chapter authors give the 
unique and new way of analyzing the data which is time 
dependent. We construct three types of matrices called Average 
Time Dependent data matrix (ATD matrix), Refined Time 
Dependent Data matrix (RTD matrix) and Combined Effective 
Time Dependent Data matrix (CETD matrix). The authors 
describe the working of this new type of matrix model by 
working with the real world transportation problem. It is proved 
this new model is effective and elegant. At this juncture the 
authors deeply acknowledge the unostentatious help rendered by 
Dr. Mandalam. 

In chapter three the authors for the first time define the new 
concept of matrices with linguistic variables. These linguistic 
matrices are used in the place of fuzzy matrices which serve as 
the dynamical system for FCM and FRM models. By this 
method the solution themselves are linguistic terms so 
mathematical interpretation is not needed. 


In chapter four authors have used super matrices in the 
following ways. They are used in super fuzzy models like super 
FCMs, super FRMs super FAMs etc. Super linear algebra using 
these super matrices pave way to super eigen values and super 
eigen vectors, infact in almost all applications of matrices in 
linear algebra. The special mention is made to the Leontief 
open production super model and closed super model. By this 
method easy comparison and bulk working is possible. 

We have also used super matrices with entries from finite 
fields in the construction of several types of super codes. 

In chapter five we define interval matrix and matrix interval 
using the natural class of intervals, that is intervals of the form 
[a, b] where a < b or a> b or a and b not comparable. Using 
these interval matrices we have built interval linear algebra and 
these interval matrices when used as stiffness matrices of any 
mechanical problem easily yield a solution. 

Suggestions and literature how we have used and 
constructed DSm matrices and DSm super matrices is given. 
However the authors have constructed DSm vector spaces and 
DSm super vector spaces. The eigen values and eigen vectors of 
these spaces will be refined labels and refined label vectors, 
which comprises chapter six of this book. 

In the final chapter we have defined the new notion of 
bimatrices and n-matrices. Using these new concepts we have 
constructed linear algebra of type I and linear algebra of type II. 
Using these matrices we can build n-eigen values and n-eigen 
vectors. 

Further these n-matrices are used in the construction of 
multi expert fuzzy models like n-FCMs, n-FRMs, n-FAMs etc. 
Finally we have built several types of n-codes using these n- 
matrices whose entries are from finite fields. 

We thank Dr. K.Kandasamy for proof reading and being 
extremely supportive. 


W.B.VASANTHA KANDASAMY 
FLORENTIN SMARANDACHE 
INDRA VENKATBABU 


Chapter One 


INTRODUCTION 


This book gives the contributions of the authors in the 
innovative applications of matrices to various fields. 

The authors were first to use matrix theory in the analysis 
of raw data. This is described in chapter two of this book. The 
method of applying matrices in this way makes one know the 
results and conclusions by just looking at it. Thus the final data 
is given a graph representation. 

Secondly the authors have constructed matrices where 
entries are just linguistic terms which are orderable. For more 
about this refer [53]. Using these linguistic values we construct 
matrices. Infact these linguistic matrices are used in studying 
fuzzy linguistic models [53]. Thus we have in the first place 
constructed linguistic matrices and used them in mathematical 
models to study socio economic problems. 

The authors have used super matrices in the first place to 
build super linear algebras. Almost all classical results of linear 
algebras are derived for the super linear algebras [45]. Using 
these super matrices we have built several types of super fuzzy 
models and these super fuzzy models can yield the result in a 
shorter time and comparision of the results are ready made [48]. 
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Further these matrices are used in the construction of 
algebraic super codes of different types. For more about 
supercodes please refer [50]. 

We build matrix of refined labels and using these matrices 
we construct DSm vector spaces of refined labels. Several 
properties about DSm vector spaces build using matrix of 
refined labels in elaborately discussed. This is yet another 
means of applying matrix theory to DSm vector spaces of 
refined labels [46]. 

We also construct DSm super vector space of refined labels. 
For these new structures the notion of supermatrices is used. In 
the supermatrix if we replace the real entries or any other entries 
by refined labels L we call them as supermatrix of refined labels 
and also they form vector spaces. This is carried out in chapter 
six of this book [52]. 

Finally we build a natural class of intervals {[a, b] | a > b 
a<bora=b,a,b € QorC or Z or C(Z,), Z, or R}. Using 
these natural classes of intervals we build interval matrices. 
These interval matrices behave like ordinary matrices as the 
entries are from this special natural class of intervals. Further 
study of stiffness matrix, load vector, mass matrix and damping 
matrix using interval for in this case lower and upper bound 
matrices is based on interval algebra. Also we find for square 
interval matrices, interval eigen values and interval eigen 
vectors [55]. 

Finally we define the notion of bimatrices and n-matrices. 
This is yet another innovative method of applying matrices. 
These n-matrices (n > 2) can be used in the construction of 
fuzzy n-models like n-FCMs, n-FRMs, n-BAMs etc., and their 
mixed structures [54, 56]. Finally we use these n-matrices in the 
construction of different types of n-codes [49]. 


Chapter Two 


AVERAGE TIME DEPENDENT (ATD) 
DATA MATRIX 


The raw data under investigation is classified under four 
broad heads viz; total number of passengers, total collection, 
number of trips, and hourly occupancy. These four broad heads 
form the columns of the matrices. The time periods of the day 
are represented by the rows of the matrices. Estimating the 
utility rate of a route is a five-stage process. In the first stage, 
we give a matrix representation of the raw data. Entries 
corresponding to the intersection of rows and columns are 
values corresponding to a live network. The initial M x N 
matrix is not uniform i.e., the number of individual hours in 
each time period may not be the same. So, in the second stage, 
we in order to obtain an unbiased uniform effect on each and 
every data so collected, transform this initial matrix into an 
Average Time Dependent Data matrix (ATD matrix). To make 
the calculations easier and simpler, we in the third stage, using 
the simple average techniques, convert the above time 
dependent data matrix into a matrix with entries e, where, ej € 
{-1, 0, 1}. 
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We name this matrix as the Refined Time Dependent Data 
matrix (RTD matrix). The value of the e, corresponding to each 
entry is determined in a special way. At the fourth stage, using 
the refined time dependent data matrices, we get the Combined 
Effect Time Dependent Data matrix (CETD matrix) which gives 
the cumulative effect of all these entries. In the final stage, we 
obtain the row sums of the combined effect time dependent data 
matrix. We have written a Java program, which estimates all 
these five stages. [41] 


Using the raw data available for any transport corporation, 
we analyze the raw data via matrices and (i) predict the 
maximum utilization time period (peak hours) of a day; and (11) 
estimate the overall utility rate of the routes. Thus, to be more 
precise our chief problem here is prediction of the peak hours 
and estimation of the most utilized routes using the raw data 
available from any transport organization. We have established 
that the results which we have predicted using the raw data on 
the route 18R which ply from Parrys to Dharmaraja koil 
coincides with the peak hours and the utilization rate of the 
route which can be obtained from the estimation of the observed 
data. Thus, our analysis not only predicts the peak hours and the 
maximum utility routes but also estimates for each route, time 
periods where buses need not be operated as operation of buses 
in those time periods will result in total loss. 


By total loss, we mean that the money collected from the 
passengers at that specified time period may not be even 20% of 
the cost of petrol and service charges spent on plying that 
service. So, even if the concern using the predictions on peak 
hours gets a profit of say 70%, if they do not comply with 
stopping of buses at the non-utilized hours, the transport 
corporation will still result only in loss. Hence, our approach 
via matrices not only suggests the peak hours and the maximum 
utilized routes but also predicts time periods and the routes in 
which the services should be totally avoided or reduced from 
their usual schedule in order to save themselves from loss. 
Thus, by adhering to our predictions the corporation might 
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curtail the loss incurred by the operation of services in the non- 
utilized time periods. 


The predictions can be explicitly seen even by a layman just 
by the observation of the combined effect time dependent data 
matrix mentioned earlier, which will be described and discussed 
in detail in the following. 


Hence, using these results the Transport Corporation can 
operate more number of buses at the peak hours and on the 
routes that are highly utilized and also stop operating buses at 
odd hours or reduce the frequency at the odd hours of the day. 
Thus, the Transport Corporation can achieve a meager gain or 
curtail the loss by using these predictions. 


We have considered the city of Madras, where the Pallavan 
Transport Corporation (PTC) is one of the transport 
organizations that caters to the demands of the community. The 
Pallavan Transport Corporation is a_ service oriented 
organization run by the state government for the people of the 
city and the neighbouring panchayats in the immediate vicinity 
of the city. It has a fleet of 1343 buses, and covers an overall of 
283 routes and satisfies the needs of more than 75% of the 
urban transit commuters. Apart from the regular services, it has 
the partial services, the night services and the special services. 


The regular services are of four types: normal, limited stop 
service, point-to-point service and the express service. Although 
the source, the destination and the path of travel is common to 
all the regular services, these services considerably differ in the 
travel time, the travel fare, the comfort and the convenience of a 
passenger. The basic difference observed in the above types of 
services is primarily due to the number of intermediate stops the 
vehicles halt in their course of travel. 


The city transport corporation connects a vast network inter- 
linking various routes. The network has a set of nodes and a set 
of links connecting these nodes. Buses connect these nodes and 
ply over the links and help in passenger and _ freight 
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transportation. In the city network, we observe that some routes 
have a high patronage almost throughout the day, while other 
routes are busy only at particular hours of a day and there are 
still some routes, which have only very meager passenger 
support throughout the day. Here, we not only predict the high 
utility and peak hours but also suggest the stoppage of meager 
passenger supported services. 


To test the efficiency of the program, we have applied the 
algorithm to a few and established that the observed peak hours 
coincide well with the estimated peak hours for each route 
individually. 


In the first stage of the problem, we split the total hours of a 
day into various time-periods viz. morning, morning peak, 
evening. Each of these time periods consists of individual hours 
of the day viz. {H,, Ho, H3, ..., H22, H23, H24}, where H; refers to 
the hour ending of the day. The working hours of the transport 
corporation under consideration is from {Hs, Ho, Hy, ..., Hai, 
Hy}. The data obtained from the transport corporation is for 
each hour ending. Hence, the time interval between H; and H;+, 
is 60 minutes. If we have all the related data viz. the total 
collection, the total number of passengers, the occupancy rate 
and the number of trips etc., for say every 30 minutes or for 
every 15 minutes, then the accuracy in the identification of the 
peak hours will further increase. 


The finer the intervals are divided the better is the 
prediction. This is also established in this chapter. Each time 
period consists of a set of H; individual hours, for example: the 
morning peak time period consists of {Hg, Ho, Hio}. These 
time-periods can vary from one analyst to another and from one 
route to another. Even, if we have each individual hour as a 
time period, still the total number of time periods must not 
exceed to total working hours of a day. 


Let M represent the total number of time periods under 
consideration. The utility rate of a time period is influenced by 
the various attributes acting on them. Let N represent the 
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number of attributes considered for the analysis. The various 
time periods are treated as rows and the various attributes are 
treated as columns of a matrix. The raw data corresponding to 
the intersection of each time period and attributes are treated as 
the entries of a matrix. Thus, from the raw data, we obtain the 
initial M x N matrix. 


The number of individual hours in each time period might 
not be the same. Hence, in the second stage, in order to obtain 
consistency in the initial matrix we use the “hourly concept”. 
The entries corresponding to the intersection of each of the time 
periods and the attributes of the initial matrix are transformed so 
that each new entry corresponds to the hourly rate viz. total 
passengers per hour, total collection per hour, total number of 
trips per hour, total occupancy per hour etc. 


Thus, we convert the above initial matrix into the Average 
Time Dependent Data matrix (ATD matrix) i.e., [aij }Mav. 


In the third stage, we use the average and the standard 
deviation to convert to above average time dependent data 
matrix into a matrix with entries ej, ej; € {-l, 0, 1}, where 1 
represents the i" row and j represents the j" column. We call 
this newly formed matrix as the Refined Time Dependent Data 


matrix (RTD matrix) i.e., [ej ]Mxn. 


The value of the entry e, corresponding to each intersection 
is determined from an interval. This interval is obtained strictly 
by using the average and the standard deviation calculated from 
the raw data. The choice of the interval made by us might not 
be able to identify the accurate peak hour of the route; hence, 
we introduce and define a parameter a, which enables us to get 
a better solution. 


We calculate the mean py, and the standard deviation 6; for 
each attribute j, j = 1, 2, .... N using the data given in the 
average time dependent data matrix. For varying values of the 
parameter a, where a € [0, 1]. We follow a rule and determine 
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the value of the entry ej in the refined time dependent data 
matrix. 


For each of the attributes j (j = 1, 2, ..., N) we have the rule: 


If aj < (uj -a* 0}) then ei = —l; 
else, 

ita; ¢ (yj — &.* oj, > & *0;): then‘e; = 0; 
else 

if ay > (uj +o * o}) then ej = 1; 


here ‘*’ denotes the usual multiplication. 


Thus, for different values of a, we obtain different refined 
time dependent data matrices. The main purpose of introducing 
the refined time dependent data matrix is only to minimize the 
time involved in performing the simple arithmetic calculations 
and operations on the matrix. 


In the fourth stage, we bring in the notion of the Combined 
Effect Time Dependent Data matrix (CETD matrix) 1.e., [¢i]mxn, 
which gives the combined effect of all the refined time 
dependent data matrices obtained by varying the parameter a. 
In the final stage, we add up the rows of the combined effect 
time-dependent data matrix. The overall time period utility of a 
route is obtained by inferring the row sums of a combined effect 
time-dependent data matrix. The highest positive value is taken 
as the highly utilized time period of a route and next lower 
value is taken as the next peak hour for the same route. Thus, 
for a particular route, we grade the utilization rate of the 
different time periods. 


The computation starts by computing the ATD matrix from 
the initial M x N matrix. Then mean and standard deviations 
are computed for the ATD vector. Based on mean, standard 
deviation and ATD vector, the RTD and CETD matrices are 
computed either physically or using a programme in Java or 
C++. 
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It is pertinent to mention here this technique of analysing 
raw data can be used in case of any problem which is time 
dependent. This is a very new way of applying matrix theory to 
find solution for a collected raw data. 


Parameters play a vital role in estimating the peak hours of 
a route. This is clearly illustrated for the route 18JJ and the 


route 17D by varying the parameter a. First, consider the route 
18JJ. We have the following: 


The initial M x N matrix obtained from the raw data is as 
below: 


[ 81.00 63.00 3.00 147.20 | 
925.00 157.00 16.00 1394.55 
1040.00 181.00 15.00 1585.50 
1159.00 271.00 13.00 1916.05 |. 
1596.00 400.00 18.00 2615.30 
1673.00 451.00 16.00 2889.60 
| 322.00 160.00 6.00 522.50 | 


The corresponding average time dependent data matrix is 


40.50 31.50 1.50 73.60 | 
308.33 52.33 5.33 464.85 
346.67 60.33 5.00 528.50 
386.33 90.33 4.33 638.68 |. 
532.00 133.33 6.00 871.77 
557.67 150.33 5.33 963.20 
161.00 80.00 3.00 261.25 | 
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The refined time dependent data matrix corresponding to a 
single value of a = 0.1 is 


ee es ee 
“1.1 1-1 
0-11 0 
f £04 
1111 
fe, Ae As el 
|-1 -1 -1 -1] 


The corresponding row sum of the above refined time 
dependent data matrix is given by 


Row sum for the first row is -4 
Row sum for the second row is —2 
Row sum for the third row is 0 
Row sum for the fourth row is 3 
Row sum for the fifth row is 4 
Row sum for the sixth row is 4 
Row sum for the seventh row is -4. 


From the row sums of the above refined time dependent 
data matrix, we observe that time period corresponding to the 
fifth row i.e., time period {Hi6, Hi7, Hig} and the sixth row i.e., 
time period {Hjo, Hx, H2:} are the peak hours of the route 
followed by the fourth row viz. time period {Hj3, Hy4, His}. 
Thus, graphically, we have this as below: 
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Graph depicting the peak periods of the 
Route 18JJ fora =0.1 


Row sum of the CETD matrix 
ak woe OA NWR 


Hours of the day 


The refined time dependent data matrix corresponding to 
the value of a = 0.4 follows: 


-1-1-1-l 


0-1 1 0 
0-11 0 
00 0 0 
1 111 
1111 
-1 0 -1-1 


From the above, we have the row sums as given below: 


Row sum for the first row is -4 
Row sum for the second row is 0 
Row sum for the third row is 0 
Row sum for the fourth row is 0 
Row sum for the fifth row is 4 
Row sum for the sixth row is 4 
Row sum for the seventh row is —3. 


From the row sums of the above refined time dependent 
data matrix, we observe that time period corresponding to the 
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fifth row i.e., time period {Hj., Hj7, Hig} and the sixth row Le., 
time period {Hj 5, Hx, H2;} are the peak hours of a day. Time 
periods {H), Hg, Ho}, {Hho, Hi, Hy} and {Hy3, Hya, Hys} have a 
row sum of zero indicating that these time periods are neither 
the peak hours of a day nor the non-peak hours of a day. Hence, 
the number of trips made in these time periods should be kept 
unaltered for the concern to maintain the profit. 


Thus, graphically, we have this as below: 


Graph depicting the peak hour for the 
Route 18JJ for a = 0.4 


Ny 
wo 
- 
ua 
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Row sum of the CETD matrix 
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Hours of the day 


The refined time dependent data matrix corresponding to 
the value of a = 0.7 follows: 


[—fal—t =] 
0-10 0 
000 0 
000 0 
[ 1-24 
1101 
-1 0 -1-1|] 
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From the above, we have the row sums as given below: 


Row sum for the first row is 4 
Row sum for the second row is —1 
Row sum for the third row is 0 
Row sum for the fourth row is 0 
Row sum for the fifth row is 4 
Row sum for the sixth row is 3 
Row sum for the seventh row is —3. 


From the row sums of the above refined time dependent 
data matrix, we observe that time period corresponding to the 
fifth row i.e., time period {Hj6, Hi7, Hig} is the first peak hour 
of the day followed by the sixth row 1.e., time period {Hj9, Ha, 
H)}}. Time periods {H,, Hg, Ho}, {Hho, Hu, Hy} and {H3, Hy4, 
H,;} have a row sum of zero indicating that these time periods 
are neither the peak hours of a day nor the non-peak hours of a 
day. Thus, graphically, we have this as below: 


Graph depicting the peak hour for the route 
18 JJ for a = 0.7 


Row sum of the CETD matrix 
oa kd © & Oa DY WRG 


Hours of the day 
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The refined time dependent data matrix corresponding to 
the value of a = 0.9 follows: 


ee ie ees 
0000 


0 
0 
1 
1 
0 


From the above, we have the row sums as given below: 


Row sum for the first row is -4 
Row sum for the second row is 0 
Row sum for the third row is 0 
Row sum for the fourth row is 0 


Row sum for the fifth row is 4 
Row sum for the sixth row is 3 
Row sum for the seventh row is —3. 


The row sums of the above matrix indicates time periods 
corresponding to the fifth row i.e., {Hj6, Hi7, Hig} is the first 
peak hour of the day followed by the sixth row ie., {Hi9, Hx, 
Hp}. Time periods {H), Hg, Ho}, {Hho, Hi, H)>} and {Hy3, Hua, 
Hs} have a row sum of zero indicating that these time periods 
are neither the peak hours of a day nor the non-peak hours of a 
day. Thus, graphically, we have this as below: 
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Graph depicting the peak hour for the 
route 18 JJ for a = 0.9 


Row sum of the CETD matrix 
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From the above analysis, we observe that the peak hours of 
a route vary from one time period to another with the change in 
the value of the parameter from 0 to 1. The row sums matrix 
obtained in the above cases were specific to only one value of 
the parameter. But the combined effect time dependent matrix 
for all the values of a € [0, 1] is given below: 


The combined time dependent data matrix for all a € [0, 1] 
is 


—10 —10 —10 —10 
-l1 -8 6 -2 
-6 4 0 
1 0O 3 
10 10 10 10 
10 10 6 10 
-9 -1 -9 -9 


The row sum of the above combined time dependent data 
matrix is 
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Row = 1, row sum = —-40 
Row = 2, row sum = —5 
Row = 3, row sum = —2 
Row = 4, row sum = 7 
Row = 5, row sum = 40 
Row = 6, row sum = 36 
Row = 7, row sum = —28. 


From the row sums of the above matrix, we observe that the 
time period corresponding to the fifth row i.e., time period {Hy., 
H)7, His} is the first peak hours of the day followed by the sixth 
row ie., time period {Hj 9, Hx, H2,} and the fourth row 
corresponding to the time period {H)3, Hy4, His}. All the other 
time periods have a row sum of negative value indicating that 
these time periods are the non-peak hours of a day. The first 
time period has the least negative value indicating that services 
can be curtailed to the greatest extent in this time period. 


Thus, graphically, we have the peak hour as depicted below: 
From the graph below, we observe that the peak hours as 
obtained from the combined time dependent data matrix, given 
by the cumulative effect of all the values of a € [0, 1], gives the 
true peak hours of a day. These results coincide well with the 
observed hours of the route 18JJ. 


Graph depicting the peak hour for the 
Route 18 JJ for a € [0, 1] 


Row sum of the CETD Matrix 
Oo 


Hours of the day 
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We give a grading of the peak hours since they are essential 
and helpful to the transportation sector in a number of ways: 


1. It is efficient in identifying the highly utilized time 
periods along each route and enables the transport 
sector to operate more number of buses during these 
hours. 


2. The transport sector identifies the poorly utilized time 
periods and makes necessary arrangements to overcome 
the loss incurred due to the services scheduled in these 
time periods. They might either discard all the trips 
made during these hours or operate minimum number 
of services. 


3. The row sum of the combined time dependent data 
matrix of each route helps in identifying the overall 
utility rate of each route. In other words, if all the row 
sums of the combined time dependent data matrix are 
positive, then we conclude that the particular route is 
highly utilized throughout the day. 


4. A positive value of the row sum of the combined time 
dependent data matrix indicates the maximum 
utilization of that specific time period. 


5. A negative value of the row sums of the combined time 
dependent data matrix indicates that the passenger 
patronage and the collection in those time periods are 
very less. Hence, the transport sector can curtail the 
number of services totally, partially or operate mini 
buses to overcome the loss. 


6. A zero value of the row sum of the combined time 
dependent data matrix indicates that the corresponding 
time periods are neither the peak hours of a day nor the 
non-peak hours of a day. 
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It is to be noted here that as the total number of time periods 
(M) increases, the accuracy in identifying the utilization rate of 
the different time-periods also increases. We substantiate our 
claim by taking two different values for the total number of time 
periods and this is carried out in the illustration I. 


Illustration I (Route 18R) 
The number of time periods is five (M = 5). 


Here, we take five time periods viz. early morning C;:{Hg, 
H;}; morning C): {Hg, Ho, Hho, Hirt; noon C3: {Hp, H33, Hy, 
H)s}; evening C4: {Hio, Hy7, Hig, Hyo}; and night Cs: {H9, Hp, 
H»}. We obtain the row sum of the combined time dependent 
data matrix as follows: 


The combined time dependent data matrix corresponding to 
M = Sis given below: 


-10 —8 -10 -10 
6 —2 6 5 
2 -1 6 1 
10 10 6 10 
4-6 0 -3 


The row sums of the corresponding combined time 
dependent data matrices is given by: 


Row sum for the time period early morning = —38 
Row sum for the time period morning = 15 

Row sum for the time period noon = 8 

Row sum for the time period evening = 36 

Row sum for the time period night = —13 


The number of time periods is seventeen (M = 17). 
Here, we take 17 time periods viz. early morning C, : {He}; 


morning C2: {H7}; C3 : {Hg}, C4: {Ho}; early noon Cs: {Hjo}; 
Co: {His}, Cr: {Hi}; noon Cs: {His}; Co: {His}, Cyo: {His}; 
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evening Cy;{Hi6}; Cio: {Hi7}, Cis : {His}; late evening Cy,: 
{Hy}, Cys: {Hoo}, Cio: {Ho} and night C7: {H >}. The row 
sums of the combined time dependent data matrices are 
obtained as follows: 


The combined time dependent data matrix corresponding to 
M = 17 is as follows: 


-10 -7 O- -10 
0 -5 O- -5 
4 2 0 2 
7 -1l O 5 
0 —2 O 0 
1 -3 10 2 
3 0 0 2 
2 -l1 O 0 


10) 272-0 =10) 


The corresponding row sum of the above combined time 
dependent data matrix is given as below: 


Row sum for the time period C;: {Hg} =—32 
Row sum for the time period C,: {H7} =—27 
Row sum for the time period C3: {Hg} =—10 
Row sum for the time period C4: {H»} = 4 
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Row sum for the time period Cs: {Hjo} = 11 
Row sum for the time period C,: {H,;} =—2 
Row sum for the time period C7: {Hj2} = 10 
Row sum for the time period Cg: {H,3} =5 


Row sum for the time period Co: {H,4} = 1 
Row sum for the time period Cio: {Hjs} =—32 
Row sum for the time period C,,: {Hj5} = 31 
Row sum for the time period Cy2: {Hj7} = 11 


Row sum for the time period C)3: {Hs} = 23 
Row sum for the time period C,4: {Hj9} =-8 
Row sum for the time period C15: {H29} = 30 
Row sum for the time period Cys: {H2;} =—9 
Row sum for the time period C,7: {H22} = —27. 


From the above Illustration I, on the analysis of the row 
sums of the combined time dependent data matrices of M = 5, 
we notice that the time period corresponding to row = 4 namely 
the evening hours {His, Hi7, His, Hio} are the busy hours 
followed by row = 2 namely the morning hours {Hg, Ho, Hho, 
Hy}. 


On analyzing the row sums of the combined time dependent 
data matrices of M = 17, we observe that row = 11 viz. evening 
{H6} corresponds to the busiest hour of the route followed by 
row = 15 viz. late evening {H20}, row = 13 viz. evening {Hjs}. 


Thus, we observe that as the number of time periods 
increases, the accuracy in identifying the busier hours of a route 
also increases. Also, we can observe the added effects of one or 
more attributes on the various time periods. 


Thus, for a € [0, 1] with a step function of a = 0.1, the 
peak hours of the route are Hi, Hoo and Hig. Time period 
corresponding to Hs is poorly utilized and hence can be 
discarded. 
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The same problem is analyzed and the combined time 
dependent data matrix for varying values of a € [0, 1] with a 
step function of a = 0.05 is given below: 


The combined time dependent data matrix, the row sum and 
the corresponding graph is given as below: 


The combined time dependent data matrix 


[-10 -4 -10 -10] 
-10-10 0 -10 
=1'=10° 0. =10 
G4 0) 5 
10 -3 0 10 
a] er <0, «if 
2: =F 10." § 
7 21-0 A 
A 23° OY 4 
=10 =4.=10.<10 
10 2 10 10 
10 -2 0 10 
10 7 O 10 
-§ 10 -10 -9 
10 -1 10 10 
—9 -5 0 -7 
|-10 -10 0 -10 
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The row sums of the above combined time dependent data 
matrix is 


Row = 1, row sum = —34 
Row = 2, row sum = —30 
Row = 3, row sum = -21 
Row = 4, row sum = 9 


Row = 5, row sum = 17 
Row = 6, row sum = —5 
Row = 7, row sum = 10 
Row = 8, row sum = 10 


Row = 9, row sum = 2 
Row = 10, row sum = —34 
Row = 11, row sum = 32 
Row = 12, row sum = 18 


Row = 13, row sum = 27 
Row = 14, row sum =—17 
Row = 15, row sum = 29 
Row = 16, row sum = —21 
Row = 17, row sum = —30. 


We illustrate and describe the problem using 18R service. 


We have applied our algorithm to a single route and 
obtained a comparison of the observed and the calculated data. 
We have taken the route 18R plying from Parrys to Dharmaraja 
Koil and the data was obtained from the Pallavan Transport 
Corporation Limited, Madras, India. 


Here, we take 17 time periods viz. early morning C;:{H¢}; 
morning C2: {H7}; C3: {Hs}, Ca: {Ho}; early noon Cs: {Hijo}; Co 
> {Hii}, C7: {Hiz}; noon Cg : {Hi3}; Co: {His}, Cio : {His}; 
evening C;,: {Hie}; Cio: {Hi7}, Ci3: {Hig}; late evening Cy,4: 
{Ho}, Cis {Hoo}, Ci6 4 {Ho}; and night Ci7 : {H»} and also 
four attributes viz. A; : total population; A, : hourly occupancy; 
A3: number of trips and Ay: total collection. 
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We view this problem from the aspect of both the operators 
as well as the users. 


The initial M x N Matrix is as follows: 


96.00 64.00 1.00 156.95 | 
71.00 21.00 2.00 103.55 
222.00 39.00 2.00 244.90 
269.00 63.00 2.00 353.50 
300.00 70.00 2.00 392.80 
220.00 60.00 2.00 328.30 
241.00 53.00 3.00 353.65 
265.00 78.00 2.00 348.65 
249.00 70.00 2.00 327.85 |. 
114.00 63.00 1.00 151.60 
381.00 91.00 3.00 526.40 
288.00 73.00 2.00 407.05 
356.00 112.00 2.00 515.15 
189.00 389.00 1.00 252.75 
376.00 78.00 3.00 569.40 
182.00 59.00 2.00 261.90 
67.00 20.00 2.00 81.35 | 
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The average time dependent data matrix is as follows: 


| 96.00 
71.00 
222.00 
269.00 
300.00 
220.00 
241.00 
265.00 
249.00 
114.00 
381.00 
288.00 
356.00 
189.00 
376.00 
182.00 
| 67.00 


64.00 
21.00 
39.00 
63.00 
70.00 
60.00 
53.00 
78.00 
70.00 
63.00 
91.00 
73.00 
112.00 
389.00 
78.00 
59.00 
20.00 


1.00 
2.00 
2.00 
2.00 
2.00 
2.00 
3.00 
2.00 
2.00 
1.00 
3.00 
2.00 
2.00 
1.00 
3.00 
2.00 
2.00 


156.95 
103.55 
244.90 
353.50 
392.80 
328.30 
353.65 
348.65 
327.85 |. 
151.60 
526.40 
407.05 
515.15 
252.75 
569.40 
261.90 
81.35 


The refined time dependent data matrices corresponding to 
different values of alpha are calculated and from them we obtain 
the combined time dependent data matrix. 


Since, the refined time dependent data matrices are obtained 
from the average time dependent data matrices they have a size 
equivalent to that of the average time dependent data matrices 
but, they have entries only in {—1, 0, 1}. For different values of 
the parameter we have different refined time dependent data 


matrices. 
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The refined time dependent data matrix corresponding to 
the value of a = 0.1 is 


-1 -1 -1 -l 
-l1 -1 0 -l 
0 -l O -l 
1 -l 0 1 
1 -1 0 1 
0 -l1 0 O 
1 -1l 1 1 
1 0 0 1 
1 -1 0 0O 
-1 -1 -1 -l 
I t 1 1 
1 -l1 0 1 
1 1 0 1 
-l1 1 =-1l -l 
1 0 1 1 
-1 -1 0 -1 
-1 -l1 0 -l 


The refined time dependent data matrix corresponding to the 
value of a = 0.2 is 


| 
_ 
| 
_ 
| 
— 
| 
— 


-1 -1 0 -1 
0 -l O -l 
1 -l1 0 1 
1 0 0 1 
0 -l1 0 O 
0 -1 1 1 
1 0 0 1 
1 0 0 0O 


| 
_ 
| 
_ 
| 
= 
| 
_ 
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1 0 1 1 
1 0 0 1 
1 1 0 1 
-l1 1 -l1 -l 
1 0 1 1 
-l1 -l1 0 -1l 
-l1 -1 0 -l 


The refined time dependent data matrix corresponding to the 
value of a = 0.3 is 


| 
— 
=) 

| 
fa 

| 
— 


| 
— 
oO > > = = 


or Oo OF KF SO 
i=) 


— 
oorrcoooCoC lO 
- | 
— 
= — | 
— 


| 
— 

| 
_ 

| 
nan 
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The refined time dependent data matrix corresponding to the 
value a = 0.4 is 


-1 0 -1 -l 
-1 -1 0 -l 
0 -l O -l 
1 0 0 
1 0 1 
0 0 0 
0 1 0 
0 0 0 
0 0 0 


| 
— 
ooroceocjcoooocooo 
| 
— 
| 
— 


| 
— 
| 
— 
oor 
i) 


The refined time dependent data matrix corresponding to the 
value a = 0.5 is 


| 
ue 
° 
SOR S oe 6: oO 1. 
| 
ne 


| 
— 

| 
aan 


34 | Innovative Uses of Matrices 


1 0 1 1 
1 0 0 1 
1 0 0 1 
0 1 -1l 0 
1 0 1 1 
0 0 0 0 


The refined time dependent data matrix corresponding to the 
value of a = 0.6 is 


| 
oO 

| 

| 


-1 -1 0 -l 
0 0 0 0 
0 0 0 0 
1 0 0 0O 
0 0 0 0 
0 0 1 0 
0 0 0 0 
0 0 0 0 
=l O° =L =1 
1 0 1 1 
1 0 0 1 
1 0 0 1 
0 1 -l 0 
1 0O 1 
0 0 0 


I 
— 
| 
— 
oo Fr 
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The refined time dependent data matrix corresponding to the 
value of a = 0.7 is 


| 
es 
° 
corooooo | 
| 
un 


ooroceocjcoooocoocoo 
| 
— 
| 
— 


| 
— 

| 
— 
oor 


The refined time dependent data matrix corresponding to the 
value of a = 0.8 is 


So So He 16.0 o:o ) 
| 
es 


cooocooococoocoo 


| 
_ 
— 
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1 01 1 
0 0 0 0O 
1 0 0 1 
0 1 -l O 
1 0 1 1 
0 0 0 O 
-1 0 0 -l 


The refined time dependent data matrix corresponding to the 
value of a = 0.9 is 


-1 0 -1 -l 
-1 0 0 -!l 
0 0 0 O 
0 0 0 O 
0 0 0 O 
0 0 0 O 
0 0 1 O 
0 0 0 O 
0 0 0 O 
-1 0 -1 -l 
1 01 1 
0 0 0 O 
1 0 0 1 
0 1 -l O 
1 01 1 
0 0 0 O 
-1 0 0 -l 
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The refined time dependent data matrix corresponding to the 
value of a = 1.00 is 


-1 0 -1 -l 
-1 0 0 -l 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 1 0 
0 0 0 0 
0 0 0 O 
-1 0 -1 -l 
1 0 1 1 
0 0 0 0 
1 0 0 1 
0 1 -l 0 
1 0 1 1 
0 0 0 0 
-1 0 0 -l 
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The combined time dependent data matrix 


-10 -7 O- -10 
0 -5 O  -5 
4 2 0 2 
7 -l1 O 5 
0 2 O 0 
1 -3 10 2 
3 0 O 2 
2 -l O 0 


The row sums corresponding to the above combined time 
dependent data matrix are: 


Row sum for the time period C,: {H.} =—32 
Row sum for the time period C2: {H7} =—27 
Row sum for the time period C3: {Hg} =—10 
Row sum for the time period C4: {Ho} = 4 

Row sum for the time period Cs: {Hjo} = 11 


Row sum for the time period Cg: {H;,} =—2 
Row sum for the time period C7: {H)2} = 10 
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Row sum for the time period Cg: {H)3} =5 
Row sum for the time period Co: {H,4} = 1 


Row sum for the time period Cjo:{H)5} =—32 
Row sum for the time period C;;:{Hj6} = 31 


Row sum for the time period Cj2:{H)7} = 11 
Row sum for the time period C;3: {Hig} = 23 
Row sum for the time period C4: {Hj9} =—8 
Row sum for the time period C15: {H9} = 30 
Row sum for the time period Cy6:{H2,} =—9 
Row sum for the time period C17: {H22} =—27 


We observe that peak hours of a route differ for varying 
values of alpha. Hence, we form the combined effect time 
dependent data matrix which gives the combined effect of all 
the values of alpha in the interval [0, 1]. Thus, the highest value 
of the row for the corresponding time period identified from the 
row sums of the combined time dependent data matrix is 
observed as the best peak period of the particular route. The 
next value of the row sum in the small matrix is observed to be 
the next peak period of the same route. In this manner, we rank 
the different time periods of a route. 


In the above estimated method, we notice that the time- 
period corresponding to row = 11 viz. evening {Hi6} 
corresponds to the busiest hour of the route followed by row = 
15 viz. late evening {Hy}, and row = 13 viz. evening {Hjs}. 
Thus, we rank the peak periods for the route 18R to be as 
evening C,,:{Hjo6}; late evening C;5: {Hx} and evening Cj; : 
{Hjs} in the order of first, second, third respectively. 


We now compare the predictions with the ridership method. 


To establish the efficiency of our method, we compare our 
results with the ridership, approach followed by Dr. J. 
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Vasudevan in his research on the rationalization of city bus 
routes. 


He had taken Madras as a case study where, he had 
discussed about the rationalization of the routes and had studied 
for a total of 90 routes. A small section of it deals with the 
average hourly bus passenger load for the city. For each hour 
ending, say 5 am. 6 am., ..., 21 p.m., 22 p.m., he had 
calculated the total number of passengers who had traveled 
during that particular hour. 


He had also calculated the percentage of the passenger load 
at any hour against the total strength of passengers who had 
traveled that day. This is clearly depicted by the following table. 


Determination of the peak hour based on the 
Passenger load per hour 


Hour ending | Passengers /hour | Percent to total 
6 96 2.47 
7 71 1.83 
8 222 =e 
S 269 6.92 
10 300 Telz 
1 220 5.66 
12 241 6.20 
13 265 6.82 
14 249 6.41 
15 114 293 
16 381 9.80 
17 288 7.41 
18 356 9.16 
19 189 4.86 
20 376 9.68 
21 182 4.68 
ze 67 Liz 

Total 3886 100 
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We have applied the matrix method to this data in order to 
obtain a comparison and establish that our method is better than 
the ridership method. From the last column of the above table, 
we identify the first peak period as the row corresponding to the 
greatest percent to the total, i.e., 9.80 corresponding to the hour 
ending 16, 1.e., evening 4 p.m., followed by 9.68 corresponding 
to the hour ending 20 i.e., night 8 p.m. 


The third peak period is the percent 9.16 corresponding to 
the hour ending 18 i.e., evening 6 p.m. In this fashion, we rank 
the peak periods of the route based on the passenger load. 


Thus, we have established that our matrix method gives us a 


better result that the ridership method of Vasudevan. 


Comparison Table depicting the peak hours 
as estimated by the three methods 


Data as Identification of the peak and non-peak 
obtained periods from the ridership method, the 
from the revenueship method and comparison of the 
transport two methods with the matrix method 

corporation 
Hour Observed Observed Matrix 
ending Method I: Method IT: Method: 
Passengers | Revenue/hour | tow sum of 
/hour the CETD 


Matrix 
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376 (II)* 


13 265 348.65 5 

14 249 327.85 1 

15 114 151.6 32 
16 381(1)* 526.4 (II)* 31 (D* 
17 288 407.05 11 

18 356 (III)* 515.15 (II1)* 23 (III)* 
19 189 252.75 


569.4 (I)* 


30 (II)* 


*denotes the ranking of the time periods as given within parenthesis 


From the above table, we obtain a comparative study of the 
observed methods and the estimated method. From the last 
column of the above table, we observe that there are rows that 
have a negative value. 

These negative values indicate that the corresponding time 
periods are comparatively poor in performance i.e., in the 
passenger usage and in the collection. Hence, the services in 
these time periods can be discarded or the frequency of the 
services can be reduced. 

This is one of the innovative means of applying matrices in 
analysis problems which has only raw data. Thus we have not 
given the graphs of them. 


Chapter Three 


FUZZY LINGUISTIC MATRICES AND THEIR 
APPLICATIONS 


In this chapter we introduce the new notion of fuzzy 
linguistic matrices and use them in fuzzy linguistic models and 
apply them to social problems. We need for this the operations 
on the fuzzy linguistic matrices. We briefly describe the 
operations on them. Let M be any matrix if the entries of M are 
taken from the fuzzy linguistic set we call such matrix M to be a 
fuzzy linguistic matrix for more refer [53]. We describe 
operations on them. Let L denote the collection of fuzzy 
linguistic terms. 


Let x = (aj, a2, a3, a4, as) be a fuzzy linguistic row matrix 
(vector) with entries from L. 


t—. te, 
X = (a, a2, a3, a4, as) = | a, 


is the fuzzy linguistic column matrix (vector). 
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That is if x = (0, bad, good, fair, very bad) be the fuzzy 
linguistic row matrix, then transpose of x is denoted by 


0 
bad 
x' = (0, bad, good, fair, very bad)'=| good 


fair 


very bad | 


which is a fuzzy linguistic column matrix. Suppose we have 
product of a fuzzy linguistic row matrix x; it can be operated 
with min or max operation with a fuzzy linguistic column 
matrix y if and only if x is a 1 x t row matrix then y must be a 
t x 1 column matrix. 


We show how the ‘min’ operation looks like. 


Suppose 
y=| . | and x = (by, by, ..., b) 


then min {min {x, y}} = min {min {a;, b,}, min {a, bp}, ..., 
min {a,, b,}}. 


Likewise we can have {max; max} operation; {max; min} 
operation and {min; max} operation. All these will be 
illustrated by examples. 


Let x = (Xj, X2, ..., X20) and y = (yi, yo, ..., Y20) with Xi, yi € 
L; 1 <i<20. Lbea fuzzy linguistic matrix. 


max {max (x, y)} = max {max (X;, yi), Max {X, Yo}, ..., 
max {X20, Y20} }. 
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Let x = (Xi, ». Ora X15) and i (yi, Y25 se ee5 Yis) 
max (min {x, y}) = max {min {x), y;}, min {Xx2, yo}, ..., min 
{X15, Vis} }- 


Likewise min (max {x, y}) = min {max {x,, yj}, 
max {Xo, Y2}, -+-, Max {X15, Vis} }. 


We will see how for a typical problem the four operations 
give different sets of answers. 


Let 


x= and 


| good 


y = (bad, 0, good, very bad, best, fair, best, better) 
be two fuzzy linguistic matrices. 
min {min (y, x)} 
= min {bad, 0, fair, 0, best, bad, 0, better} = 0 
max {min (y, x)} 


= max {bad, 0, fair, 0, best, bad, 0, better} = best, 
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min {max (y, z)} 


= min {good, bad, good, very bad, best, fair, 
best, good} 


= very bad and 
max (max {y, x}) 
= max {good, bad, good, very bad, best, fair, best, good} 
= best. 


Suppose we are interested in finding min {x, y} then we 
have 


min {x, y} = 


bad 0 good very bad good fair good better | 
bad 0 bad very bad bad bad bad bad 
bad 0 fair verybad fair fair fair fair 
00 O 0 i | oe 0 
bad 0 good very bad best fair best better | 
bad 0 bad very bad bad bad bad bad 
00 0 0 0 O O 0 
bad 0 good very bad good fair good better | 


This is the way a fuzzy linguistic column matrix of order 
8 x 1 is multipliced with a 1 x 8 fuzzy linguistic row matrix 
with min operation gives a 8 x 8 fuzzy linguistic matrix. 
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We can also find max {x, y} = 


[good good good good best good good good | 
bad bad fair bad best bad _ bad better 
fair fair fair fair best fair fair better 
bad 0 fair 0 best bad 0 _ better 

best best best best best best best best | 
bad bad fair bad best bad _ bad better 
bad 0 fair 0 best bad 0 better 
[ good good good good best good good good | 


We see clearly max {x, y} # min {x, y}. 


Further 0 dominates in min and best dominates in max. 
When an expert wants to boost the results on the positive side 
he can use the max operation. If min operation is used it gives 
the worst state of affairs. 


Thus we see we get different results for these four types of 
operations. According to need one can use any one of the 
operations. 


Now we find the transpose of a 6 x 5 fuzzy linguistic matrix 
M. 


[ good _ best bad fair best 
0 bad — worst 0 very bad 
bad = fair ~=veryfair bad better 
Let M= 
worst 0 good better 0 
0 fair 0 worst best 


| good fair best worst 0 | 


Now the transpose of this fuzzy linguistic 6 x 5 matrix M 
denoted 
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good 0 bad worst 0 good | 
best bad fair 0 fair fair 
M'=| bad worst very fair good 0 best 
fair 0 bad better worst worst 
best very bad better 0 best 0 


Likewise if 


| best 0 bad 
fair bad 0 

0 very fair good 
N = | good 0 best 
best — worst 0 

0 best good 
[good better best | 


is any fuzzy linguistic 7 x 3 matrix, to find N‘. 


best fair 0 good best 0 good 
N'=| 0. bad very fair 0 worst best better |. 
bad 0 good best 0 - good best 


We see N' is a3 x 7 fuzzy linguistic matrix. 


Now we can find the product of two rectangular linguistic 
matrices M and N if M is a n x t matrix then N must be at x m 
matrix then only MN is defined how ever NM is not defined in 
this case and MN is an x m fuzzy linguistic matrix. 


However product of any n x n matrix with itself is always 
defined. 


We will illustrate these two situations by some examples. 
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Let 


good bad 0 worst 
bad 0 good best 
0 good bad 0 
good bad good bad 


M= 


000 0 
in (min (M, M)) 000 0 
= min (min (M, = 
00 0 0 
000 0 


We see if ‘0’ occurs atleast once in every row and atleast 
once in every column then min (min (M, M)) = (0). 


Now we find 
good bad bad bad 
good good good bad 


bad bad good good] 
good good good good 


max (min (M, M)) = 


We now find the value of 


good good good best 

best best best best 
max (max (M, M)) = , 

good good good best 


good good good best 
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0 bad bad 0 
bad 0 bad bad 
bad bad 0 bad] 
bad bad bad _ bad 


min (max (M, M) = 


We see max {max (M, M)} gives an extreme or better 
values and min {min {M, M}} gives an extreme low values, 
where as min max non negative values and max min more 
positive values. 


Such four types of operations can be used as per need of the 
problem. 


Now if M is an x m matrix and N is a m x t matrix we can 
find min (min {M, N}), min (max {M, N}), max {max (M, N)} 
and max (min {M, N}). 


We will illustrate these four types of operations. 


Let 
low 0 high 0 very low 
_| 0 high 0 low high 
high OO medium 0 low 
low low 0 medium 0 
and 
| low 0 medium 0 
high low 0 high 
N= 0 very high low 0 
medium 0 high low 
| very low medium 0 high 
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be two fuzzy linguistic matrices associated with temperature of 
an experiment. 


00 0 0 
Clean Gane NY) 00 0 0 
early min (min {M, = : 
f 00 0 0 
00 0 0 

low high low v.low 


{min (M, N)} high medium low _ high 
max {min (M, = . ; 
low medium medium low 


medium low medium low 


high veryhigh high high 
high veryhigh high high 
high veryhigh high high | 
high veryhigh high high 


max {max (M, N)}= 


vlow 0 0 low 
; 0 low 0O 
Now min {max (M, N)} = i 9 0 1 
ow 


We see this four types of operations gives four types of 
4 x 4 fuzzy linguistic matrices. 


Now we can also find linguistic row matrix with a square or 
a rectangular fuzzy linguistic matrix which is compatible. 


Consider X = (fast, slow, very slow, just fast, very fast) to 
be the fuzzy linguistic row matrix. 
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Take 
slow medium fast 
very slow slow medium 
M= fast slow fast 
fast medium slow 
| just fast slow __ very slow | 


be a5 x 3 fuzzy linguistic matrix. Consider min {min {X, M}} 
= (very slow, very slow, very slow). 


Consider max (min {X, M}) = (just fast, medium, fast). 
Consider min (max {X, M}) = (slow, slow, medium). 


To find max (max {X, M}) = (very fast, very fast, very fast). 


Here we describe the notion of Fuzzy Linguistic Cognitive 
Models (FLCM). 


For Fuzzy Linguistic Cognitive maps please refer [53]. 


Example 3.1: Suppose we are interested in studying the child 
labour problem. Let (C;, Co, ..., Cs) be six attributes / concepts 
associated with it. 


C,; - Child Labour 

C, - Good Teacher 

C3; - School Drop out 

Cy, - Poverty 

Cs - Public encouraging child labour 


Let L = {0, often, + often, -often, very much, much, not that 
much, little, very little, more etc.}. 
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So the state vectors as well as the related fuzzy linguistic 
matrices take their values from the set L. Also the vertices of 
the fuzzy linguistic graph take their values from L. 


Now using the experts opinion we have the following fuzzy 
linguistic graph. 


+ often 


Now using this linguistic graph we have the following fuzzy 
linguistic matrix M. 


C, Cs C, C, C. 

Cc 0 0 0 +often + very much 
M= C, 0 0 -often 0 0 
C,| +often -—much 0 0 0 
C, 0 0 +often 0 0 
C,|+verymuch 0 0 0 0 


Now our state vector takes values from L. The min and 
max operations on L are defined as follows: 
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Let min {0, a;} = 0 and max {0, a;} =a; for all a; € L. 
a; € L min {aj, a;} = a; and max {aj, a;} =a}. 
min {a;—a;} =— a; and max {a;, —a;} = a; for all a; € L. 
For instance min {often, very often} = often 
max {often, very often} = very often 
min {much, often} = often, 
max {much, often} = much 
min {—much, often} =— much 
and max {—much, often} = often. 
Like this operations on L are performed. 


Now we find x|M using as before min {min (a;, m,j)} where 
X1 = (a1, a, ..., a6) and M = (mj); my, a € L, 1 <i,j <6. 


Thus x;M = (+ often, 0, + often, + often, 0) x 


0 0 0 +often + very much 
0 0 -often 0 0 
+often -—much 0 0 0 
0 0 +often 0 0 
+verymuch 0 0 0 0 


= (+ often, 0, + often, + often, + often) 
leading to a fixed point. 


Example 3.2: Consider the problem of finding the seven transit 
system which includes the level of service and the convenience 
factors. We have the following eight attributes. 


C,; - Frequency of the service along a route 
C, -  In-vehicle travel time along the route 
C; - Travel fare along the route 


C, - Speed of the vehicles along the route 
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Cs - Number of intermediate points in the route. 

Cs - Waiting time 

C; - Number of transfers in the route 

The fuzzy linguistic terms associated with C;, C2, ..., C7 as 


well as the problem are L = {0, often, always, a little, much, 
very much, usually, some times}. 


We now give the fuzzy linguistic graph whose vertices are 
Cy, Co, ..., C7 and the edges take values from L are as follows. 


C,[0 often 0 0 0 often 0 
C,}0 0 OO verymuch often 0 0 
M=C,}0 0  O- verymuch much 0 some 
C,}0 0 much 0 0 0 0 
C,)0 0 often 0 0 0 much 
C.j0 0 0 0 0 oO 0 
G00 0 0 much 0 0 
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Let X = (often, 0, some, 0, often, some, 0) is the fuzzy 
linguistic state vector. 
To study the effect of X on M using max {max {X, M}} 
after updating max {max {X, M}} is (often, often, some, very 
much, often, some, much) = X,. 


Now max {max {X,, M}} after updating is say X); 


X> = (often, very much, some, very much, often, some, very 
much). 


Now we find max {max (X2, M)} after updating we get say 
X3; 


X3 = (often, very much, some, very much, often, some, very 
much). 


Now suppose another expert wants to use for the same 
fuzzy linguistic state vector X and the same dynamical system, 
to study the effect using the ‘max min’ operation. 

Now max {min (X, M)} after updating we get say Xj, 

X, =(often, often, some, 0, often, some, often). 

We find max {min (X, M)} after updating we get say Xp. 
This is the way the operations are performed on the fuzzy 
linguistic matrix. 

X> = (often, often, some, 0, often, some often). 

For about Fuzzy linguistic relational map-model refer [53]. 

We can use any of the four types of operations {max, max} 
or {min, max} or {max, min} or {min, min} in these fuzzy 


linguistic matrices [53]. 


We will illustrate this situation by an example. 
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Example 3.3: Let us study the employee - employer fuzzy 
linguistic relational model. Suppose we have the following 
fuzzy linguistic concepts / attributes associated with the 
employee taken as the domain space 


D,; - Pay with allowances and bonus 


D, - Only pay to employee 


D; - Pay with allowance to employee 

D, - Best performance by the employee 

Ds; - Average performance by the employee 

De - Employee works for more number for hours. 


Suppose the following nodes / concepts are taken as the 
range space of the employer. 
R, - maximum profit to the employer 
R, - Only profit to the employer 
R; - Neither profit nor less to the employer. 
The fuzzy linguistic terms associated with the fuzzy 
linguistic domain and range spaces be taken as L. 


L = {0, gain, loss, no loss no gain, just gain, just loss, gain, 
heavy loss, good gain}. 
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good gain 


no loss\no gain 


y 


We give the associated fuzzy linguistic matrix N of the 
fuzzy linguistic graph. 


R, R, R, 
D,[ good gain 0 just loss | 
D, gain 0 0 
N= D,| no loss no gain 0 0 
D, gain 0 0 
D,; 0 nolossno gain 0 
D,|. just loss 0 0 | 
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Now we find the resultant of any fuzzy linguistic vector on 
N, the dynamical system associated with the problem. 


Let X = (gain, 0, loss, 0, gain, loss) be the given fuzzy 
linguistic state vector. The effect of X on the fuzzy linguistic 
dynamical system N is as follows: 


Max {max (X, N)} = ((good gain, gain, gain)} = Y 


We find max (max {Y, N‘}) = (gain, gain, loss, gain, gain, 
loss) = X, (say) after updating. 


We find max (max {X,, N}) 
= (good gain, gain, gain) = Y; = Y. 


This is the way the operations using linguistic matrices are 
performed on these models. 


Finally for the first time we introduce the notion of fuzzy 
linguistic relation equations and describe some of the properties 
related with them. 


Let L; and L, be any two fuzzy linguistic sets that is both L, 
and L, contain fuzzy linguistic terms or L; = L, otherwise. Let 
R be a fuzzy linguistic relation that is to each fuzzy linguistic 
term of L; two or more fuzzy linguistic terms in L, are 
assigned, that is to each x; € L, to the domain space (L, is the 
domain linguistic space) we associate a y2 € L, and here the 
degree of membership is not a value between [0, 1] but a fuzzy 
linguistic value from L. 


This we will first illustrate by some example. 
Let L; = {x1, Xo, X3, X4, Xs, Xo, X7} be seven fuzzy linguistic 


terms and L2 = {y1, y2, y3, Y4, Ys} be some five fuzzy linguistic 
terms with the following relation. 
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We have fuzzy linguistic membership matrix R associated 
with the above map. 


y Y2 Y3 V4 Ys 
x,| good fair 0 0 0 
x,| 0 0 worst 0 0 
x;| 0 fair 0 0 0 
x,| 0 0 0 bad good 
x;| best good 0 0 0 
X,| 0 0  v.bad 0 best 
x,| 0 0 0 far 0 


Thus for fuzzy linguistic sets L; and L2 we can have a fuzzy 
linguistic membership function R. The representation by the 
diagram will be known as the fuzzy linguistic sagittal diagram 
and the fuzzy linguistic matrix will be known as the fuzzy 
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linguistic membership matrix where the membership values are 
also fuzzy linguistic terms. 

We can compose two fuzzy linguistic binary relation 
provided both of them take their fuzzy linguistic membership 


values from the same fuzzy linguistic set L. 


We will just illustrate how such compositions are made. 


Suppose 
good bad _ fair best 
p= bad fair good good 
0 good fair good 
good bad good fair 
and 
good bad good fair 0 best 
Q= bad good best good bad 0 


best fair best 0 good bad 
0 fair. good bad good good 


be two fuzzy linguistic membership matrices of a fuzzy 
linguistic binary relation with the values taken from the same 
fuzzy linguistic space L. 


Thus we have used the innovative methods of applying 
matrices by constructing in the first place fuzzy linguistic 
matrices and secondly using them in fuzzy models like Fuzzy 
Cognitive Maps, Fuzzy Relational Maps and Fuzzy Relational 
Equations. 
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Now 

good 
Pod bad 

(@) = 

0 

good 
good _ bad 
bad = good 
best fair 
0 fair 


where R is found in the following way. 


L 
fee: .? 


bad 

fair 
good 
bad 


good 
best 
best 
good 


fair _ best 
good good 
fair good 
good fair 
fair 0 best 
good bad 0 
0 good bad 


bad good good 


tT; = max [min (good, good), min (bad, bad), 
min (fair, best), min (best, 0)] 
= max [good, bad, fair, 0] = good 


tf = max [min (good, bad), min (bad, good), 
min (fair, fair), min (best, fair)} 
= max {bad, bad, fair, fair} = fair 


r3 =max {min {good, good}, min {bad, best}, 
min {fair, best}, min {best, good} } 


= max {good, bad, fair, good} = good 
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t4 =max {min {good, fair}, min {bad, good}, 
min {fair, 0}, min {best, bad}} 
= max {fair, bad, 0, bad } = fair. 


and so on. 
We have the following fuzzy linguistic binary relation 


which is described by the fuzzy linguistic sagittal diagram and 
the related fuzzy linguistic membership matrix. 
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The fuzzy linguistic membership matrix associated with X is 


x 3 X4 Xs Xe X, 
x,/ medium fast slow 0 0 0 0 
Xs 0 0 fast just fast 0 0 0 
X; 0 slow 0 0 0 fast 0 
X4 0 0 O 0  veryslow 0 fast |. 
Xs 0 0 slow 0 0 0 0 
X¢ 0 0 O 0 fast 0 O 
x,| 0 0 O 0 0 slow fast 


Operations can be performed on these matrices. 


Chapter Four 


SUPERMATRICES AND THEIR 
APPLICATIONS 


The concept of super matrices was first introduced by [17]. 
Later the authors developed a new algebra using these matrices 
called super linear algebras [45]. Thus for the first time we use 
only super matrices to build super linear algebras. This 
application of super matrices is not only innovative and has lot 
of applications. Further one can find bulk eigen values using 
super square diagonal matrices. This method helps in easy 
comparison of eigen values and also bulk performances takes 
lesser time [45]. 


Apart from this super fuzzy matrices are used in the 
construction of super fuzzy models like New Super Fuzzy 
Relational model, New Super Fuzzy Cognitive Maps model, 
New Super Fuzzy Associative Maps model, we make use of 
them in the construction of super codes [50]. 


In this section for the first time we introduce the new notion 
of Super Fuzzy Relational Maps (SFRMs) models and they are 
applied to real world problems, which is suited for multi expert 
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problems. When we in the place of fuzzy relational matrix of 
the FRM (using single expert) use fuzzy supermatrix with multi 
experts we call the model as Super Fuzzy Relational Maps 
(SFRMs) model. 


We just recall the definition of Domain Super Fuzzy 
Relational Maps (DSFRMs) model. 


DEFINITION 4.1: Suppose we have some n experts working on a 
real world model and give their opinion. They all agree upon to 
work with the same domain space elements / attributes / 
concepts; using FRM model but do not concur on the attributes 
from the range space then we can use the special super fuzzy 
row vector to model the problem using Domain Super Fuzzy 
Relational Maps (DSFRMs) Model. 


The DSFRM matrix associated with this model will be given 
by Su 


tpveel, Ty oust, ty tp wt, 
D, 
D 
Sy= 
BR 
=| 8e | Sia Lael | 


where each S\, is am * t matrix associated with a FRM 
given by the i” expert having Dj, ..., Dm to be the domain 
attributes and (t,t, ... ti ) to be the range attributes of the i” 


expert, i = I, 2, ..., n and Sy the DSFRM matrix will be a 
special super row vector / matrix (1 Si Sn). 
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However if n is even a very large value using the mode of 
programming one can easily obtain the resultant vector or the 
super hidden pattern for any input supervector which is under 
investigation. 


These DSFRMs will be known as Domain constant 
DSFRMs for all the experts choose to work with the same 
domain space attributes only the range space attributes are 
varying denoted by DSFRM models. 


We see these super row vectors (super row matrices) find 
their applications in the DSFRMs model. 


Next we proceed on to recall the notion of super FRMs with 
constant range space attributes and varying domain space 
attributes. 


DEFINITION 4.2: Let some m experts give opinion on a real 
world problem who agree upon to make use of the same space 
of attributes / concepts from the range space using FRMs but 
want to use different concepts for the domain space then we 
make use of the newly constructed special super column vector 
as the matrix to construct this new model. 


The column attributes i.e. the range space remain constant 
as Rj, ..., Rs for all m experts only the row attributes for any i” 


expert is Di, D}, ...,D'; i = 1, 2, ..., m; vary from expert to 


expert. This system will be known as the Range constant fuzzy 
super FRM or shortly denoted as RSFRM model. 


Thus the associated special super column matrix S™ is 
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These super row vectors (super row matrices) find 
applications in the RSFRM models. 


DEFINITION 4.3: Suppose we have m experts who wish to work 
with different sets of both row and column attributes i.e. domain 
and range space using FRMs, then to accommodate or form a 
integrated matrix model to cater to this need. We make use of 
the super diagonal fuzzy matrix, to model such a problem. 
Suppose the first expert works with the domain attributes 


Deda D, and range attributes Rj, ..., Ri The second expert 


works with domain attributes D,, wile and with range 


attributes i ..., R2 and so on. Thus the m" expert works with 


ny 


Di", ....D" domain attributes and Ry",...,R" range 


attributes. We have the following diagonal fuzzy supermatrix to 
model the situation. We are under the assumption that all the 
attributes both from the domain space as well as the range 
space of the m experts are different. 
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The super fuzzy matrix S associated with this new model is 
given by 


| Ki RovcR RR eRe. ie RR RE 
Dif 7 
7} OM, (0) (0) 
D} 
D, 
: (0) (0) 
De 
(0) (0) 
Dr 
(0) (0) M, 
D” 


where each M, is a t; X n; matrix associated with the FRM, we 
see except, the diagonal strip all other entries are zero. 


We call this matrix as a special diagonal super fuzzy matrix 
and this model will be known as the Special Diagonal Super 
FRM Model which will be denoted briefly as (SDSFRM). 


The super diagonal matrices are used in the SDSFRM 
models. 


DEFINITION 4.4: Suppose one is interested in finding a model 
where some mn number of experts work on the problem and 
some have both domain and range attributes to be not 
coinciding with any other expert and a set of experts have only 
the domain attributes to be in common and all the range 
attributes are different. 
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Another set of experts are such that only the range 
attributes to be in common and all the domain attributes are 
different, and all of them wish to work with the FRM model 
only; then we model this problem using a super fuzzy matrix. 
We have mn experts working with the problem. 


Let the t; expert wish to work with domain attributes Fs 


Feige cae and range attributes q;, 3, - Tiny 
The t) expert works with P', P., ..., Pes as domain 
attributes and the range attributes q, ; qs. his Gans and so on. 


Thus for the t; expert works with Pe, oe Oe 


m(z,) 28 domain 
space attributes and q,, 43, +» Yy,)» 48 range attributes (I Si 


Sm(t)) andi Sn (t,)). 


So with these mn experts we have an associated super FRM 
matrix. Thus the supermatrix associated with the Super FRM 
(SFRM) model is a supermatrix of the form 


11 12 In 
Anceyncay Ancetynca) Ae cistinl 
21 22 2n 
Aie2yn(tt) Arit2)n(02) Antesnciey 
S(m) = 
A”! A”? Am™ 
m(tm)n(tl) m(tm)n(t2) m(um)n(tn) | 


where 
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J Jj Jj 
qi qh eae Gnit;) 
i 
Ff 
i 
A! = P, 
m(t;)n(t;) (a! ) 
. m(t; )n(t;) 
i 
Pug) 


1 Sis<mand 1 Sj Sn. S(m) is called the super dynamical FRM 
or a super dynamical system. 


, fil 
This matrix Anaya 


, corresponds to the FRM matrix of the 


(ij)"" expert with domain space attributes P', P}, ..., P\.,, and 
range space attributes qj, qj, ..., Gas , 1 <i <m and 
IlsSj sn. 


Super matrices are used as the dynamical system of SFRM 
models. 


Next we proceed onto show how super row vectors (super 
row matrices) are used in the construction of SDBAM models. 
We recall the definition of SDBAM model. 


DEFINITION 4.5: Suppose a set of n experts choose to work with 
a problem using a BAM model in which they all agree upon the 
same number of attributes from the space F,, which will form the 
rows of the dynamical system formed by this multi expert BAM. 
Now n distinct sets of attributes are given from the space F, 
which forms a super row vector and they form the columns of 
the BAM model. 


Suppose all the n experts agree to work with the same set of 
t-attributes say (x) X2 ... X,) which forms the rows of the synaptic 
connection matrix M. Suppose the first expert works with the p; 


set of attributes given by (y, y} ... y,, ), the second expert with 
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P2 set of attributes given by (y\ yy ... y, ) and so on. Let the i 
expert with p; set of attributes given by (yj, y, ... Vy ) for i = 1, 
2, ..., n. Thus the new BAM model will have its elements from F,, 
where any element in F, will be a super row vector, T = 
(WY; sis | yy yy cy, | Viv ---¥, J. Now the synaptic 
projection matrix associated with this new BAM model is a 
special row supervector M,. given by 


Lal 1 De sec 2 Nix n 

My Ve + Vp. Vi Vn 0 Vp, ee. Vy Var Vp: 
%, 
Mx, 
Xx, 


Here the elements /attributes from F, is a simple row rector 
where as the elements from F, is a super row vector. 


We call this model to be a multi expert Special Domain 
Supervector BAM (SDBAM) model and the associated matrix is 
a special row vector matrix denoted by M,. Let X = (x; X2 ... X)) 


n 


Bie VA 17 ol, 1 Vee | ee SE 
If X = (x; X2 ... xX) © F, is the state vector given by the expert we 


find 


XM, = YieF, 
We = Kies... 


and so on. This procedure is continued until an equilibrium is 
arrived. Similarly if the expert chooses to work with Y = 


LR Vecved | ee a - | WV, --¥;,] © Fy then we 
find the resultant by finding 
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YM,’ > X, then find XM, and proceed on till the system 
arrives at an equilibrium state. This model will serve the 
purpose when row vectors from F,, are a simple row vectors and 
row vectors from F, are super row vectors. 


In the following we describe how the super column vectors 
(super column matrices) we used in the construction of synaptic 
matrix of the SRBAM models. 


DEFINITION 4.6: Suppose we have a problem in which all m 
experts want to work using a BAM model. If they agree to work 
having the simple vectors from F, i.e., for the columns of the 
synaptic connection matrix i.e. there is no perpendicular 
partition of their related models matrix. 


The rows are partitioned horizontally in this synaptic 
connection matrix i.e., the m experts have distinct sets of 
attributes taken from the space F’, i.e. elements of F, are super 
row vectors. The resulting synaptic connection matrix M, is a 
special super column matrix. 


Let the I“ expert have the set of row attributes to be 
0G Re ox ee ), the 2" expert have the set of row attributes given 
by (xX; x3 ....x7,) and so on. Let the i" expert have the related 


row attributes as (x, x; 10%); i=1,2,..., m. 


Let the column vector given by all them is [y; ... Yn]. The 
related super synaptic connection matrix 
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Vi V2 e+ Vn 


M. is a special super column vector / matrix. 


Suppose an expert wishes to work with a super row vector X 


from F,, then X = [x) x} el, IE aes aed, [ sen] Hy Ay oe 


we find Xo M. > Y €F,, YM.'= X, € Fy we repeat the 
same procedure till the system attains its equilibrium i.e., a 
fixed point or a limit cycle. 


This model which performs using the dynamical system M. 
is defined as the Special Super Range BAM (SRBAM) model. 


We now proceed onto describe the SDSBAM model which 
makes use of the super diagonal matrices. We describe mainly 
these models to show how these models makes use of the super 
matrices. This SDSBAM model synaptic connection matrix is a 
diagonal super matrix. 
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DEFINITION 4.7: Suppose we have n experts to work on a 
specific problem and each expert wishes to work with a set of 
row and column attributes distinct from others using a BAM 
model. Then how to obtain a suitable integrated dynamical 


system using them. Let the first expert work with (x, x, ... x,,) 


attributes along the row of the related synaptic connection 
matrix of the related BAM and (yj y}... y, ) the attributes 
related to the column, let the second expert give the row 
attributes of the synaptic connection matrix of the BAM to be 
(x; x3 ...x° ) and that of the column be (y, y; ... y,,) and so 


Ny 


on. Let the i” expert give the row attributes of the synaptic 
connection matrix of the BAM to be (x x}... x’) and that of 


the column to be (yi ys iV) for i = 1, 2, ..., n, the 


supermatrix described by 


Vi Va Vp Vie see is. De eed a: 
ee ] 
A (0) (0) 

1 
x 
: (0) A (0) 
Mp = x 
(0) (0) (0) 
x} 
: (0) (0) 4, 
si L 4 
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where A! is the synaptic connection matrix using the BAM 


model of the ri expert, i = I, 2, ..., n where 
VY Vo oe Vy, 
oe 
ae | 
ApS 
x! 


(0) denotes the zero matrix. 


Thus this model has only non zero BAM synaptic connection 
matrices along the main diagonal described by Mp. The rest are 
zero. 


The dynamical system associated with this matrix Mp is 
defined to be the Special Diagonal Super BAM (SDSBAM) 
model. 


Yet another application (or use) of super matrices is in the 
super BAM models which is described below. 


DEFINITION 4.8: Suppose we have mn number of experts who 
are interested in working with a specific problem using a BAM 
model; a multi expert model which will work as a single 
dynamical system is given by the Super BAM (SBAM) 
model. Here a few experts have both the row and columns of 
the synaptic connection matrix of the BAM to be distinct. Some 
accept for same row attributes or vectors of the synaptic 
connection matrix but with different column attributes. Some 
accept for same column attributes of the synaptic connection 
matrix of the BAM model but with different row attributes to 
find the related supermatrix associated with the super BAM 
model. The supermatrix related with this new model will be 
denoted by M, which is described in the following: 
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Vee, ek. ee Ba 
Es 

A A! A 
X, 

my 

a a aL 
¥ 

x 

=| a aS A 
Xi, 


where A, is the synaptic connection matrix of an expert who 
chooses to work with (xi x}... a ) along the row of the BAM 
model and with (yj y} ... a ) along the column of the BAM 


model i.e., 


He ee 


I <$i<mand1 <j Sn. Thus for this model both the attributes 
from the spaces F,, and F,, are super row vectors given by 


2 
P2 


X= [x} x} side IB ge eer ol oe aa eee, 
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in F. and 


2 n 


PE Tee On Be ed, | De Ps 
Jrom the space or the neuronal field F,,. 


The supermatrix M, is called the synaptic connection 
supermatrix associated with the multi expert super BAM model 
(SBAM model). Now having defined the multi expert super BAM 
model we proceed on to describe the functioning of the super 
dynamical system. 

Let X = [x; x} sess, | a a a | set | Ay Re ee ef FF, 
be the super row vector given by the expert, its effect on the 
multi super dynamical system M,. 


Pet ic eee 


n n n 


= 1 ie ie, | Br Pell We Wt Sy 


Yo M! > X, €F,. 
X,0 M, + Y, eF;; 


and so on and this procedure is repeated until the system attains 
a equilibrium. 


Finally we just indicate how the super matrices are used in 
super fuzzy associative memories. For more about this refer 
[48]. 


A fuzzy set is a map nw : X — [0, 1] where X is any set 
called the domain and [0, 1] the range i.e., » is thought of as a 
membership function i.e., to every element x € X, p assigns a 
membership value in the interval [0, 1]. But very few try to 
visualize the geometry of fuzzy sets. It is not only of interest but 
is meaningful to see the geometry of fuzzy sets when we discuss 
fuzziness. Till date researchers over looked such visualization 
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[24], instead they have interpreted fuzzy sets as generalized 
indicator or membership functions; i.e., mappings p from 
domain X to range [0, 1]. But functions are hard to visualize. 
Fuzzy theorist often picture membership functions as two- 
dimensional graphs with the domain X represented as a one- 
dimensional axis. 


The geometry of fuzzy sets involves both domain X = 
(X),..., Xn) and the range [0, 1] of mappings up : X — [0, 1]. The 
geometry of fuzzy sets aids us when we describe fuzziness, 
define fuzzy concepts and prove fuzzy theorems. Visualizing 
this geometry may by itself provide the most powerful argument 
for fuzziness. 


An odd question reveals the geometry of fuzzy sets. What 
does the fuzzy power set F(2*), the set of all fuzzy subsets of X, 
look like? It looks like a cube, What does a fuzzy set look like? 
A fuzzy subsets equals the unit hyper cube I" = [0, 1]". The 
fuzzy set is a point in the cube I". Vertices of the cube I" define 
a non-fuzzy set. Now with in the unit hyper cube I" = [0, 1]" we 
are interested in a distance between points, which led to 
measures of size and fuzziness of a fuzzy set and more 
fundamentally to a measure. Thus within cube theory directly 
extends to the continuous case when the space X is a subset of 
RK: 


The next step is to consider mappings between fuzzy cubes. 
This level of abstraction provides a surprising and fruitful 
alternative to the prepositional and predicate calculus reasoning 
techniques used in artificial intelligence (AI) expert systems. It 
allows us to reason with sets instead of propositions. The fuzzy 
set framework is numerical and multidimensional. The AI 
framework is symbolic and is one dimensional with usually 
only bivalent expert rules or propositions allowed. Both 
frameworks can encode structured knowledge in linguistic form. 
But the fuzzy approach translates the structured knowledge into 
a flexible numerical framework and processes it in a manner 
that resembles neural network processing. The numerical 
framework also allows us to adaptively infer and modify fuzzy 
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systems perhaps with neural or statistical techniques directly 
from problem domain sample data. 


Between cube theory is fuzzy-systems theory. A fuzzy set 
defines a point in a cube. A fuzzy system defines a mapping 
between cubes. A fuzzy system S maps fuzzy sets to fuzzy sets. 
Thus a fuzzy system S is a transformation S: I" > I’. The n- 
dimensional unit hyper cube I" houses all the fuzzy subsets of 
the domain space or input universe of discourse X = {X, ..., 
Xn}. I? houses all the fuzzy subsets of the range space or output 
universe of discourse, Y = {yj, ..., yp}. X and Y can also denote 
subsets of R" and R?. Then the fuzzy power sets F (2*) and F 
(2*) replace I" and I’. 


In general a fuzzy system S maps families of fuzzy sets to 
families of fuzzy sets thus S: I" x...xI"" >I?! x...x I> Here 
too we can extend the definition of a fuzzy system to allow 
arbitrary products or arbitrary mathematical spaces to serve as 
the domain or range spaces of the fuzzy sets. We shall focus on 
fuzzy systems S: I" — I’ that map balls of fuzzy sets in I" to 
balls of fuzzy set in I”. These continuous fuzzy systems behave 
as associative memories. The map close inputs to close outputs. 
We shall refer to them as Fuzzy Associative Maps or FAMs. 


The simplest FAM encodes the FAM rule or association (Aj, 
B;), which associates the p-dimensional fuzzy set B; with the n- 
dimensional fuzzy set A;. These minimal FAMs essentially map 
one ball in I" to one ball in I’. They are comparable to simple 
neural networks. But we need not adaptively train the minimal 
FAMs. As discussed below, we can directly encode structured 
knowledge of the form, “If traffic is heavy in this direction then 
keep the stop light green longer” is a Hebbian-style FAM 
correlation matrix. In practice we sidestep this large numerical 
matrix with a virtual representation scheme. In the place of the 
matrix the user encodes the fuzzy set association (Heavy, 
longer) as a single linguistic entry in a FAM bank linguistic 
matrix. In general a FAM system F: I" > I° encodes the 
processes in parallel a FAM bank of m FAM rules (Aj, Bj), ..., 
(Am Bm). Each input A to the FAM system activates each stored 
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FAM tule to different degree. The minimal FAM that stores (Aj, 
Bj) maps input A to B;’ a partly activated version of B;. The 
more A resembles A;, the more B,’ resembles B;. The 
corresponding output fuzzy set B combines these partially 


activated fuzzy sets B),B},..., Bi. 


m 


B equals a weighted 


average of the partially activated sets B = w Bi +...+W,,B 


where w; reflects the credibility frequency or strength of fuzzy 
association (Aj, B;). In practice we usually defuzzify the output 
waveform B to a single numerical value y; in Y by computing 
the fuzzy centroid of B with respect to the output universe of 
discourse Y. 


More generally a FAM system encodes a bank of compound 
FAM rules that associate multiple output or consequent fuzzy 
sets B';, ..., BS with multiple input or antecedent fuzzy sets Ay. 
.... Aj. We can treat compound FAM rules as compound 
linguistic conditionals. This allows us to naturally and in many 
cases easily to obtain structural knowledge. We combine 
antecedent and consequent sets with logical conjunction, 
disjunction or negation. For instance, we could interpret the 
compound association (A', A’, B), linguistically as the 
compound conditional “IF X' is A' AND X? is A’, THEN Y is 
B” if the comma is the fuzzy association (A', A’, B) denotes 
conjunction instead of say disjunction. 


We specify in advance the numerical universe of discourse 
for fuzzy variables X', X° and Y. For each universe of discourse 
or fuzzy variable X, we specify an appropriate library of fuzzy 
set values A’, ..., A,’ Contiguous fuzzy sets in a library 
overlap. In principle a neural network can estimate these 
libraries of fuzzy sets. In practice this is usually unnecessary. 
The library sets represent a weighted though overlapping 
quantization of the input space X. They represent the fuzzy set 
values assumed by a fuzzy variable. A different library of fuzzy 
sets similarly quantizes the output space Y. Once we define the 
library of fuzzy sets we construct the FAM by choosing 
appropriate combinations of input and output fuzzy sets 
Adaptive techniques can make, assist or modify these choices. 
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An Adaptive FAM (AFAM) is a time varying FAM system. 
System parameters gradually change as the FAM system 
samples and processes data. Here we discuss how natural 
network algorithms can adaptively infer FAM rules from 
training data. In principle, learning can modify other FAM 
system components, such as the libraries of fuzzy sets or the 
FAM-rule weights wj. 


We now define the notion of super fuzzy associative 
memories of different types. 


DEFINITION 4.9: We have a problem P on which n experts 
wishes to work using a FAM model which can work as a single 
unit multi expert system. Suppose all the n-experts agree to 
work with the same set of attributes from the domain space and 
they want to work with different and distinct sets of attributes 
from the range space. Suppose all the n experts wish to work 
with the domain attributes (x; x2 ... x, from the cube 
I' =[0, I]x...x[0, 1]. Let the first expert work with the range 
eee 
t—times 


attributes (y, y, ... y, ) and the second expert works with the 
range attributes (y; y; ... Ye, ) and so on. Thus the i'" expert 


works with the range attributes (y y, hs y, ) t= 1,2, 0, N. 
Thus the range attributes 
VED Vy Tie en el ee 
are taken from the cube "*"?» =[0, 1]x[0, I]x...x[0, /]. 
Dit Dy+..+P, times 


This we see the range attributes are super row fuzzy vectors. 


Now the matrix which can serve as the dynamical systems 
for this FAM model is given by Fr. 
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VV Wee Dy,. SB Way, 
x, 
A A, A 


n 


Clearly Fp is a special super row fuzzy vector. Thus F: I' UO 


[P*P2*-"Pa Suppose using an experts opinion we have a fit 


vector, A = (a), Q, ..., a) ; a; €{0,1}, then A o Fr = max min 
(ai, fi); a; € A and fi; © Fr. Let A o Fr = B = (bj), then Fr 0 B = 
max min (fj, b') and so on, till we arrive at a fixed point or a 
limit cycle. The resultant fit vectors give the solution. This Fr 
gives the dynamical system of the new model which we call as 
the Fuzzy Special Super Row vector FAM model (SRFAM 
model). 


DEFINITION 4.10: Suppose we have n experts working on a 
problem and they agree upon to work with the same range 
attributes and wish to work with distinct domain attributes 
using a FAM model. We built a new FAM model called the 
special super column fuzzy vector FRM model (SCFAM) and its 
related matrix is denoted by F,. The fit vectors of the domain 
space are simple fit vectors where as the fit vectors of the range 
Space are super row fit vectors. 


Now we describe the special column fuzzy vector FAM, F. 
by the following matrix. The column attributes of the super fuzzy 
dynamical system Fare given by 


(V1.2. vy) € L =[0, I]x...x[0, 7]. 
~~" 


s—times 
The row attributes of the first expert is given b 
P & a 
x!,x),...,x. ), the row attributes of the second expert is given 
19% D; DP &§ 
by (x, x3, sede ). Thus the row attributes of the i” expert is 


given by (x; x Le a ), t= 1, 2, .., Nn. 
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We have 
VY, Vn oe Vs 
1 ay 
x 
3 A 
1 
Xp, 
2 
x 
A, 
_ 2 
F.= Xp 
n 
x, 
0 
x” 


to be a special super column fuzzy vector / matrix, where 


Vy V2 vee Vs 


i=1/, 2, ..., nis a fuzzy p; x s matrix. Suppose the expert wishes 


to work with a fit vector (say) [x} x, ... x [Se 26 ee a eer | 


xi X;...x), J. Then X o F.= B where B is a simple row vector 
we find F,o B=[0, 1] x...x[0, I]=I 
sf 
P,t+p2t...+p, times 


to work till we arrive at an equilibrium state of the system. 


; we proceed on 


Pite--+Py 
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Next we proceed on to define FAM when n expert give 
opinion having distinct set of domain attributes and distinct set 
of range attributes. 


DEFINITION 4.11: Let n experts give opinion on a problem P 
and wish to use a FAM model, to put this data as an integrated 
multi expert system. 


Let the first expert give his/her attributes along the column 
as (Yi Vy... Yn ) and those attributes along the row as 


1 Jl 1 
(X, Xz ...X, ). 


| a ee ee ee, and (x, KS 263 * 0) be the column and 


row attributes respectively given by the second expert and so 
on. 


Thus any i" expert gives the row and column attributes as 
pte eee x) and (yj yi, ... y,) respectively, i =1, 2, 3, ..., n. 


So for any i” expert the associated matrix of the FAM would 
be denoted by A’ where 


Now form the multi expert FAM model using these n FAM 
matrices A’, A’, ..., A" and get the multi expert system which is 
denoted by 
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VII, WII. ae ViVi, 
x [ 
Al (0) (0) 
Xp, 
x 
(0) (0) 
Fo= xt 
(0) (0) 
pa 
: (0) (0) A’ 


This fuzzy supermatrix Fp will be known as the diagonal 
fuzzy supermatrix of the FAM and the multi expert system which 
makes use of this diagonal fuzzy supermatrix F'p will be known 
as the Fuzzy Super Diagonal FAM (SDFAM) model. Now the 
related fit fuzzy supervectors of this model F, and Fy are fuzzy 


. 1 1 1 2 22 2 
super row vectors given by X= (X, xX, ...X, |X, Xz --. Xp, | 


Xi X)...x, ) € Fy and Y= (VW. V, a Eee i Vee | me 
Vie Ve wv Vg) Cy 


Now this new FAM model functions in the following way. 


Suppose the expert wishes to work with the fuzzy super state 


fit vector X = (x, x) ne en eee | 


n n n 
vee Ky | vee | Mp Hy oe By 


gt I= (Pi De Ia NW Da Ig DI) EP 
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Now Fpo Y = X; € F, and X;0 Fp = Y;€ F, and so on. 


This procedure is repeated until the system equilibrium is 
reached. 


DEFINITION 4.12: Let us suppose we have a problem for which 
mn experts want to give their opinion. Here some experts give 
distinct opinion both for the row attributes and column 
attributes. Some experts concur on the row attributes but give 
different column attributes and a few others have the same set 
of row attributes but have a different set of column attributes. 
All of them concur to work using the FAM model. To find a 
multi expert FAM model which can tackle and give solution to 
the problem simultaneously. 


To this end we make use of a super fuzzy matrix F, which is 
described in the following. Let the mn experts give their domain 
and column attributes as follows. 


The first expert works with the domain attributes as 
(4.35 ee) and column attributes as (y, y, seg Ve Dh The 


second expert works with the same domain attributes viz 
(ope eee x) and but column attributes as (y, y, ... yi). The 
i” expert, 1S i <n works with (x} x)... xy ) as the domain 


attributes and (y; V% ... y, ) as the column attributes. 


The (n + 1)" experts works with the new set of domain 
attributes (Xx; X; ... x ) but with the same set of column 


attributes viz. (y| y, ... Vad: 


th . . De beats 2 
Now the n + j” expert works with using (x; X; ...X,,) as 
the domain attribute and (y\ y%, ... Y;, ) as the column attribute 
I <j <n. The (2n + 1)" expert works with (te he ees Ga) as the 


row attribute and (y; yy ... yi, ) as the column attribute. 
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Thus any (2n + k)" expert uses ( a eS as a ) to be the row 
attribute and (y\ y) ... v) to be the column attribute Is k S 
n. Thus any (t, + 1)" expert works with (xi x‘... xi, ) as the row 


attribute (I St Sm) and (y, y, ... y,,) as the column attribute 


lsrsn. 


Now as 1 St S<mand 1 Sr Sn we get the FAM matrices of 
all the mn experts which is given by the supermatrix F,. 


Vie Ve Ge Ve eV. is DE VE we 
x; [ 
| 4 4, 4, 
x», 
a 
A A& 
Ps 7 a 
x 
A” A Ar 
x, L 


i; ° ° 5 th 
where Aj, is a fuzzy matrix associated with ij” expert 
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yw 


1 <i Smand 1sj Sn. This model is known as the multi expert 
fuzzy Super FAM (SFAM) model. The fit vectors associated with 
them are super row vectors from F, and Fy. 


The fit super row vector X from F,, is 


X=[x} x, aa, [ae ee ara ass ua ae ee | 


"P2 


and the 


Xe F,= [?**? =[0, 1]x[0, 1]x...x[0, Z]. 
—_— 


Pi +P2+...+ Py times 


The fit super row vector Y from F, is 


n n 


VSO Ib se De || YE Me soe Ve, | Sen | D8 se Ht 


Y eF, _ [Uta t +40 =[0, 1|x[0, Ix... [0, 7]. 
OOF 


Gi +G2+-..+Qm times 


Thus if 


= 1 1 1 2 2. 2, m m m 
X= [Ky Hy vee X pq | My Ay vee Xp, | ove | AY My --- HE 
is the fit vector given by an expert; its effect on F, is given 
by Xo F, = Y € F,; now F,o Y= X' e€ F, then find the effect of 
Monk eX ora € F,, and so on. 
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We repeat this procedure until we arrive at a equilibrium 
state of the system. 


For exact use of the super matrices please refer [17]. We 
have shown how the abstract model is associated with a super 
matrix. 


We have used super matrices to construct super linear 
algebra [45]. We have developed super linear algebras using 
super matrices is an analogous way like linear algebra and 
derived several classical theorems. Here we recall a few of 
them for more about the innovative use of super matrices refer 
[45]. 


DEFINITION 4.13: Let V = (V|| ... |V,,) be a super vector space 
and T, = (T; | ... | T,) be a linear operator in V. IfW = (W, | ... | 
W,,) be a super subspace of V; we say that W = (W, | ... | W,,) is 
super invariant under T if for each super vector a@ = (a | ... 
|a,) in W = (W, | ... | W,,) the super vector T,(a) is in W = (W, | 
.. | W,) i.e. if TW) is contained in W. When the super subspace 
W = (W, | ... | W,) is super invariant under the operator T, = 
(T, | ... | T,) then T, induces a linear operator (T,)y on the super 
subspace W = (W,| ... | W,,). 


The linear operator (T,) w is defined by (T,)w (a) = T;(a) for 
ainW=(W, |... | W,) but (T,)w is a different object from T, = 
(T; | ... | T,) since its domain is W not V. 


When V =(V,|... | Vn) is finite (nj, ..., n,) dimensional, the 
invariance of W = (W, |... | Wn) under T, = (T, | ... | Ty) has a 
simple super matrix interpretation and perhaps we should 
mention it at this point. Suppose we choose an ordered basis B 


=4Bi |: | Ba) (Gy 20, (2/07 22a, tor V = WV] ee 
V,,) such that B’ = (a; ... a, |...| Qj... @f) is an ordered basis 
for W = (W, | ... | Wy); super dim W = (mn, ..., 1). Let 


A =[T,], so that 
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n n, 
Ta; = Sa, au; | ae | Ssas.at 
i,=l 


Since W = (W, | ... | W,) is super invariant under T, = (T, | 
... | Tn) and the vector T,a,; = (Ta, |...| T,a; ) belongs to W = 
(W, |... | Wn) for j, <1. This means that 


i=l i, =1 


i th 
ae Sau | Sate 


je<ty t=, 2, ..., n. In other words Aj, = (Ai. 


Uh 


[..[A2, = © 


|... | 0) ifj, <r, and 1, > rn, 


1 C¢, 
0 |o}| 0 
0 D, 
> |B & ' 
= Oo: “DB 
0 
0 0 


where B, is an r X r, matrix, C; is ar, X (nm, — r,) matrix and D, is 
an (n, — 1) X (nm — 1) matrix t= 1, 2, ..., n. 


In view of this we prove the following interesting lemma. 
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LEMMA 4.1: Let W = (W, | ... | W,) be an invariant super 
subspace for T, = (T; | ... | Ty). The characteristic super 
polynomial for the restriction operator 


(Ty = (Ty, |---|.) divides the characteristic super 


polynomial for T,. The minimal super polynomial for 
(Z,),.= (Q)y, | ---\(Z.),, ) divides the minimal super polynomial 


for T;. 
Proof: We have [T;]3 = A where B = {B; ... B,} is a super basis 
for V =(V,|... | Vn); with B; = {on he a, } a basis for Vj, this 


is true for each i, i = 1, 2, ..., n. A is a super diagonal square 
matrix of the form 


where each 
B.C. 
Aj= 
0 D, 
fori=1, 2, ...,n; Le. 
B, C, 
0 0 0 
0 OD, 
B 
0 : 0 
A= 0 
0 0 
n C, 
0 0 0 
0 OD, 
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and B=((T,),,|, where B’ is a basis for the super vector 


subspace W = (W;, | ... | W,) and B is a super diagonal square 
matrix; 1.e. 


Now using the block form of the super diagonal square matrix 
we have super det (xI — A) = super det (xI — B) x super 
det (xI — D) 


i.e. (det (xl; — Aj) |...| det (xI, -A,)) 
= (det (xI; —B,)det(xI/ —D,)|...| 
det (xI/, -B,) det(xI, —D,)). 


This proves the restriction operator (T,)w super divides the 
characteristic super polynomial for T,. The minimal super 
polynomial for (T;)w super divides the minimal super 
polynomial for T,. 


It is pertinent to observe that Ij , I/,1,,...,1, represents 


different identities i.e. of different order. 
The K" row of A has the form 
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where Ce is some r; X (nm — r,) matrix; true for t = 1, 2, ..., n. 


Thus any super polynomial which super annihilates A also 
super annihilates D. 


Thus our claim made earlier that, the minimal super 
polynomial for B super divides the minimal super polynomial 
for A is established. 


Thus we say a super subspace W = (W, | ... | W,) of the 
super vector space V = (V, | ... | Vn) is super invariant under 
T, =(T, |... | Ta) if TW) C W ie. each T\(W;) C Wj; for i = 1, 
2,....n Le. if a =(a)|... | Gn) e W 


then T,a = (T, a, | ... | TxOtn) where 
St pole tl ici © 
Oy = Kj, Fie +X, 
2 2 2 1 
Oy =X, OA) +...+X, OL, 


and so on 


=) nn nn 
O,, =X, 0, +..,+%, a, . 


ad 114 bot A nonin non n 
Ta = (t, x,0,+...+1,x,0, |-..[t xpa;+...+t, x, a, ). 


Now B described in the above theorem is a super diagonal 
matrix given by 
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Thus the characteristic super polynomial of B i.e. (T;)w) is 


g=(gi|---| Sn) = 
((x—cl)* gOC= ee: yx |....(x—¢; )* (KCK re) 


where e; = dim W, fori=1,2,...,K,andt=1,2,...,n. 


Now we proceed onto define T, super conductor of any a 
into W =(W|... | Wa). 


DEFINITION 4.14; Let V = (V;| ... | V;,) be a super vector space 
over the field F. W = (W, | ... | W,) be an invariant super 
subspace of V for the linear operator T, = (T; | ... | T,) of V. Let 
a=(a;|... | &) be asuper vector in V. The T-super conductor 


of « into W is the set S,(a@,W) =(S,(a@3W,) |...|S; (@,3W,)) 


which consist of all super polynomials g = (g; | ... | Zn) (over 
the scalar field F) such that g(T,)a is in W, i.e. (g)(T))@ | ... | 
£n(T) &)) €W=(W,| ... | W,). i.e. g(T) a) © W; for every i. Or 
we can equivalently define the T, — super conductor of a in W is 
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a T; conductor of a; in W; for every i = 1, 2, ..., n. Without loss 
in meaning we can for convenience drop T, and write the super 
conductor of a into W as S(a;W) = (S(a,;W,)|...|S(@,;W,)) - 

The collection of polynomials will be defined as super 
stuffer this implies that the super conductor, the simple super 
operator g(T,) = (gi(T)) | ... | Sn(T,)) leads the super vector a 
into W. In the special case W = (0 | ... | 0), the super conductor 
is called the T, super annihilator of a. 


The following important and interesting theorem is proved. 


THEOREM 4.1: Let V = (V;| ... | V,) be a finite dimensional 
super vector space over the field F and let T, be a linear 
operator on V. Then T, is super diagonalizable if and only if the 
minimal super polynomial for T, has the form 


P= Or) <P) =). OH NEHGS ee) 


1 1 
where  (C, ...Cg |---|€f -- CK are such that each set 


a steer are distinct elements of F for t = 1, 2, ..., n. 


Proof: We have noted that if T, is super diagonalizable, its 
minimal super polynomial is a product of distinct linear factors. 
To prove the converse let W = (W, | ... | W,) be the super 
subspace spanned by all of the characteristic super vectors of T, 
and suppose W = (W, |... | Wn) # (Vi | ... | Wn) ie. each W; # 
V;. By the earlier results proved there is a super vector @ not in 
W = (W, | ... | Wy) and a characteristic super value 
C, =(c; ,...C; ) of T; such that the super vector B = (T—c,loa 


ie. (B, |...[B,) =(T,-¢;,1,)a,|...\(1, -ef1,)a,) lies in W = 
(W, | ... | W,). Since (B; | ... | Ba) is in W, 


B= (4... +B, [Bites Bg [BT ++ BR) 
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where B, = Bi tsnce By, for t= 1, 2, ..., n with T.B; = c;B;; 1 <i< 
K ie. (TB; |..-T,B)) = (c, B;,|---le2 Bi); (lL < i: S K,) and 


therefore the super vector 


h(T, )B = (hy (cB; +... +h; (Cy, )Bx, |. 
h, (c;' )B; +...+b, (CK, )Bx, ) 
= (h, (7, )B,|--.1h, (TB) 


is in W =(W|]|... | W,) for every super polynomial h = (h | ... | 
h,). 


Now (x — c;) q for some super polynomial q, where p = (p, | 
... | Pn) and q = (qi |... | qn). 


Thus p = (x — cj) q implies 
Pp = (pil... | Pa) 
= (&-c; )al-..—e} )a,) 


ie. (q, —4,(¢; )I-..1d,—d, Ce; = (x=; )hy|...(—eF Yh,). 
We have 


q(T, a—q(c,)a=(q,(T, a, —4,(c; Jo |...| 

G(T, On — In (Cj 0) 

= h(T,)(T,-¢,Da=h(T,)p 

= (6,(0)C,—¢,1,)o, |-[h, (1G, <4 1) @,) 


=(h,(T))B, |...|h, (T, )B,)- 
But h(T,)B is in W =(W,| ... | Wa) and since 


0 = p(T, a = (p, (7, a, |...[p, (T,, Ja, ) 
= (T, -¢,Dq Ta 


= (7, =e) q, (T,)a,|...(Z, —¢; I,)q, (7, on); 
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the vector q(T,)a is in W. Therefore q(cj)a is in W. Since a is 
not in W we have q(c;) =(q,(c;, )I---14, (Cc; ))=(0l...|0). Thus 
contradicts the fact that p = (pi |... | Pn) has distinct roots. 

If T, is represented by a super diagonal square matrix A in 
some super basis and we wish to know if T, is super 


diagonalizable. We compute the characteristic super polynomial 
f=(f; |... | f,). If we can factor 


f 


CE | 0: "Ta 


1 a n 
= ((x-€})" Ke |e?) (x8) 


we have two different methods for determining whether or not T 
is super diagonalizable. One method is to see whether for each 1 


= 1, 2, ..., n we can find di independent characteristic super 
vectors associated with the characteristic super values c; . The 


other method is to check whether or not 


(T, =e). =eDie. (=e). -e) 
|..-1(E, - ep I, ..-(T, — ex, 1) 


is the super zero operator. 


Several other interesting results in this direction can be 
derived. Now we proceed onto define the new notion of super 
independent subsuper spaces of a super vector space V. 


DEFINITION 4.15: Let V = (V;| ... | V,) be a super vector space 
over F. Let W, =(W) |...|W,"), W, =(W2 |...|We')... 
W, =(W,|....W2) be K super subspaces of V. We say W,, ..., 
Wx are super independent ifa; + ... + Ax = 0; a © W; implies 
each a;= 0. 


a, = (a |...|a')eW,=(W} |...|W"); 
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true for i = 1, 2, ..., K. If W; and W> are any two super vector 
subspaces of V = (V; | ... | Vx), we say W, =(W;' |...|W,") and 
W,=(W;|...|W;') are super independent if and only if 
WOW, = WOW. IW OWE) = (0| 0 |... 10. If Wr Wo 
., Wx are K super subspaces of V we say Wi, W>, ..., We are 
independent if W, 1 W227... OA We = WLOAW, O...We|..| 
Wl AWS O...0We) = (0 | ... | 0) . The importance of super 
independence in super subspaces is mentioned below. Let 


Wo = Wht. + W's 
= (Wi+...4We || +..4W2) 
= WW |..1W,) 
W,' is a subspace V; and W/=W/'+...+W, true for i = 1, 2, ..., 


I 


n. Each super vector a in W can be expressed as a sum 
a=(al|...,a')=((a)+...+a)|...a7t...¢a%)) ie. each 
a=a+...ta,; ad © W, If Wi, Ws ..., We are super 
independent, then that expression for a is unique; for if 


a=(B,+...+ Be) = (Bi t+...+ By |B +...4+ 82) 
B<W,;i=1,2,..,K. B = Bo +...+ B then 


a-a=(0]...|0) = (a -B))+...+ (ay - BI. 
(a) —B')+...+( Be — a) 


hence each a; - 8; =0; 1 Si SK; t=1, 2, ..., n. Thus W,, W2, 
., Wr are super independent so we can operate with super 
vectors in W as K-tuples (a shiv Vitek (Oh ect )s ai eW,; 1 
SiskK;t=1, 2, ..., n. in the same way we operate with R* as 
K-tuples of real numbers. 


LEMMA 4.2: Let V = (V; | ... | V,) be a finite (nj, ..., My) 
dimensional super vector space. Let W), ..., Wx be super 
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subspaces of V and let W =(W;\'+...W,|...|W" +...+WZ). The 
following are equivalent 


(a) Wi, ..., Wx are super independent. 

(b) For each j; 2 sj SK, we have W, 0 (W, + ... + Wj-) = 
{(0|...| 0} 

(c) If B; is a super basis of W;, 1 <i SK, then the sequence 
B=(B, ... Bx) is a super basis for W. 


The proof is left as an exercise for the reader. In any or all of the 
conditions of the above stated lemma is true then the supersum 


W=W, +... + We = (W, +...+ Wy |...)W +...+W2) where 
W, =(W! |...|W,") is super direct or that W is a super direct 


sum of W,, .... We ie. W = W, @ ... © We ie. 
(W; ©...8W, |...| W ®...@ Wz) . If each of the W; is (1, ..., 
1) dimensional then W = W, @® ... © Wy, 


=(W; ®...8W!|...|W" ®...8 W"). 


DEFINITION 4.16: Let V = (V;| ... | V,) be a super vector space 
over the field F; a super projection of V is a linear operator E, 
on V such that E?=E, ie. E, = (E; | ... | E,) then 


E? =(E)|...|E{) = (E; | ... | Ey) ie. each E;is a projection on 
V;;i = 1, 2, ..., n.. Suppose E, is a projection on V and R = (R; | 
... | R,) is the super range of E, and N = (N_ | ... | N,) the super 
null space or null super space of E;. The super vector B = (f; | 
... | B) is in the super range R = (R; | ... | R,) if and only if E,B 
= B ie. ifand only if (E)f, | .... | Enfi) = (81 \ ... | B) Le. each 
Ef; = 8 fori =1, 2, ..,n.. If B= Ea ie. B=(B; |... |B) = 
(E; a | ... | EnQ,) where the super vector a = (a, | ... | Q,) then 
EB = E?a = E,a= p. Conversely if B= (i | | B) = EB 
=(E/f, | .. | Enf,) then B= (B, | ... | B) is in the super 
range of Es. Thus V=R @N i.e. V = (V;| .... | Vi) = (Ri ON, | 
won| ge Ng 


Further the unique expression for = (a) | ... | Q,) as asum 
of super vectors in R and N is a = E,a + (a— E;@) i.e. @ = 
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E,a@, + (a; -— E;a,) for i = 1, 2, ..., n. From what we have stated it 
easily follows that if R and N are super subspace of V such that 
V=RQ@N i.e. V=(V;| ... | Vi) = (Ri ON] | ... | Rn ON, then 
there is one and only one super projection operator E, which 
has super range R = (R; | ... | R,) and null super space N = (N; 
| ... | N,). That operator is called the super projection on R 
along N. 


Any super projection E, is super diagonalizable. If 
{(@4.G, |...|@'...a@) is a super basis for R = (R; | ... | Ry) 
1 1 n 
teem [ol Oe et 


| N,) then the basis B = (a)... G,, )|--sepne, ) = (B, |... | B,) 


super diagonalizes E, . 


and (a ...@” ) is a super basis for N = (N_ | ... 


(E,), = 
= (B,J)... Ey.) 
where I, is ar, * r,; identity matrix; t = I, 2, ..., n. Thus super 
projections can be used to describe super direct sum 
decompositions of the super vector space V = (V;| ... | Vn). 


THEOREM 4.2: Let F be the field of real numbers or the field of 
complex numbers. Let A be a super diagonal matrix of the form 
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A, 
0 
A= 
0 
be a(n; X nj, ..., Nn X Ny) matrix over F. The super function g = 


(g1 | ... | Sn) defined by g(X, Y) = Y'AX is a positive superform 
ie | F"”") if and only if there exists 


an invertible super diagonal matrix 


on the super space (F"”" 


Each P; is an; xn; matrix i = 1, 2, ..., n with entries from F 
such thatA =P P ;i.e., 
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DEFINITION 4.17: Let 


be a superdiagonal matrix with each A; an; * n; matrix over the 
field F; i =1, 2, ..., n. The principal super minor of A or super 
principal minors of A (both mean the same) are scalars 


A,(A) = (A (A) - [Ar (4) 
defined by 
Ait F P 
yee 
Are: Ay 
A,(A) = superdet 0 : 0 
AP Ae 
Al acd 
0 0 : : 
Ak Aga 
Ai, Ai, Al, Ai, 
—|det| : ai det | : : 
Ay Au 4, Ay, 


for 1 Sk, Sn, and t = 1, 2, ..., n. 
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Several other interesting properties can also be derived for 
these superdiagonal matrices. 


We give the following interesting theorem and the proof is 
left for the reader. 


THEOREM 4.3: Let f = (f; | ... | fi) be a superform on a finite 
(n1, ..., Nx) dimensional supervector space V = (V; | ... | V,) and 
let A be a super diagonal matrix of fin an ordered superbasis B 
= (B,| ... | B,). Then fis a positive superform if and only if A = 
A* and the principal super minor of A are all positive. 

bes 


Note: The principal minor of (A; | ... | An) is called as the 
principal superminors of A or with default of notation the 
principal minors of {Aj, ..., An} is called the principal super 
minors of A. 


T, = (T, | ... | T,) a linear operator on a finite (nj), ..., Dn) 
dimensional super inner product space V = (V | ... | Vn) 1s said 
to be super non-negative if T, = T. 


1.€., 
(T1|... | T= (FG 


a 


n 


ie. T, = T fori=1,2,...,n and 
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(Ts & | &) = (Ti oy | 1) |... | Tn On | On) 2 (0 ... 10) 
for all a = (a; |... | Qn) in V. 
A super positive linear operator is one such that T, = T. and 
(To | o) = (Ty | 01) |... | (Tn On | On)) > (0 |... 10) 
for all a = (a; |... | Gn) #(O| ... |0). 

Several properties enjoyed by positive operators and non 
negative operators will also be enjoyed by the super positive 
operators and super non negative operators on super vector 
spaces, with pertinent and appropriate modification. 


Throughout the related matrix for these super operators T, will 
always be a super diagonal matrix A of the form 


. . . * 
where each A; is a n; X n; square matrix, | <i<n, A=A and 
the principal minors of each A; are positive; 1 <i<n. 


Now we just mention one more property about the super 
forms. 


THEOREM 4.4: Let f = (f; | ... | f,) be a super form on a real or 
complex super vector space V = (V; | .. | V,) and 


n 


ia ane \ a super basis for the finite dimensional 


super subvector space W = (W, | ... | W,) of V = (V; | ... | Vi). 
Let M be the super square diagonal matrix where each M; in M; 
is ar; X vr; super matrix with entries (I < i < n). 
Mi, = F(Q, @), Le. 


1 1 
{ay se, 
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and W'= (W/|...|W,) be the set of all super vectors B = (B; | 
.. |B) of Vand WI\W'=(W, WI |...|WNW =|... | 
0) if and only if 


is invertible. When this is the case, V=W+W'i.e.V=(V,| ... 
|V =, +h! |... [+0 ) 


The proof can be obtained as a matter of routine. 


The projection E, = (E; | ... | E,) constructed in the proof 
may be characterized as follows: 


E.B = a; 
(E:B, |... | ExBn) = (1 |... | On) 


is in W and B — a belongs W'= (W/ | ... | W{ ). Thus E, is 
independent of the super basis of W = (W, | ... | W,) that was 
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used in this construction. Hence we may refer to E, as the super 
projection of V on W that is determined by the direct sum 
decomposition. 


V=Wew; 
(Vil...1Viy= (W,eW, |... |W, OW, ). 


Note that E, is a super orthogonal projection if and only if 
W'= W* =(W,'| ... |W). Now we proceed onto develop the 


analogous of spectral theorem which we call as super spectral 
theorem. 


Finally these super matrices are used to develop the notion 
of Leontief economic super models. 


Matrix theory has been very successful in describing the 
interrelations between prices, outputs and demands in an 
economic model. Here we just discuss some simple models 
based on the ideals of the Nobel-laureate Massily Leontief. Two 
types of models discussed are the closed or input-output model 
and the open or production model each of which assumes some 
economic parameter which describe the inter relations between 
the industries in the economy under considerations. Using 
matrix theory we evaluate certain parameters. 


The basic equations of the input-output model are the 
following: 


a, AyD ain P; P; 
4, An» Bon |) | Po Wi] Pe 
ayy ano aie" Aan P, P, 


each column sum of the coefficient matrix is one 


i. = pj20,1=1,2,...,n. 
il. a Og = ly 2g te 
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forj=1,2,...,n. 


P, 
are the price vector. A = (aij) is called the input-output matrix 
Ap = p that is, (I— A) p=0. 


Thus A is an exchange matrix, then Ap = p always has a 
nontrivial solution p whose entries are nonnegative. Let A be an 
exchange matrix such that for some positive integer m, all of the 
entries of A™ are positive. Then there is exactly only one 
linearly independent solution of (I — A) p = 0 and it may be 
chosen such that all of its entries are positive in Leontief open 
production model. 


In contrast with the closed model in which the outputs of k 
industries are distributed only among themselves, the open 
model attempts to satisfy an outside demand for the outputs. 
Portions of these outputs may still be distributed among the 
industries themselves to keep them operating, but there is to be 
some excess some net production with which to satisfy the 
outside demand. In some closed model, the outputs of the 
industries were fixed and our objective was to determine the 
prices for these outputs so that the equilibrium condition that 
expenditures equal incomes was satisfied. 


x; = monetary value of the total output of the i" industry. 


dj = monetary value of the output of the i" industry needed to 
satisfy the outside demand. 
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6i; = monetary value of the output of the i" industry needed by 
the j'" industry to produce one unit of monetary value of its own 
output. 


With these qualities we define the production vector. 


xX, 
ma 
x= ; 
X 
the demand vector 
d, 
d 
2 
d=]. 
d, 
and the consumption matrix, 
91, Op Oi, 
0o,, O Oo 
21 22 2k 
C= : 
Ox, Ox Orn 


By their nature we have 
x>0,d>0andC20. 


From the definition of oj and x; it can be seen that the quantity 
Oj X1 + O72 X2 +...4+ Oix Xk 


is the value of the output of the i" industry needed by all k 
industries to produce a total output specified by the production 
vector x. 
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Since this quantity is simply the i” entry of the column 


vector Cx, we can further say that the i" entry of the column 
vector x — Cx is the value of the excess output of the i" industry 
available to satisfy the outside demand. The value of the outside 
demand for the output of the i" industry is the i” entry of the 
demand vector d; consequently; we are led to the following 
equation: 


x—Cx=dor 
(I-C)x=d 


for the demand to be exactly met without any surpluses or 
shortages. Thus, given C and d, our objective is to find a 
production vector x = 0 which satisfies the equation (I — C)x = 
d. 

A consumption matrix C is said to be productive if (1 — C)' 
exists and (1—C)'>0. 


A consumption matrix C is productive if and only if there is 
some production vector x > 0 such that x > Cx. 


A consumption matrix is productive if each of its row sums 
is less than one. A consumption matrix is productive if each of 
its column sums is less than one. 


Now we will formulate the Smarandache analogue for this, 
at the outset we will justify why we need an analogue for those 
two models. 


Clearly, in the Leontief closed Input - Output model, 


pi = price charged by the i” industry for its total output in reality 
need not be always a positive quantity for due to competition to 
capture the market the price may be fixed at a loss or the 
demand for that product might have fallen down so badly so 
that the industry may try to charge very less than its real value 
just to market it. 
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Similarly aj 2 0 may not be always be true. Thus in the 
Smarandache Leontief closed (Input-Output) model (S-Leontief 
closed (Input-Output) model) we do not demand p; = 0, p; can 
be negative; also in the matrix A = (aj), 


aij =F ap; +... taj #1 


so that we permit aj;'s to be both positive and negative, the only 
adjustment will be we may not have (I — A) p = 0, to have only 
one linearly independent solution, we may have more than one 
and we will have to choose only the best solution. 


As in this complicated real world problems we may not 
have in practicality such nice situation. So we work only for the 
best solution. 


Here we introduce a input-output model which has some p 
number of input-output matrix each of same order say n x n 
functioning simultaneously. 


We shall call such models as input — output super row 
matrix models and describe how it functions. Suppose we have 
p number of n x n input output matrix given by the super row 
matrix A = [A,; | ... | A,] where each A; is a 
n X n input output matrix which are distinct. 


A=[A\|...|Agl 
1 1 p 
ay ain an ain 
a a) a a? 
21 2n | | 21 2n 
1 1 p 1 
ant Aon ant aon 
where a,+a,+...¢a,=Lt = 1, 2).20, prand: j= 1, Ze seg: 


Suppose 
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oA fs a) heer 
P, |P, |o | Pa 


be the super column price vector then 


A * P=P, the (product) * is defined as A * P = P that is 


FAP fess AGP a) LPF fees P2] 
A*P=P 
that is 


(I-A) P=(0]... | 0) 
ive., (I— Aj) P; | ... | 1— Ap) Pp) =(0] ... | 0). 


Thus A is an super-row square exchange matrix, then AP = 
P always has a row column vector solution P whose entries are 
non negative. 


Let A = [A, | ... | An] be an exchange super row square 
matrix such that for some positive integer m all the entries of 
A™ i.e. entries of each A? are positive for m; m = 1, 2, ..., p. 
Then there is exactly only one linearly independent solution of 


(I-A) P=(0}...| 0) 
ie., (I— Aj) P; | ... | Z—Ap) Pp) = (0 ... | 0) 


and it may be choosen such that all of its entries are positive in 
Leontief open production super model. 


Note this super model yields easy comparison as well as this 
super model can with different set of price super column vectors 
and exchange super row matrix find the best solution from the p 
solutions got from the relation 


(I-A) P=(0}...| 0) 
Ke =Ad) Pilon USA PO =O ac); 
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Thus this is also an added advantage of the model. It can 
study simultaneously p different exchange matrix with p set of 
price vectors for different industries to study the super 
interrelations between prices, outputs and demands 
simultaneously. 


Suppose one wants to make a study of interrelation between 
prices, outputs and demands in an industry with different types 
of products with different exchange matrix and hence different 
set of price vectors or of many different industries with same 
type of products its interrelation between prices, outputs and 
demands in different locations of the country were the economic 
status and the education status vary in different locations, how 
to make a single model to study the situation. In both the cases 
one can make use of the super input-output model the relation 
matrix which is a input-output super diagonal mixed square 
matrix, which will be described presently. 


The exchange matrix with p distinct economic models is 
used to describe the interrelations between prices, outputs and 
demands. Then the related matrix A will be a super diagonal 
mixed square matrix 


Aj, ..., Ap are the exchange matrices describing the p-economic 
models. Now A acts as integrated models in which all the p 
entities function simultaneously. Now any price vector P will be 
a super mixed column matrix 


P, 
P= 
P, 
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where each 
pe 
Beet h.2.) (eg 
Py, 
fort=1,2,..., p. 


Here each A; is a n; X m exchange matrix; t= 1, 2, ... 


= P is given by 


,p. AP 
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whose entries in each P; are non negative; 1 <t<p. 


Let A be the super exchange diagonal mixed square matrix 
such that for some p-tuple of positive integers m = (mj, ..., m,), 


A‘ is positive; 1 < t < p. Then there is exactly only one 
linearly independent solution; 


and it may be choosen such that all of its entries are positive in 
Leontief open production super model. 


Next we proceed on the describe super closed row model 
(or row closed super model) as the super closed model (or 
closed super model). 


Here we have p sets of K industries which are distributed 
among themselves i.e. the first set with K industries distributed 
among themselves, the second set with some K_ industries 
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distributed among themselves and so on and the p set with some 
K industries distributed among themselves. It may be that some 
industries are found in more than one set and some industries in 
one and only one set and some industries in all the p sets. This 
open super row model which we choose to call as, when p sets 
of K industries get distributed among themselves attempts to 
satisfy an outside demand for outputs. 


Portions of these outputs may still be distributed among the 
industries themselves to keep them operating, but there is to be 
some excess some net production with which they satisfy the 
outside demand. In some super closed row models the outputs 
of the industries in those sets which they belong to were fixed 
and our objective was to determine sets of prices for these 
outputs so that the equilibrium condition that expenditure equal 
income was satisfied for each of the p sets individually. 


Thus we will have 


xX = monetary value of the total output of the i™ 

industry in the t™ set 1 <i<Kand1<t<p. 
gd’ r= monetary value of the output of the i™ 
industry of the t" set needed to satisfy the 


outside demand, | <t<p,I=1,2,...,K. 


Oo, = monetary value of the output of the i” 


industry needed by the j" industry of the t™ 
set to produce one unit of monetary value of 
its own output, 1 <i<K;1<t<p. 


With these qualities we define the production super column 
vector 
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xX, 
X, : 
Th 1 
x= X, = : 
: YP 
ee ay 
a : 
p 
Xk 
The demand column super vector 
1 
d, 
Sy) |s3 
: 1 
eal lp 
P 
= d, 
d, : 
P 
dy 
and the consumption super row matrix C = (C, | ie ie) 
1 1 1 p p 
01, Op OK 91, Op OK 
1 1 1 p p 
a Sar 5 Oox | | Oy, Or OoK 
1 1 1 p p Pp 
Oxi Ox2 Oxk Oxi Ox2 Oxk 
By their nature we have 


‘oS 
i) 


X>/i]sd>] i] andC>(0|...|0). 
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For the t” set from the definition of 0; and x; it can be seen 


that the quantity 
Oj, X, +0), X) +...4+ 0) Xq 


is the value of the i" industry needed by all the K industries (of 
the set t) to produce a total output specified by the production 
vector X;. Since this quantity is simply the i” entry of the 
column vector C,X; we can further say that the ie entry of the 
column vector X, — XC, is the value of the excess output of the 
i” industry (from the t™ set) available to satisfy the outside 
demand. 


The value of the outside demand for the output of the i™ 
industry (from the t™ set) is the i” entry of the demand vector d,; 
consequently we are lead to the following equation for the t™ set 
X, — CX; = d; or (I — C,)X; = d; for the demand to be exactly met 
without any surpluses or shortages. Thus given C, and d; our 
objective is to find a production vector X, = 0 which satisfies the 
equation 


(I = Ci) Xt =d, 


so for the all p sets we have the integrated equation to be 


d-C)X =d 
i.e., [I —C)X: |... | T— Cy) Xp] 
=(dy|! 20 da) 
The consumption super row matrix C = (C; | ... | Cp) is said to 
be super productive if 
(I-C)' =[-C,)" |...) -C,) 


exists and 


(I-C)" =[(I-C,)" | ... | d-C,)“] 2 [0]... 0]. 
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A consumption super row matrix is super productive if and 
only if for some production super vector 


X, 


xX 


n 


Ol: 


such that X > CX ie. [X1 | ... | Xp] > [CiX |... | CoXpl. 


A consumption super row matrix is productive if each of its 
row sums is less than one. A consumption super row matrix is 
productive if each of its column super sums is less than one. 
The main advantage of this super model is that one can work 
with p sets of industries simultaneously provided all the p sets 
have same number of industries (here K). This super row model 
will help one to monitor and study the performance of an 
industry which is present in more than one set and see its 
functioning in each of the sets. Such a thing may not be possible 
simultaneously in any other model. 


Suppose we have p sets of industries and each set has 
different number of industries say in the first set output of K, 
industries are distributed among themselves. In the second set 
output of K, industries are distributed among themselves so on 
in the p" set output of K,-industries are distributed among 
themselves the super open model is constructed to satisfy an 
outside demand for the outputs. Here one industry may find its 
place in one and only one set or group. Some industries may 
find its place in several groups. Some industries may find its 
place in every group. To construct a closed super model to 
analyze the situation. 


Portions of these outputs may still be distributed among the 
industries themselves to keep them operating, but there is to be 
some excess some net production with which to satisfy the 
outside demand. 
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Let 

xX; = monetary value of the total output of the i® 
industry in the t" set (or group). 

d; = monetary value of the output of the i" industry 
of the group t needed to satisfy the outside 
demand. 

om = monetary value of the output of the i" industry 


needed by the j" industry to produce one unit 
monetary value of its own output in the t™ set or 
group, 1 <t<p. 


With these qualities we define the production super mixed 
column vector 


1 

xX, 

X, : 

“eS 1 

Xx, 

xX = X, — : 
. x 

x : 

xP 


dj 
d, di, 
d = = : 
d,j | 4 
dy 
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and the consumption super diagonal mixed square matrix 


C= 
where 
t t t 
01, Op 1K, 
t t t 
Or, On Ork 
C =_ t e 
t * ° ’ 
t 1 
Ox1 Sx,2 Ox k, 


true fort = 1, 2, ..., p. 


By the nature of the closed model we have 


| >< 
a) 
i) 


X=] = [=|7], d= 


~ 

Ol: 
jon 
ol 


and 


From the definition of o;, and x; for every group (set t) it can 


be seen the quantity 6, X,+...0,X, is the value of the 


output of the i” industry needed by all K, industries (in the t™ 
group) to produce a total output specified by the production 
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vector X; (1< t <p). Since this quantity is simply the i" entry of 
the super column vector in 


=[CiX |... | CpXp]' 
we can further say that the i” entry of the super column vector 
X,— CX; in 
Ky OOK, 
X-CX =| : : 
Xx, ~ C,X, 


is the value of the excess output of the i industry available to 
satisfy the output demand. 


The value of the outside demand for the output of the i® 


industry (in t™ set / group) is the i" entry of the demand vector 
d,; consequently we are led to the following equation 


X,— CX, = d, or (I; — C,) Xi = d, (1 <t< P), 


for the demand to be exactly met without any surpluses or 
shortages. 


Thus given C, and d, our objective is to find a production 
vector X; = 0 which satisfy the equation (I, — C,)X;= d. 


The integrated super model for all the p-sets (or groups) is 
given by X -CX =d_Le., 
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X, ~ CX, d, 
X, C,X, d, 
X,-C,X,] |, 
or 
X, d, 
x1. (ed, 
1.€., 
(I, z C,)X, dq, 
(I, - C,)X, d, 


where I is a K, x K; square identity matrix t = 1, 2, ..., p. 


Thus given C and d our objective is to find a production 
super column mixed vector 
Ke 0 


ye eae 


x 


Pp 


Ole 


which satisfies equation (I— C) X =d 


(I, —C,)X, dq, 
ie. : = ee 


(I,-C,)X, } {d 


P 


A consumption super diagonal matrix C is productive if 
(I-C)' exists and ie., 
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exists and 


(I,-C,)" 


IV 


A consumption super diagonal matrix C is super productive if 
and only if there is some production super vector 


X, 0 
Ke >] 
X, 0 
such that 
X C,X, 
X>CXie.| : | > ; . 
X, CX, 


A consumption super diagonal mixed square matrix is 
productive if each row sum in each of the component matrices is 
less than one. A consumption super diagonal mixed square 
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matrix is productive if each of its component matrices column 
sums is less than one. 


The main advantage of this system is this model can study 
different sets of industries with varying strength simultaneously. 
Further the performance of any industry which is present in one 
or more group can be studied and also analysed. Such 
comprehensive and comparative study can be made using these 
super models. 


Finally we develop the notion of super special codes. 


Here we define two new classes of super special row codes 
using super row matrix and super mixed row matrix. 


DEFINITION 4.18: Suppose we have to transform some n set of 


non 


| ome | 1 222: 2 n 
ky, «1, ky message symbols a, d}...d, 5 Aj Qs ..Aj, 5.45 A; A5 Ay , 
t : . . 
a, €F,; 1 St Snand 1 Si Sk; (q a power of a prime) as a set 


of code words simultaneously into n-code words such that each 
code word is of length n;, i = I, 2, .... n and n;— k; = n2— kp = 

. = m,—k, =m say i.e., the number check symbols of every 
code word is the same i.e., the number of message symbols and 
the length of the code word may not be the same. That is the 
code word consisting of n code words can be represented as a 
super row vector; 


-_ y iaey 1 2). 2. yi non n 
Ky =] pH yee Xq, | Xp Xy vey, | vee | Hp Ay oo Hy 


1 


ni> k;, 1 Si <n. In this super row vector ie - a, Ilsjsk; i= 


se i i i 
I, 2, ..., n and the remaining n; — k; elements x, ., X49 ++ Xy, 


are check symbols or control symbols; i = 1, 2, ..., n. 


These n code words denoted collectively by x, will be 
known as the super special row code word. 
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As in case of usual code, the check symbols can be obtained 
in such a way that the super special code words x, satisfy a 


super system of linear equations; H,x; =(0) where H, is a 


super mixed row matrix given by H, = [H; | H2| ... | Ha] where 
each Hj is a m x nj matrix with elements from F,, i= 1, 2, ..., n, 
1.€., 


T 
H,x? =[Hp|Ho| ss: | He] |x, x2 x? | 


H,(x2)' |... | (x2)" | 


= [I(0) | (0) |... |] 


i.e., each H; is the partity check matrix of the code words a a 


= 1,2, ...,n. Hs = [H;| H2| ... | Hs] will be known as the super 
special parity check super special matrix of the super special 
row code C,. C, will also be known as the linear [(n; nz ... ny), 
(k; ky ... k,)] or [(m1, ky), (mo, kp), ..., (On, Kn)] super special row 
code. 


If each of the parity check matrix H; is of the form 
(A,,1 ies pee 


i? ~n;—k; 


H, = [H;| H2|... | Aa] 
= (Ait. ) | (Aula, ) | Ge | (Ady ce )] ral (1) 


C, is then also called a systematic linear ((n; m2 ... Ny), (ki ko ... 
k,)) super special code. 


If q = 2 then C, is a super special binary row code; (k; + ... 
+ k,) | (n; + nm) + ... + ny) is called the super transmission(or 
super information) rate. 


It is important and interesting to note the set C, of solutions 
x, of H,x! = (0) i.e., known as the super solution space of the 
super system of equations. Clearly this will form the super 
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special subspace of the super special vector space over F, of 
super special dimension (k, ky ... k,). 


C, being a super special subspace can be realized to be a 
group under addition known as the super special group code, 
where H, is represented in the form given in equation I will be 
known as the standard form. 


Now we will illustrate this super special row codes by some 
examples. 


Example 4.1: Suppose we have a super special binary row code 


given by the super special parity check matrix Hs = [H; | H2| Hs] 
where 


and 


1.e., the super row matrix associated with the super special code 
is given by 
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= [(Ai, Is) | (Ao, Js) | (As, I; 


i.e., we have given the super special code which is a binary row 
code. The super special code words are given by 


Pb ot oh ab Sah po Fe 2 oP 
x, =[a) A, a; X, X; X_| a; A, a; ay Xs X_ x; | 


x, Ee 


8: 28 pO oli S ipo! eee gest fe 1 
Ay G5 Ae A de x} |= [ x! 


H,x! = (0) gives 3 sets of super linear equations ie., H,x! = 
(0) is read as 


[H: | H2| Hs] BS Bowles | -|H, (x!)’ | H, (x?)" 
= [(0) | (0) | (OI: 


H, (x!)' =@) 


results in linear equations solving which we can get the code 
words. 


Thus the number of super special row code words in this 
example of the super special code C, is 8 x 16 x 32. 


The super transmission rate is 12/21. Thus this code has 
several advantages which will be enumerated in the last chapter 
of this book. We give yet another example of super special code 
in which every super special code word is a super row vector 
and not a super mixed row vector. 


Now we proceed on to define the notion of super special 
row repetition code. 


DEFINITION 4.19: Let C,=|C/ | C? 


| Cr | be a super 


special row code in which each of the C is a repetition code, i 
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= 1, 2, ..., n, then we define C, to be a super special repetition 
row code. Here if H, = [H;|H>| ...|H,] is the super special parity 
check matrix of C,, then each H; is a t— 1 xt matrix that is we 
have 


A, =A)=... =H, =|. 
ESO: Ease. uh 
t-I1 xt 
is the parity check matrix. The super special code words 
associated with C, are just super row vectors only and not super 


mixed row vectors. The number of super special code words in 
Geis 2 


We illustrate a super special row repetition code by the 
following example. 


Example 4.2: Let C, tC, ee? Cc 


repetition code with associated super special row matrix H, = 
[Hi | Hp | H3| Ha] 


ont be a super row 


ot: OO Ot Oe 
es Re 8 Pe 
Tr Os Oath. OTL Ord: 70 
10001/10001 
1 O26: [LR 2 O00 
1010 0/1010 0 
1001 0/1001 0 
iO OO 1/120) O <0 4 


DEFINITION 4.20: Let C, be a super special parity check mixed 


row code i.e., C. =[C! ecg a) where C, is obtained 


Ss 
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using the super special mixed row matrix H, = [H,| H2| ... | Hn] 
where each H; is a unit row vector having t; number of elements 
1.e., 


E Tos ak los ele i J 
HS) nd ee 


t, times t, times t, times 


where at least one t; #t; for i #j. Any super special code word in 
C, would be of the form 


nn ni} 1 
Xx, X>5 Xy | = E 


22: 2 
Xj X71 Xi, E 


oa [oe I 
a =| Ky wos Xy, 
‘i T 5 i ‘ 
with H,x, = (0); i.e., each x. would contain only even number 
of ones and the rest are zeros. 


C, = [C, | Co | ... | Cy] is defined to be super special parity 
check row code. C, is obtained from the parity check row matrix 
/ vector H, = [H; | H. | ... | H» ] where H; = H,=... = H, 


= E 1 --- 1]. Here a super special codeword in C, would be 


m times 


a super row vector of the form [ x! lesa ee xt | with each 
<= Es re a where only even number of x; are ones and 
the rest zero, 1 <j <mandi=1, 2,...,n. 


Now we will illustrate the two types of super special parity 
check (mixed) row codes. 


Example 4.3: Let C, =|c c 


check code having H, = [H;|H2|H3J=[1 111 1/111]/1111 
1] to be the super special parity check matrix associated with it. 


c}| be a super special parity 


DEFINITION 4.21: Let C,=[C! | C? 


ia) oa be a super 


special row code. Suppose x, =(x! i ee | x" jis a 
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transmitted super code word and y, = [ y, 


Pelee get [as 
the received supercode word then eg = yy, — XxX, = 
[y’ — x! sella -x"|=[e! e" | is called 
the super error word or the super error vector. 


De SB 2 
Vs —% @; 


RY 


We first illustrate how the super error is determined. 


Example 4.4: Let C, eke: | C2 |e? | (en be a super special 


code with associated super parity check row matrix H, = [H; | H 
| H3| Hy) 


Clearly y, + e, = Xs. 


DEFINITION 4.22: The super Hamming distance d,(x,, y,) 
between two super row vectors of the super special vector space 
V,, where x,=[x' |x? |... | x" Jandy, =[y' |? |... | x7] 


is the number of coordinates in which x and y'. differ for i = 


1, 2, ..., n. The super Hamming weight w,(x;) of the super vector 


Xs = [ x! pal er | x] in V, is the number of non zero 
coordinates in each ba ,i=1, 2, .., 0. In short w(x) = d(xs, 
(0). 


As in case of usual linear codes we define super minimum 
distance d*,, of a super special linear row code C, as 


min 
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so : 
Gam — oe d, (u,v, ) ? 
“Uy, s 


d, (Us, Vs) = ds (Us — Vs 5 (0)) = Ws (Us — Vs). 


Thus the super minimum distance of C, is equal to the least 
super weights of all non zero super special code words. 


Now the value of 


d\n = min d, (Us, Vs) 
u, ,V, € C, 
U, #V, 
= ’ 1 2 n 1 2 n 
= min. d,([u! ts u® |,[ v! Me lo vj) 
eee 


= min [d(ul, vi )+d(uz, vz)+...+d(u?, vt) ]. 


Now di‘, of the super special row code given in example 
3.1.2 is 7 verified using the fact in C,=[C! | C? | C2]; diinC: 


=3, di,,C; =2 and d},,C; =2. Hence d',,(C,) =3+2+2= 


7. So we will denote d*,, =min d,(u,,v,) by d;,(C,), Us. Vs 


min 


€ Cy, Us # Vs. 


Fl CE | sal C2 ah, | (CHC et (CP) 
= min, d(x, yi)+ ymin d(x:, yz) +...+ min d(x?, ys). 
Xs.¥s EC. X50Ys EC. xP oYs EC! 


| er a noyyn 
Xs¥#Ys Xs #Ys Xe #Ye 


Now we proceed on to define the dual of a super special 


row code. 
DEFINITION 4.23: Let C,=[C!|C?|...| C"] be a super 
special row [(n}, ..., My), (ki, ..., Ky)] binary code. The super 


special dual row code of C, denoted by 


eal(ey (ey ||) 
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where (cy = {ui |ui-vi=0 forall vi eC }i=1, 2, .., 0. 


Since in C, we have n; — kj = no — kp = ...= ny, — ky, i.e, the 


number of check symbols of each and every code in C' is the 
same for i = 1, 2, ..., n. Thus we see n = 2k; alone can give us a 
dual, in all other cases we will have problem with the 
compatibility for the simple reason the dual code of C' being 
the orthogonal complement will have n;— k; to be the dimension, 
where as C' will be of dimension k;, i = 1, 2, ..., n. Hence we 


can say the super special dual code would be defined if and only 
ifn; = 2k; and such that n; = nz = ... = Mp. 


We can define the new notion of super special syndrome to 
super code words of a super special row code which is 
analogous to syndrome of the usual codes. 


DEFINITION 4.24: Let C, be a super special row code. Let H, be 
the associated super special parity check matrix of C, the super 
special syndrome of any element y, € V; where C, is a super 
special subspace of the super special vector space V, is given by 


Sv) = H,y>. S(v) = (0) if and only if ys € C;. 


Thus this gives us a condition to find out whether the 
received super code word is a right message or not. Suppose ys 


is the received super special code word, we find S(y,) = H,y’ ; 
if S(v,) = (0) then we accept y, as the correct message if S(ys) = 
a y #(0) then we can declare the received word has error. 


We can find the correct word by the following method. 
Before we give this method we illustrate how the super special 
syndrome is calculated. 


Example 4.5: Let C,=[C! | C? | C 


row code. Let H, = [Hj | H2 | H3| H4] be the super special parity 
check matrix of C,. 


en be a super special 
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Let 
1 0 0 010100 
H, =| 1 1 0 10101 0 
0 000 1j;0 121 00 ~=1 
0 1 00;0 110100 
1 0 0/1 10410 1 0 
1 1 01;0 0000 0 1 


Now we have to shown how to find whether the received 
super special code word is correct or otherwise. It is important 
to note that what ever be the super special code word x, € C, 
(i.e., it may be a super special mixed row vector or not) but the 
syndrome S(x,) = H,x! is always a super special row vector 


which is not mixed and each row vector is of length equal to the 
number of rows of H,. 


Now we know that every super special row code C, is a 
subgroup of the super special vector space V, over Z, = {0, 1}. 
Now we can for any x, € V, define super special cosets as 

x, + C,= {x, +c, | c, € C,}. 


Thus 


V5 = {Zp x Zy x...x Zy| Zo x... X Zy| ... |Z x Zp x... X Zp} 
=C,U [x, $C] Ui. U0 [x4 CG] 


where G=(C |e 


and 


II 
7 
can 


and 
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metC,= [xl +C! |x? +C? |... 


n n 
x" +2]. 


Now we can find the coset leader of every x! +C! as in 


case of usual codes described in chapter one of this book. Now 
if 


aly] 


ys=[ys | yz 


e? +(0)|...[e" +(0) | isa 
special super coset leaders then using the relation y, — e, we get 
ys — €, to be super special corrected code word. It is interesting 
to note that each e! + (0) has a stipulated number of coset 


is the received message and [e! +(0) 


leaders depending on n;, i= 1, 2, ..., n. 


We will illustrate this by the following example. 


Example 4.6: Let C, =¢ e | be a super special row code. 


Suppose H, = [H; | H»] be the super special row matrix 
associated with C,. Let 


1041 0/1 041 41 +0 
H,= [H;| Ha] = 


L 2° O° 1/0 1-4) @ 4 
Now C=| Cle | with C! = {0 000), (1011),(0101), 
(111 0)}and C? = {(00000),(10010),(01001),(011 
1),(11011),(1110),(10101),(11100)}. 


C, = {[0000|00000],[1011|00000}, 
0101/00000},[1110/00000}, [0000|/10010 


[ [ ] 
(0000|01001],[1110|01001],[0101|0100 1], 
[1011|10001],[0000|00111),,f1110|0011 1], 
(0101/0011 1],[1011|00111],[0000|1101 1, 
[1110|1101 1], [1011/1101 1],f0101|1101 1], 
(0000|01110],[1011|01110],[0101|0111 0], 
(0000|10101],[1011|10101],[0101|1010 1], 
(0000|11100],[1110|11100],f1011|11100], 
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[0101 |11100] and so on}. 


Clearly |C,| = 32 . Now the coset table of C! is given by 


Message code words 
0 0 1 0 0 1 ie i 
00 0 0 1011 010 1 111 


Other cosets 


10 0 0 0011 1101 01410 
010 0 1111 000 1 101 0. 
0010 100 1 01411 110 0 
fe iw 8g 
coset leaders 
Now the coset table of C? is given by 
message 0 0 0 100 
codewords 00000 41001 9 
1? 0 OO 2 DO TO 
othercosets <0 1 00 0 171 101 0 
Oo) OE ie as Ge ah aly 
message 01 0 00 1 
codewords 0 1001 9011 41 
1 10041 10o041 1 1 
other cosets +1 0001 04111 41 
01101 000141 
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— 


message 1 0 01 1 
codewords 1 1011 0141 4110 
01011 411110 
othercosets 41 0011 0041410 
111i1t1 04310140 
message 1 0 1 ly I 1 
codewords 1 0101 11100 
O20: We Oe o> par ST OO 
othercosets 41 1 101 7310100 
10001 11000 


Suppose y, = [1 1 1 1|1 111 1] is the received word then 
S(ys) = H,y! ¥ [(0) | (0)]. e, = [0 1 00|00 1 0 0] is the super 
set coset leader. Thus x,=y,t+e,=[1011]/11011] eC, 


With the advent of computers calculating the super special 
coset leaders is not a very tedious job. Appropriate programs 
will yield the result in no time. 


Now we proceed on to describe/define the super special row 
cyclic code. 


DEFINITION 4.25: Let C, = [C; | C2 | ... | C,/ be a super special 
row code. If every C! is a cyclic code in C, we call C, to be a 


super special cyclic row code. H, = [H,| H>| ...|H,] denotes the 


super special parity check row matrix of the super special cyclic 
code. 


We illustrate this by the following example. 


Example 4.7: Let C, = Ker | C2 


Cc | be a super special cyclic 
code with an associated super special parity check matrix 


H, = [H, | H2| Hs] 
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1 101 0 0 

=|0 1 1 01 0 

001101 
00111 i1;0 01001 
0101 1 0/0 100 1 0 
101 10 0;1 00 1 0 0 


We see each of C!, C2 and C. are cyclic codes. 


Now we see in general for any super special mixed row 
code with an associated super special parity check matrix H, 
which happens to be a super mixed row matrix we cannot define 
the super special generator row matrix G,. The simple reason 


being if the code words in each of the C! in C, where 
C=C Colas Cis i= 1, 2, ..., n happens to be of 
different length then it would be impossible to define a super 
generator matrix. So we shall first define the notion of super 


special generator row matrix of a super special row code(mixed 
row code). 


DEFINITION 4.26: Let C,=[C! | C? 


én | ce be a super 
special row code. A super special row matrix which generates 
C, exists if and only if in each C' the codes in C, have the same 
number of message symbols, that is if C, has a super special 
parity check row matrix H, = [H; | H) | ... | H,] then we 
demanded each C: must have the same number of check 
symbols. Likewise for the super special generator row matrix to 
exist we must have G, = [G; | G2 | ... | G,/ where Cc have the 


same number of message symbols which forms the number of 
rows of the super row generator matrix Gs. 


We shall first illustrate this by an example. 
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Example 4.8: C, =[C | oF | be a super special row code. 


Let G, =[G! 


G?] 


One can easily find the super code words using Gs [50]. 


DEFINITION 4.27: Let C,=[C! | C? 


ae | (oa be a super 
special mixed row code. If each of the codes C' have the same 
number of message symbols then we have the super special 
generator mixed row matrix G, =(G! | G | aed | G" | 


associated with C,. Number of message symbols in each of the 
G! are equal and is the super special mixed row matrix G,; 1 < 


in. 
We illustrate this by the following example. 


Example 4.9: Let C,=[C! | C?.| Cc 
mixed row code. Let G, =[G! | G? 


super special mixed row generator matrix given by 


ci] be a super special 


G?] be the associated 


G,=[G' |G? |G] 
10041 1 £0 
=/0 1001 1 
0011 0 41 
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oor 
oro 
a) 
oro 
a) 
oor 
a) 
oor 
oro 
a) 
— 
oro 
a) 
oor 


THEOREM 4.4: Let C, =e. |C? ee | C" |be a super special 


row code with H, = [H,| H2| ...| H,] , the super special parity 
check matrix. If each H; = (Aj, I,4), i = 1, 2, ..., n then G, = [G; 
| G) | ...| G,J with G; = (, — A"); 1 Si <n if and only if the 


length of each code word in C is the same fori = I, 2, ..., n. 


For proof refer [50]. 


For any super special row code we can Hg the parity check 
matrix and the super special generator matrix Gs such that 


Gs H§ = (zero super vector). 


Now having defined the new class of super special row 
(mixed row) codes we will now define new classes of mixed 
super classes of mixed super special row codes C, i.e., we may 
have the super special row code to contain classical subcodes as 
Hamming code or cyclic code or code and its orthogonal 
complement and so on [50]. 


DEFINITION 4.28: Let C, = [C; | C2 | ...| C,/ be a super special 
row code. If some of the C;’s are Hamming codes, some C;’s are 
cyclic codes i #j, some C's are repetition codes and some C;,’s 
are codes and C,’s are dual codes of C,’s ; 1 Sj,k, tii p<n 
then we call C, to be a mixed super special row code. 


It is important to mention here that even if two types of 
classical codes are present still we call C, as a mixed super 


special row code. 


We will illustrate them by the following examples. 
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Example 4.10: Let C, = [C, | C2 | C3 | C4] be a mixed super 
special row code. Here C; is a Hamming code, C; the repetition 
code, C3 a code of no specific type and Cy, a cyclic code. 


Let the mixed super special parity check matrix H, 
associated with C, be given by 


H, = [H;| H2 | H3| Ha] 


1 1} 1 
1 1}1 
0 1 


1 
0 
1 1 


Example 4.11: Let C, = [C; | C2 | C3] be a mixed super special 
row code. Let H, = [H; | H2 | H3] be the associated super special 
parity check mixed row matrix. C; is the Hamming code, C, any 
code and C; a repetition code. 


1 11112 21éi%idy1 01 21 0 0 0 
00011 21 21 0;/1 +1210 1 0 0 
lo 110011 0/1 00001 0 
1010101 0;0 12100 0 ~41 
110 0 0 
101 0 0 
100 1 0 
100 0 1 


is the mixed super special parity check mixed row matrix for 
which G, does not exist. 
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We define the new notion of super special Hamming row 
code. 


DEFINITION 4.29: Let Cs = [C; | C; |... | Cy] where each 


C! is a (2"— 1, 2"—1-—m) Hamming code for i = 1, 2, ..., n. 
Then we call C, to be a super special Hamming row code. If 
A= [H, | H2 | ... | Hy] be the super special parity check matrix 
associated with C, we see H, is a super special row matrix 
having m rows and each parity check matrix H; has m rows and 
2” — 1 columns, i =1, 2, ..., n. 


Further the transmission rate can never be equal to 4. If 
m > 2 then will the transmission rate be always greater than 2 ? 


For more examples refer [50]. 


We define using super column matrices the notion of super 
column codes which is an innovative means of using super 
column matrices. 


DEFINITION 4.30: Suppose we have to describe n codes each of 
same length say m but with varying sets of check symbols by a 
single matrix. Then we define it using super column matrix as 
the super code parity check matrix. Let 


fALO 


Pore 


be a set of m codes, Cj, C2, ..., Cn where all of them have the 
same length n but have n—k), n— kp, ..., n— Ky to be the number 


of check symbols and kj, kz, ..., km are the number of message 
symbols associated with each of the codes C), Cy .... Cn 
respectively. 


Let us consider 
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H! 
le 
H™ 


where each Hi is the n — k, x n parity check matrix of the code 
Ci; i= 1, 2, ..., m. We call H* to be the super special parity check 
mixed column matrix of C, and C, is defined as the super 
special mixed column code. 


The main difference between the super special row code 
and the super special column code is that in super special row 
codes always the number of check symbols in every code in C, 
is the same as the number of message symbols in C; and the 
length of the code C; can vary where as in the super special 
column code, we will always have the same length for every 
code C; in C, but the number of message symbols and the 
number check symbols for each and every code C; in C, need 
not be the same. In case if the number of check symbols in each 
and every code C; is the same. Then we call C, to be a super 
special column code. 


In case when we have varying number of check symbols 
then we call the code C, to be a super special mixed column 
code. 


In the case of super special column code C, = [C; | C2 | ... | 
Cm]' we see every code C; in C, have the same number of 
message symbols. Thus every code is a (n, k) code. It may so 
happen that some of the C; and C; are identical codes. 


For examples refer [50]. 
DEFINITION 4.31: Let C, = [C;| C2 | ... |C,]' where C,’s are 


codes of same length m. Suppose each C; is generated by a 
matrix Gj, i = 1, 2, ..., n, then 
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generates the super special column code C,. We call G, the 
super special generator column matrix which generates C,. If in 
each of the codes C; in C,, we have same number of message 
symbols then we call G, to be a super special generator column 
matrix; i = 1, 2, ..., n. If each of the codes C;’s in C, have 
different number of message symbols then we call G, to be a 
super special generator mixed column matrix. 


We say G, is in the standard form only if each G; is in the 
standard form. Further only when G, is in the standard form and 
G, is a super special column matrix which is not a mixed matrix 
we have H, the super special parity check column matrix of the 
same C, with 


Joo 


G,H! = 


Ol-:: 


For examples refer [50]. 


We proceed onto define classical super special column 
codes. 


DEFINITION 4.32: Let C, = [C; | C2 | ...| C,/‘ be a super special 
column code if each of the code C; is a repetition code of length 
n then C, is a super special repetition column code with C; = 
Co peared Cy. 
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The super special column parity check matrix 


where 


S 


an) 
ll 
a) 


LP QoG Oy 4 “Ff 
It is important to note that unlike the super special 
repetition row code which can have different lengths the super 


special repetition column code can have only a fixed length and 
any super special code word 


2 
x; 


n 
[xt] 


where x! = (1 1 ... 1), n-times or (0 0 ... 0) n-times only; 
1 Sj Sn. 


_— v4 
Xs = [ x! 
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DEFINITION 4.33: Let C, = [C;| C> | ...| C,]' be a super special 
parity check column code . Let the super special parity check 
column matrix associated with C, be given by 


Hi, 
eee 
H, 
Where 
H,=H,=... =H, Chr: 1): 
uU“«.“—“ 


m-—times 


Thus we see we cannot get different lengths of parity check 
codes using the super special column code. However using 
super special row code we can get super special parity check 
codes of different lengths. 


DEFINITION 4.34: Let C, = [C; | C2 | ... | C,]' be a super 
special column code if each of the codes C; in C, is a (2" — 1, 2” 
— 1—m) Hamming code for i = 1, 2, ... , n then we call C, to be 
a super special column Hamming code . It is pertinent to 
mention that each code C; is a Hamming code of same length; i 
=1,2,...,N. 


DEFINITION 4.35: C, = [C; | C2 | ...| C,/' is a mixed super 
special column code if some C;’s are repetition codes of length 
n some C;’s are Hamming codes of length n, some C,’s parity 
check codes of length n and others are arbitrary codes, 1S i, j, 
k Sn. 


We illustrate this by the following example. 
DEFINITION 4.36: Let C, = [C; | C> | ... | Cn]' be a super 


special column code if each of the codes C; is a cyclic code then 
we call C, to be a super special cyclic column code. However 
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the length of each code C;; i = I, 2, ..., n will only be a cyclic 
code of length n, but the number of message symbols and check 
symbols can be anything. 


Now we illustrate this by an example. 


DEFINITION 4.37: Let 


m 
8 |e 


m m 
C; C; 


where Cc are codes I $i Smand 1 <j Sn. Further all codes 
(Gh thd ..C)" are of same length Cc. ...,C;' are of same 
length and C',C?,...,C% are of same length. C/,C},...,C/ 
have same number of check symbols, C;,C3,...,C? have same 


number of check symbols and C;",C3',...,C” have same number 
of check symbols. 


We call C(S) to be a super special code. We can have the 
super parity check matrix 


H; 
Hy 


H(S) = 


Ay 
where H; ’s are parity check matrices 1 $i <mand 1 <j Sn. 
Further Hj,H},...,H; have the same number of rows, 
H;,H3,...,H; have same number of rows and so on. 
Hy ,H;',...,H)” have the same number of rows. Likewise 


1 2 
H,,H,,...,H;" have the same number of columns, 


148 | Innovative Uses of Matrices 


1 2 
H,,H;,...,.H;' have the same number of columns and so on. 


H!,H?,...,H” have same number of columns. 

The notion of super matrices are ingeneously applied to 
super fuzzy models, super linear algebra and super special 
codes. Such type of applications is not only innovative but can 
lead to easy comparison and save the working time. 


Chapter Five 


INTERVAL MATRICES AND NATURAL 
CLASS OF INTERVALS 


In this chapter first we recall properties of intervals and then 
define the notion of matrix interval using matrices with entries 
from R. Also we define the matrix of natural class of intervals 
and illustrate them with examples we however get the relation 
between the interval matrix and matrix interval and the relation 
between the natural class of matrix intervals and the matrix of 
natural class of intervals. We will call [A, B] an interval where 
A>BorA<BorA =B or A and B are not even comparable 
but are bound by some common features. When A and B are 
integers or modulo integers we get intervals which we call them 
as natural class of intervals. For if [15, 20] where 15 and 20 € 
Zy5 is the situation when we cannot in any way compare them. 
So if we replace A and B by matrices of same order we call that 
interval as matrix interval. 


With this rough idea now we proceed onto define the row 
matrix interval. 
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By an interval [a, b] we mean all elements x such that a < x 
<b are in [a b]. [a a] is nothing but a. So when in an interval 
a= b it collapses to a single point. 


Let Q;= {[a, b] | a,b € Q}. 


Let Q/= {[a, b] | a, b € Q” U {0}, Q” denotes the set of 
positive rationals}. 


Z; = {[a, b] | a, b € Z}; denotes the set of all intervals on Z, 
i.e., the set of positive and negative integers. 


Zi = {[a, b] | a, b € Z* U {0}; Z* the set of positive 
integers}. 


R; = {[a, b] | a, b € R} ie., the collection of all intervals in 
the set of reals. 


R{= {[a, b] | a, b € R® U {03, R’ the set of positive reals}. 
Z., = {[a, b] | a, b € Zp}. 


Thus we have defined the collection of all intervals of 
different types. We first make the following observation which 
is very essential. 


We see Qc Qi (a= [a, b] € Q). 


Likewise Z c Z;,Q’ ¢ Q*, Z, ¢Z, Rr DR,R’c RY. 


Thus these which will be known as interval sets are 
generalized or contain the related sets as subsets. We shall call 
Q; the rational intervals, Z; the integer intervals, Rj the real 


intervals and Z\ the modulo integer intervals. 


Further Z c Z; CQ c Ry 
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ZECQCQ ER, ZCQOCRER,. 


[—7, 3] is an integer interval in Z). 

[3, 5] is a positive integer interval in Z; . 
We see Z’ CZ} CZ. 

[—5/2, 3/7] is a rational interval in Q). 


[9/7,5/2] is a positive rational interval in Q* ; Q°c QF <Q, 


[2 a5 | is a real interval in R; but ree V7] is a 
positive real interval in R}. ThusR° c R} CR). We see [0, 1] 


[0, 4] and modulo integer intervals of Z‘. 


It is important and interesting to note Z! has only finite 
number elements in it. 


Infact Z has only four elements viz. {[0, 0], [1, 1] [0, 1]}. 


Z;, has only 6 elements viz {[0, 0] [1, 1], [2, 2], [0, 1], [0, 
2], [1, 2]}. 


Z,, has only 10 intervals viz. {[0, 0], [1, 1], [2, 2], [3,3]. [0, 
1); (8,21. (0,3), (1; 21. (1, 31, (2.37: 

Z; has only 15 modulo intervals. {[0, 0], [1, 1], [2, 2], B. 
3], [4, 4], [0, 1], [0, 2], [0, 3], [0, 4], [1, 2], (1, 3], (1, 41, [2, 3], 
[2, 4], [3, 4]}. 


Z;, has only 28 modulo intervals. 
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Thus Z' has only n+ (n-l) +... + 2+1 = modulo 


n(n +1) 

2 
intervals. 

Zis has only {[0, 0], [1, 1], [2, 2], ..., [14, 14], [0, 1], [0, 2], 
..-, [0, 14], [1, 2], ..., [1, 14], [2, 3], ..., [2, 14], ..., (13, 14]} = 
Loe 

Number of intervals in Z}, is 120. 

Now we recall the operations on intervals given by |26] 


DEFINITIONS [26] 


=[c, d] if and only ifa= candb=d 
Cc [e, dl ifonly ife<a<b<d. 

[a, b] < [c, d] if and only if b<c. The width of an interval 
[a, b] is denoted by w([a, b]) = b — a and the magnitute of the 
interval |[a, b]] = max ([al, |b]). 


Degenerate intervals are those intervals which has zero 
width, can be identified with real numbers. Thus a = [a, a]. 


We now proceed onto give the arithmetic operations as 
given by Kuperman. 


If * denotes any one of the symbols +, -, ., | (. Product, / 
division). 


Then the interval [a, b] * [c, d] = {sey 


axx<b 
except 
c<y<d 


we do not define [a,b] | [c, d] if o € [c, d]. (This excludes the 
definition by zero possibility). Alternatively we can define the 
arithmetic operations on intervals by giving the endpoints of the 
intervals resulting from the sum, difference, product or quotient 
of two intervals. 

Equivalent to the definition in I we thus we 
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[a, b] + [c, d] = [atc, b+d] 

[a, b] — [c, d] = [a—d, b—c] 

[a, b] . [c, d] = [min (ac, ad, be, bd), max (ac, ad, bc, bd)] 
[a, b] / [c, d] = [a, b] / [1/d, 1/c], provided 0 ¢ [c, d]. 


For the denominator cannot be zero. Thus in the definition 


of division c and d are either both positive or both negative. By 
way of numerical example we thus have 


(i) [-2, 3] +[-9, -3] =[-11, 0] 
[8, 11] + [5, 7] =[13, 18]. 
[-3, 8] + [3, 4] =[0, 12] 


(ii) [-3, 8] -[-7, 2] =[-3 — 2), 8+2] 
=[-1, 10] 


(8, 17/2; 6] = (8-6, 17=©)] 
= [2, 19] 


[2, 7] - [-1, 3] = [2-3, 74-1)] 
=[-l, 8] 


(111) [3, 5]. [7,9] =[min {21, 27, 35, 45}, 
max {21, 27, 35, 45}] 
= [21, 45] 


[-2, 1]. [2,5] =[min {-4, 2, 5, -10}, 
max {—4, 2, 5, -10}] 
= [-10, 5] 


[-5, 2]. [3,7] =[min {-15, 6, 14, -35}, 
max {—15, 6, 14, -35}] 
= [-35, 14] 


[-2, 7] /[3, 8] = [min {-2/3, -2/8, 7/3, 7/8}, 
max {—2/3, -2/8, 7/8, 7/3}] 
= [-2/3, 7/3] 
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[-5, 4] / [-2,6] =[min {5/2, -5/6, 4/2, 4/6}, 
max {5/2, -2, 4/6, -5/6}] 
=([-2,5/2]} 


[-8, 9]/[3, 7] = [min {-8/3, 8/7, 3, 9/7}, 
max {—8/3, 8/7, 3, 9/7} 
= [-8/3, 3] 


Now we recall the properties of interval arithmetic given by 
Kooperman. 


For more about these notions please refer [26] 


We may note that interval arithmetic is a generalization or 
extension of real arithmetic since [a, a] is a real number and the 
definitions hold for intervals of this form as well. 


Thus laws that hold for interval arithmetic must 
automatically hold for real arithmetic. But laws that hold for 
real arithmetic may not always hold for interval arithmetic. 


For example, interval arithmetic is associative and 
commutative with respect to addition and multiplication. Thus, 
given interval numbers I, J, K we have 


(I, J)+K=1+(J+K) 
I+J=J+I 


But distributive laws does not always hold for interval 
arithmetic. 


3, 5] ({1, 2] + [-7, 2]) 
= [3, 5] [1, 2] + [3, 5] [-7, 2] 
=[min {3, 5, 6, 10}, max {3, 5, 6, 10}] + 
[min {—21, 6, —35, 10}, max {—21, 6, —35, 10}] 
= [3, 10] + [-35, 10] 
= [-32, 20] = [3, 5] [1, 2] + [3, 5] [-7, 2] 
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Consider [3, 5] ({1, 2] + [-7, 2]) 
= [3, 5] [-6, 4] = [min {-18, —30, 12, 20}, 
max {-18, —30, 12, 20}] 


= [-30, 20] 
Thus [3, 5] ([1, 2] + [-7, 2] # [3, 5] (1, 2] + 3, 5] 7, 2]. 


It is pertinent to mention here that the distributive laws hold 


good when the interval is a real number say k; i1.e., the 
degenerate interval [k, k] of zero with, then it is easy to verify 


[k, k] ([a, b] + [c, d]) = [k, k] [a, b] + [k, k] [c, d] 
= k([a, b] + [c, d]) =k [a, b] + k [c, d] 

=k [atc, b+d] = [ka, kb] + [ke, kd] 

= [k(atc), k(b+c)] = [ka + kc, kb + kd] 

= [ka + kc, kb + kd] 


Thus we have 
b [a — 6a, a + da] = b ([a, a] + [-a, da]) 

= b[a, a] + b[—da, da] 
[ba, ba] + [-|b] 6 a, |b| da] 

= [ba — |b|da, ba + |b| da] 
b[a — 6a, a+ da] = [ba — |b] da, ba + |b|da] 


(We may note that da is non negative by implication of its 
appearance in [a — 6a, a+ da] for in any interval [c, d], c < d). 


We now proceed onto recall the concept of rounded interval 
arithmetic. An important use of interval arithmetic is to 
determine upper bounds for the errors due to rounding. 


In rounded interval arithmetic working in either single or 
double precision, we regard each real number as an interval. An 
in the first place we round the lower endpoint down and the 
upper end point up. 
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In this way, the maximum possible effect of rounding is 
taken into account, so that the interval for the final result 
contains the true value. 


In using rounded interval arithmetic it should be clear that 
the final result can be obtained as sharply as desired, i.e., the 
width of the final interval can be as small as desired, by using a 
sufficiently long wordlength. By way of a numerical example. 
Suppose we multiply two number. a = 0.2310 | 581 and b = 
0.8351 | 621 is a decimal coded machine using a precision of 
wordlength of four decimal digits. 


Then, in the first place we write down the intervals 
containing a and b. 


a € (0.2310, 0.2311] and b € [0.8351, 0.8352] 


Then carrying out the multiplication we have 
ab € [0.2310, 0.2311] [0.8351, 0.8352] 

= [0.1929 | 18, ..., 0.1929 | 92 ...] 

= [0.1929, 0.1930]. 


Thus, rounded interval arithmetic leads to the following 
interval which contains ab : [0.1929, 0.1930], the width of this 
interval being 0.1930 — 0.1929 = 0.0001. 


To obtain the product ab more sharply we clearly have to 
use a longer wordlength for the computation. 


The actual rounding down and up can most easily be 
achieved in machine coding. The lower endpoint should be 
rounded down so that the rounded value is less than or equal to 
the value being rounded. Thus if the value being rounded down 
is positive we need merely truncate while if the value is 
negative we truncate, examine the truncated portion, and, if this 
is not zero, we subtract one unit in the least significant position. 
And the upper end point must be rounded up if a computation is 
being carried out in rounded interval arithmetic, an interval 
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being said to be properly rounded if the endpoints are rounded 
as described above. 


We may mention that the rounding up and down can be 
coded in Fortan necessary although not as efficiently. 


Thus the final interval obtained by rounded interval 
arithmetic takes account of the errors due to rounding and 
contains the true value. 


Interval vectors and matrices and vectors and matrices 
whose elements are interval numbers, the superscript I being 
used to indicate such a vector or matrix. 


Thus given matrices B = (bj) and C = (cj) of order n such 
that 

bi < cj, i,j = 1, 2, ..., n, then the interval matrix A= [B, C] 
is defined by 

A' = [B, C] = {A = (aj) | by < ay < cy i, j = 1, 2, ..., np. 
Kuperman [26] calls B the lower endpoint matrix of the interval 
matrix A' = [B, C], C the upper endpoint matrix and M = 
Y2 (B+C) the midpoint matrix. 


DEFINITION 5.1: Let A = (aj, ..., d,) and B = (bj, ..., by) a;, 6; 
€ R and Z, or O or Z or C where A and B are of some order, 
viz; 1 xn. We define the row matrix of natural class of interval 
as 


[A, B] = [(ap, «1, An), (b1, «.., by]. 


If each a; > b;, 1 < i <n we say [A, B] such that A > B if 
each a; < b;, 1 < i < n we say the row matrix natural class of 
interval [A, B] with A < B. For more about properties and 
working with natural class of intervals refer [57-58]. 


If some aj > b; and some a; < b; and some a, = by we cannot 
say A>BorA<B. We will illustrate the row matrix of natural 
class of interval or in short row matrix interval by some 
examples. 
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Example 5.1: Let X = [(5, 3, —1, 0, 2, 4), (2, 4, 5, 7, -3, 1)] = 
[A, B] is the row matrix interval clearly A is not comparable 
with B. 


Example 5.2: Let Y =[A, B] = [(5, 2, 0, 7, 9, 6, 10), (7, 3, 4, 5, 
9, 10, 8, 12)] be a row matrix interval where A < B. 


Example 5.3: Let M = [(8, 3, 4, 9, -2, 7, 18, -9), (2, 1, 3, 5, — 
10, 2, 9, -14) = [A, B] be the row matrix interval where A > B. 


Example 5.4: Let S = [(3, 4, 2, 1, 5), (3, 4, 2, 1, 5)] = [A, B] be 
the row matrix interval, here A= B. So S = (3, 4, 2, 1, 5). We 
call this row matrix interval as degenerate interval. 


We will now proceed onto define the notion of addition, 
multiplication, subtraction and division of these row interval 
matrices. 


First we define the notion of order in a row interval matrix. 


Let X = [A, B] be a row interval matrix where order of A 
and B is 1 x n. We define the order of X, the row interval 
matrix to be 1 x n. Thus X is a 1 x n row interval matrix. 


We can add two row interval matrices if and only if they are 
of same order otherwise addition is not defined; for we cannot 
add a 1 x 5 row matrix with a 1 x 9 row matrix. 


DEFINITION 5.2: Let X = [A, B] = [(a1, a2, ..., An), (b1, by «.., 
b,)] and Y = [C, D] = [(c1, C2 «5 Cn), (di, Az «.., d,)] be any two 
1 xn row matrix interval. We define X + Y = [A, B] + [C, D] 

=[A+C,B+D] 

= [(A}, 0005 An) + (Ch voy Cn), (By «05 Dn) + (dy, «5 Ay] 

= [(a1 FCI, ase5 An F Cy); (b; + di, ach b, + d,)] 


X + Yis again a 1 xn row matrix interval. 


We will first illustrate this by some examples. 


Interval Matrices and Natural Class of Intervals | 159 


Example 5.5: Let X = [A, B] = [(3, 1, 5, 0, —3, 2), (7, 8, 9, 11, 
0, 1)] and Y = [C, D] = [(-5, 2, 3, -7, 8, 10), (2, 11, 10, -4, 5, - 
3)] be any two 1 x 6 row matrix intervals. 


We see X + Y = [A, B] + [C, D] 
10, 1) 


2, 11, 10, -4, 5, -3)] 
8, 10), (7, 8,9, 11, 0, 1) 


= [(3, 1, 5, 0, —3, 2), (7, 8, 9, 11, 

[(-5, 2, 3, -7, 8, 10, (- 

= [@, 1, 5, 0, -3, 2) + 5, 2, 3, -7, 
+ (-2, 11, 10, -4, 5, -3)] 

=[(3 + CS), 142, 5+3, 0+(—7), -3 + 8, 2+10), 
(7-2, 8+11, 9+10, 11+(-4) 0+5, 1 + (-3)) 

= [C2, 3, 8, -7, 5, 12), (5, 19, 19, 7, 5, -2)] 

is again a | x 6 row matrix interval. 


Example 5.6: Let P =[A, B] = [(3, 7, 1), (5, 1, -3)] and 


R = [C, D] = [(2, 1, 0), 7, 5, 2)] be two 1 x 3 row matrix 
intervals. 
P+R =[A, B]+[C, D] 
=[A+C,B+D] 
= [(, 7, 1) + C2, 1, 0), (5, 1, -3) + (7, 5, 2)]. 
Now as in case of natural class of intervals we may have 
A+B>C+DorA+B<C+D orA+B not comparable with 
C+DorA+B=C+D. 


Example 5.7: Let 
S = [(0, 0, 0, 0, 0, 0), (3, 1, 2, -5, 7, 2)= [S), S2] and 

= [(7, 2, 1, 0, 5, 1), (0, 0, 0, 0, 0, 0)] = [Pi, P2] be any two 
1 x 6 row matrix intervals. 


Now S+ P= [Si, S>] a [Pi, P>| 

= [(Si + Pi), (S2 + P2)] 

= [(0, 0, 0, 0, 0, 0) + (7, 2, 1, 0, 5, 1), (3, 1, 2, —5, 7, 2) 
+ (0, 0, 0, 0, 0, 0)] 

=([(7, 2, 1, 0, 5, 1), (3, 1, 2, -5, 7, 2)] 


is again a | x 6 row matrix interval. 
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We observe that addition of 1 x n row matrix intervals is 
commutative for if A = [A;, A2] and B = [B,, B2] then A + B = 
BA. 


Further addition of 1 x n row intervals is associative, that is 
ifA= [Ai, Aol, B= [Bi, Bo] and C = [Ci, Cy] then (A + B) +C 
=A+(B+C). 


We define 0 = [(0, 0, ..., 0), (0, 0, ..., 0)] as the zero 
interval matrix or 1 x n row zero matrix interval or 1 x n row 
matrix zero interval. 0 acts as the additive identity. 


For if A = [A;, A>] and (0) = [(0), (0)] then A + (0) = 
(0) + A=A for all A; 1 x n row matrix interval. Further if 
A = [Aj, Ag] then —A = [—A;, —Ag] and A + (—A) = [(0), (0)]. 


Let V = {[A, B] = M where A and B are 1 x n row matrices 
with entries from R or Z or Q or Z, or C}; V is an additive 
abelian group of | x n row matrix interval. 


Now we can define product of two 1 x n row matrix 
intervals. Suppose A = [A;, B,| and B = [C;, D,] be two 1 x n 
row matrix intervals; we define product AB as follows: 


Let AB = [Ai, Bi] [Ci, Dj] 

=[AiC;, BiD,] 

= [(ai, A2, «+5 an)s (c1, C2, wey Cn) (b,, by, sibs b,), 
(dj, do, ..., da)] 

= [(ai, Cj, a2C2, ..., AnCn); (bd), bod>, tales b,d,)]. 


We see AB is again a 1 x n row matrix interval. 
We will illustrate this situation by some simple examples. 
Example 5.8: Let 
V= [A, B] = IG; 7; 2, 0, af 8, 9, 11), (3, 1, +2, 4, 0, 18, 3, 0)] 
1 


and S= [Si, S>] > [(3, 2, 0, l, 5; 0, -2, 0), (0, 4, 0, 2, > 0, 2; —4)] 
be two | x 8 row matrix intervals. 
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Consider VS = [A, B] [Sj, S2] 
= [AS;, BS] 


=[6, 7, -2, 0, 7, 8, 9, 11), (3, 1, +2, 4, 0, 18, 3, 0), 
(3,2 > 0, 1, 35 0, a2, 0), (0, 4, 0, 2, l, 0, 2, —4)] 
=[(5, 14, 0, 0, 35, 0, -18, 0), (0, 4, 0, 8, 0, 0, 6, 0)] 


is again a | x 8 row matrix interval. 


Example 5.9: Let V = [(3, 2, 0, 5, 1), (7, 8, 1, 0, 3)] and 
= [(0, 1, 2, 0, 5), (3, 0, 0, 7, 0)] be two 1 x 5 row matrix 
intervals. Now VM = [(3, 2, 0, 5, 1), (7, 8, 1, 0, 3), (0, 1, 2, 0, 
5), (3, 0, 0, 7, 0)] 
= [(0, 2, 0, 0, 5), (21, 0, 0, 0, 0)] is again a 1 x 5 row matrix 
interval. 


Example 5.10: Let A = [(-3, 2, 5, 1, 0, 2), (7, 0, 0, 8, 0, 1)] and 
= [(0, 0, 0, 0, 8, 0), (0, 7, 5, 0, —3, 0)] be two 1 x 6 row matrix 
intervals. 


We see AB = [(-3, 2, 5, 1, 0, 2), (7, 0, 0, 8, 0, 1)] 
[(0, 0, 0, 0, 8, 0), (0, 7, 5, 0, —3, 0)] 
aa [(-3, a, Dy 1, 0, 2); (0, 0, 0, 0, 8 0), (7, 0, 0, 8, 0, 1) 


(0, 7, 5, 0, -3 
= [(0, 0, 0, 0, 0, 0), (0, 0, 0, 0, 0, 0)]. 


Thus A # [(0), (0)] and B # [(0), (0)], 
but AB = [(0), (0)]. 


Now we see if V = {M = [A, B] | collection of all 1 x n row 
matrix intervals} then V is closed with respect to multiplication. 
The operation multiplication is commutative and associative on 
V; but we see V is only a commutative semigroup under 
multiplication. For we see V has infinitely many zero divisors. 


Example 5.11: Let M = [A, (0)] and [(0), B] = P be two 1 x n 
row matrix intervals. We see MP = [(0), (0)] is a zero divisor. 


Now we define subtraction of two 1 x n row matrix interval. 
Suppose M = [(5, 3, —1, 0, 5), (3, 2, 4, 2, 0)] and P = [(3, 0, 1, 2, 
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4), (2, -1, 5, 4, -3)] be two 1 x 5 row matrix intervals. We fine 
M — P= [(5, 3, -1, 0, 5), (3, 2, 4, 2, 0)] [(3, 0, 1, 2, 4), (2, -1, 5, 
4, -3)] =[(Q, 3, -2, -2, 1), (1, 3, -1, —2, 3)] is again a 1 x 5 row 
matrix intervals. 


The following points are pertinent to be observed. 
(1) The operation subtraction is not commutative. 
(2) The operation of subtraction is non associative. 


Consider A = [Aj, Bi] = [(3, 2, —1, 0), (8, 5, 1, 2)] and 
B = [A), By] = [(6, 3, 4, 8), (4, 3, 2, 7)] two 1 x 4 row matrix 
intervals. 


A—B=[A,, Bi] — [A2, Bz] = [Ai — A», By — Bo] 
= [(3, 2; 1, 0) (6, 3, 4, 8), (8, 5, i 2), (4, 3, 2; 7)] 
= [(-3, -1, -5, -8), (12, 2, -1, —5)]. 


Consider B — A = [Ap, Bz] — [Aj, Bi] 

= [96, 3, 4, 8), (-4, 3, 2, 7)] — [(3, 2, -1, 0), (-4, 3, 2, 7) 
— (8, 5, 1, 2)] 

= [(3, 1, 5, 8), (-12, —2, 1, 5)]. 

A-B#B-A. 


Thus the operation subtraction is non commutative. 


It is easy to check A — (B — C) = (A- B) —C for any | x n 
row matrix interval. 


Consider A = (Aj, B,), B = (A:, Bz) and C = (As, B3) be 
three | x 3 row matrix interval. 


AS (Ay, Bi) = (6,2, 0) 4-2, 159), 
B= (Ag, B) or 3, 2s 4), (3; E2 =1)] 
and C = (A3, Bs) = [(2, 4, 5), (0, 1, 2)]; 


(A -B)-C=([@, 2, 0), (2, 1, 5)]-[C3, -2, 4), (5, 3, -D] 
—[(2, 4, 5), (0, 1, 2)] 
)}- 2, 4, 5), 0, 1, 2)] 
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A-(B-—C)=[(3, 2, 0), (-2, 1, 5)] — (((3, -2, 4), (5, 3, -1)] 
“7 [(2, 4, 2) (0, l, 2))) 

= [(3, 2, 0), (-2, 1, 5)] — [-5, -6, -1), (5, 2, -3)] 

= [(8, 8, 1), 7, -1, 8)]. 

Thus (A -B)-C #A—(B—C). 


We now proceed onto define the notion of division of two 
1 x n row matrix intervals. 


Division of an 1 x n row matrix interval A = [A,, B;] can be 
divided by B = [A), B,] if and only if no coordinate in A, and B, 
is zero: 


That is if A = [Ai, B,] = [(ai, a2, .++5 an) (b;, b», sey b,)] 
and B = [(c), C2, ..., Cn), (di, da, ..., da)]. 


Then A/ B= [A\/A>, B,/Bo] 
= [(aj/C1, a2/C2, ..., An/Cn), (b1/d), b2/do, ..., b»/dn)] where c; # 
0 and dj; #0 for 1 <i<n. 


We see clearly A/B # B/ A even if none of the coordinates 
in A, and B, is zero. 


Also the division operation is non associative. 


Now we can give better structure for the 1 x n row matrix 
interval. 


THEOREM 5.1: Let V = {M = [A, B] the collection of all I xn 
row matrix intervals with entries from R or Q or Z or Z, or C or 
C(Z,)}; Vis a commutative ring with unit and has zero divisors. 


The proof is direct and hence is left as an exercise to the 
reader. 


Note V in general is an infinite ring. V is a finite ring only 
the entries of M are from Z, or C(Z,). 
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We will give one or two examples of this ring. 


Example 5.12: Let 

M={P = [A, B] = [(a1, a2, a3), (bi, be, b3)] | ai, bi € Z25 1 Si < 3} 
be the collection of all 1 x 3 row matrix intervals. Clearly M is a 
finite commutative ring with zero divisors. Further |M| = 64. 


Example 5.13: Let 

P= {M=[A, B] = [(ai, a2), (bi, b2)] where a;, b; € Z, 1 <1 < 2} 
be the collection of all 1 x 2 row matrix intervals P is a ring. P 
has no units but zero divisors. P is of infinite order. 


Example 5.14: Let K = {M = [A, B] = [(ai, ..., aio), (bi, «-., 
bio)] | ai bj € C, 1 <i < 10} be a collection of all 1 x 10 row 
matrix intervals, K is a ring of infinite order. 

We can define substructures on these structures. This is 
simple and hence is left as an exercise; these rings have subrings 


and ideals. 


Now we bring in the representation / connection between 
1 x n row matrix intervals and | x n interval matrices. 


We know N = ({ay, by], ..., [an, bn]) is a 1 x n row interval 
matrix and M = [A, B] 


= [(a, ..., an), (bi, ..., b,)] is a | x n row matrix interval. 
We see every M can be made into N and vice verse. 
Consider N = ([aj, bi], ..., [an, bn]) = [(a1, ---, an), (bi, .--5 Dn)] 


(by taking the first components of every interval in N together 
and the second components of every interval in N together. 


Consider M = [(a), ..., an), (bi, ..., Da)] = ((ar, bi], ..., [any Dn]). 
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Write M as n intervals by taking the first component of A 
and first component of B to form the first interval [a;, b;] and so 
on we see M = [A, B] becomes N and N becomes M. 


Now we proceed onto define the notion of n x 1 column 
matrix interval. 


DEFINITION 5.3: Let 


a, b, 

2 || 

M=/[A, B] = oe 
a b 


be a pair of column matrices. We define M to beam x 1 
column matrix interval. 


From now on wards if we say a m x | column matrix 
interval we mean a pair [A, B] where A and B are m x 1 column 
matrices (vectors). 


We will illustrate this situation by an example. 


Example 5.15: Let 


ND Wn BW NY 


Misa7 x 1 column matrix interval. 
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Example 5.16: Let 


io) 
[oe] 


N 


& 


be a5 x 1 column matrix interval. 


Example 5.17: Let 


be a 6 x 1 column matrix interval. 


We can add only two column matrix intervals of same 
order; just as in case of usual m x | column matrix interval. 


We will just show how addition of two m x 1 column 
matrix intervals is obtained. 
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Example 5.18: Let 


8 | | -2 2/|-3 
2 0 1}| 2 
0 5/|-l 
M=]| 4 |,] 5 and N= | | 4},| -4 
5 ||-1 0} | 3 
7 8 3] | 2 
-l}| 2 5/| 4 


be two 7 x | column matrix intervals. 


8]/-2]) /((2|[-3] 
2 || 0 1}| 2 
Ou. 5||-1 
We findM+N=]|| 4 |.) 5 |] +||4|,) 4 
5 ||-1 0|| 3 
7|| 8 eM ie 
=1]} 24) Us|. 4] 


lI 
NUON 

+ 
nw oO fh NY NY 

n 

+ 

| 

& 
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10 


10 


aN 


is again a 7 x 1 column matrix interval. 


Example 5.19: Let 


Il 
SCO wm N a 
nA BR wWN 

teh) 

5 

a. 


5 0 
4 }/-l 
P=) -—3},| 2 
2 0 
-1}| 3 


be any two 5 x | column matrix intervals. 


6] [1 

Fi 2 
M+P=//8],] 3 
9||4 

0/|5 

Ke 5 

7 4 

=|} 8 }+]-3 

9 2 

|O| [-1 


5 |[ 0 
4 ||-1 
+}|-3],| 2 
21] 0 
|-1] | 3 
1 0 | 
2\" (=1 
J 3}+] 2 
4} |0 
| 5 3 
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yf 
1} fa 
=|| 5,15 
11] } 4 
-1] | 8 


is a5 x | column matrix interval. 


We see addition of two n x | column matrix intervals is 
commutative and associative. 


0||0 
O=]] .|,| . || 1s the zero column matrix interval. 
0} |0 
Further for every 
a, b, —a, —b, 
a b —a —b 
M= 2 || 2 || we have —M = a eas 
a, b, —a, —b, 


0} |0 
and M+ (—M)=]] . |,} . 


Inview of all this we have the following theorem. 
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THEOREM 5.2: Let 
={M=|| 2 || 2 || a, b; €ZorR or C(Z,) or QO or Z, or C} 


be the collection of all n x I column matrix intervals. V is an 
abelian group under addition. 


The proof is simple and hence is left as an exercise for the 
reader. 


Clearly if M and N be two n x | column matrix intervals 
then MN cannot be defined or found. 


Example 5.20: Let 


a, b, 

a, | | b, ; 
M=;sP= : where aj, b; € Z3, 1 <i < 4} 

a, b, 

ag. || By 


be the collection of all 4 x 1 column matrix interval. M is an 
abelian group under addition of finite order. 


Example 5.21: Let 


a, b, 

as! |B, : 
W=4S= «Isl where a;, b; € Z, 1 <i <7} 

ag | |e 


be a 7 x 1 column matrix interval group of infinite order. 
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Now having seen examples of m x 1 column matrix interval 
group. We now proceed onto define m x n matrix intervals 
(m #n). 


We will illustrate this by some examples. 


ai, a, 11 b,, 

a b b 

21 2 21 2 
Let P = man Fs 
ani Finn bi bi 


where aj, b;, € Z or Q or R or C or C (Z,) or Zy. 1 $1 < mand 1 
<j<n}. We say P is am x n matrix interval. 


We will illustrate this by some examples before we define 
some operations of these sets. 


3 01 4//-3 5 1 2 
5 3 4 0/;/0 1 2 O 
LetP=|}}-1 2 3 4},) 1 2 0 44), 
7 0 8 4}|5 1 6 2 
0 10 2};0 01 -!1 


Pisa5 x 4 matrix interval. 


We can as in case of row / column matrix intervals add any m x 
n matrix interval. 


Consider 


and 


NON W 
—- fb OF 
| 
a 
i) 
| 
— 
orn s 
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—2 1 4/;/ 1 2 3 
0 3 l}/-4 5 -6 
NS ; 
2 0 4//7 -8 -9 
5 1 0]/|0 1 -2 


be two 5 x 5 matrix intervals. We find M+ WN as follows: 


3 4 olf2 8 4 
men=|{2 © 2]P9 24h, 
64 3h) 0.4 
71 0||6 4 0 
BBS Gs Ad [iis eh 
Oe ae See 
2 0 4/|7 -8 -9 
rollo 4 
54-0) 72.1 42 8 4) Tr 3 
FO at be Bee Hts 24S. 26 
te) ee ee cee 
F-41011 SL OV A OOM NO OTS Se 
i Sa ee | 
Bem ee 
Wiha sae IE) Be Ses 96 
1 30 6.52 8 


is again a 4 x 3 matrix interval. Thus we can add matrix 
intervals of same order. 


In view of this we have the following theorem. 
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THEOREM 3.3: Let V = {[[@ij] mv [bij] ma] where aj, by € Z, or 
C(Z,) or Zor OQ or RorC, 1 Sism,1 Sj Sn} be them xn 
matrix interval. V is an additive abelian group. 


If elements of V are in Z,, or C(Z,) then V is finite otherwise 
V is infinite. 


However multiplication is not defined on V _ as 
multiplication is not compatible. 


Example 5.22: Let 


a, a,\(b, b, 
V= a, a, |, b, b, ai, bj € Zs, 1 <i,j <6} 
a; as b, b, 


be the additive 3 x 2 matrix interval. V is a group of finite order 
and V is commutative. 


Example 5.23: Let 
Mz Ce: a: oi ibe Bs. “Dy. By 
Mae Ae. Hay. Ae fh De: Da? Be. D; 
ai, bj € Z, 1 i,j < 8} 


be the additive 2 x 8 matrix group interval. M is of finite order 
and is commutative. 


Example 5.24: Let 


36 de 26 3 


be two 2 x 2 matrix intervals. 
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We define 


ere HE Me Se 

sone ea ttle) 
mf e Me He OD 
Jt . np . 4) 
Ne 


aa,+be, ab,+bd,)/ee,+fg, ef, +h, 
ca,+dc, cb,+dd,)\ge,+hg, gf, +hh, 


is again a 2 x 2 matrix interval. 
a By Ve Fy a pipe sk 
LP = ; , 
ec di) \g, bh, e dj\g h 
_ |({aa,t+bc ab+bd)\fae+bg af+bh 
catde be,+dd)’lcet+dg cf+dh}|- 


Clearly PT # TP. 


Thus the product of two square matrix intervals is in general 


non commutative. 


Further we see the product of square matrix intervals can 


also be zero. 
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Example 5.25: Let 


be any two 3 x 3 matrix interval. 


| aS || 
OL mt © 
om N = Oo 
— Oo /O om 
ere SSS; 
non CN ea: ie 
oOo - © oN CO 
| aS SAN 
————— ay | a 
Aas Qon 
oN CO nan a 
-=—= ON on Oo 
Tes Fen 
omN a 
nano nNOn oO 
_—ESSSs 
I I 


is a3 x 3 matrix interval. 


Consider 
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is again a 3 x 3 matrix interval. However RP # RP. 


Example 5.26: Let 


be two 4 x 4 matrix intervals. 


is again a 4 x 4 matrix interval. 
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Consider MN = 
0 2 1 0)f-1 0 2 1\)(/4 0 1 0/0 1 0 ~+1 
4001} 0 1 0 2//0 2 0 442 0 1 0 
102 OFF 1 1 0 07;7}1 1:0 OF} 1 0 2 0 
03 4 0/0 10 i1/(0 0 1 OF0 2 0 +1 
1 3 0 4\)(1 4 2 4 
_|{-4 1 8 5)/4 8 10 4 
L! 2 a A 
4 7 0 6/\1 0 2 0 
is a4 x 4 matrix interval. 
Now NM = 
-1 0 2 1\/0 2 1 0)\/0 1 0 1)f4 0 1 +0 
0 10 214 00 1//2 0 1 0/0 2 0 4 
Lea OO a 8. 2 OUT! 2, eG VG 
0 1 0 1/0 3 4 0/10 2 0 1/0 0 1 0 
217 0)/0 2 1 4 
2 aG B Dp Qeod) 22° 0 
4211/16 21 0 
43 4 1)/\0 4 1 8 


is a4 x 4 matrix interval. 
We see MN # NM. We have the following theorem. 


THEOREM 5.4: Let V = [[dij] mon [bij]mxn] be a collection of m 
xm interval matrices with aj, bj € Z, or C(Z,) or Z or Q or R 
or C. V is a non commutative ring with respect to matrix 
interval addition and multiplication. 
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This proof is direct and simple and hence left as an exercise 
for the reader to prove. 


Clearly V is a ring with zero divisors units and idempotents. 

Further we see as in case of row matrix interval we can in 
case of square matrix interval also obtain the square interval 
matrix. Also from a square interval matrix find the square 
matrix interval. 


This we will illustrate by some examples. 


Example 5.27: Let 


52.0 i Qs 5: Tb 32 
y-|[3 5 4 1]/0 15 -1 5 
2: 3 Oe Aa OF Ae 0 
016 0J/\0 1 0 7 


be a 4 x 4 matrix interval. We have V = [Vi, V2] = [(aij), (bi)]. 
Now we can write V = ((aj, bij]). 


[5,3] [0,-5] [Ll] [2,2] 
[3,0] [5,15] [4-1] [15] 
[2,1] [3,2] [0,3] [4,0] 
[0,0] [11] [6,0] [0,7] 


[5,3] [0,-5] 1 2 

[3,0] [5,15] [4-1] [1,5] 

[2,1] [3,2] [0,3] [4,0] 
0 1 [6,0] [0,7] 


is clearly a 4 x 4 interval matrix where intervals in V are 
increasing or decreasing or degenerate intervals. 
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Thus given a 4 x 4 matrix interval we can convert it into a 
4 x 4 interval matrix. 


Now we just show how an interval square matrix M can be 
written as a square interval. Consider 


[0,1] [2,3] [5,0] [1.1] [2,-1] 
[4,2] [0,0] [7,1] [8,0] [0,5] 
M=|] [11] [2,4] [3.3] [L-1] [8,9] 
[7,1] [4,3] [1.2] [10,0] [4,4] 
18,5] [7,7] [2,0] [5,5] [3,7] 


be a5 x 5 interval matrix. 


Or 2 Se a DT, BO. tS 
A he Fe Be OH Oe Oe Ae > 
WriteM=||/1 2 3 1 8|J1 4 3 -1 9 
7 ASE 40 AS Be Oe A 
RG Deh BS &. So °F 


= [M,, Mo]; M is a5 x 5 matrix interval. Thus we have the 
following theorem the proof of which is direct and simple. 


THEOREM 3.5: Let M = ([mj, nj]) be an xn interval matrix. 
Then M = [(mj), (nj)] is an xn matrix interval. Conversely if A 
= [(aj), (by)] is an xn matrix interval then A = ([aj, bj]) is a 
n Xn interval matrix. 


The proof is left as an exercise to the reader. 


Thus we see we can have in case of square matrix intervals 
the notion of determinant of intervals. If M = [M;, M2] where 
M, and Mz; are n x n matrices then the determinant interval M is 
denoted by |M| = |[Mi, Mo]| = [|Mj], [Mal]. 


180 | Innovative Uses of Matrices 


We will first illustrate this situation by some simple 
examples. 


Example 5.28: Let 


be a 2 x 2 matrix interval. 


Determinant interval of P denoted by 


mG 3) 


=[1, 30] is again an interval matrix. 


Example 5.29: Let 


1S’) 
~ 
on 
So 


N 
So 
Nn 
=) 
a 


be a3 x 3 matrix interval. 


|A| = ||] 0 


Nn 
= 
_ 
So 
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1 0 0 0 0 1 
7 —0 +1 
0 -l 5 -l 5 0 
= [61, —12] is the interval determinant. 


Now we show how inverse of an interval matrix can be 
found. 


Suppose M = [M,, Mp] is a n x n matrix interval, then 
M'=[M,',M;' J. 
Thus MM‘ = | M,M,',M,M;"' | 
= [M,"M,,M;'M, | 
=M'M 
= [Is Tal: 


We can use any known method to determine the inverse of a 
n X n matrix interval. 


We can use Gauss - Jordan method to find inverse for a 
matrix interval. 
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Suppose 
aj, AyD Ain b, b, b,, 
— | far Aa a,, | |b, b, b,, 
M _ < ny é . 
an ano Aon b,, Dis Bea 


= [A, B] = [(aij), (b)] then the rows of A are denoted by Ri, 
Rp, ..-, Rin and the rows of B are denoted by Ro, Ra, ..., Ron 
we write M as 


ay he ae Pm 
Gigi Ags a,,|90 1 0 
ai sa is ain|9 0 1 
by Bye b,, | 1 0 .. 0 


We by row operations on A and B make A and B to Inxn 
identity matrices and all the operations are done on A and B are 
also simultaneously performed on Inyn. 


The adjoined matrices Inx, turns to be A! and B! 
respectively. That is M'=[A', B']. 


For instance if 


M =[A, B] = [(ay), (bi)] 
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2 
i to find M?. 
3 


i 


-1/3 1/3 


1/3 


0 
Now M!= | 


—1/3 "|| 


1\( 1/3 
2 || 0 
1)}\ 0 


3 1 
0 1 
0 0 


MM! = 
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ll 
oo Fr 
or Oo 
r- CO CO 
oo Fr 
or Oo 

i=) 


= [, L]. 
Also 
Ls = TB VSL PV. ae? SV Se 
MM) Oe SNe ae Bh, SBP: ae Oe tS 
0 0 10 0 1J/lo 172 o}lo 1 3 


Thus MM!=M'!M. Mis the inverse matrix interval of 
M. 


We can using Gauss Jordan method to find the inverse 
matrix interval. 


Now we can use the method of row reduction and bring a 
m X n matrix interval M to row-reduced echelon matrix interval. 


(i) If the m x n matrix interval M = [M,, Ni] is row 
reduced 


(a) If the first non zero entry in each non zero row of 
M, and N; is equal to | in M. 


(b) Each column of M, and N,; in M which contains 
the leading non zero entry of some row has all its 
other entries 0. 


(ii) Every row of M, and N; which has all its entries 0 
occurs below every row which has a non zero entry. 
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(iii) If rows Ru, Rp, ows R 


non zero rows of M, and N, respectively and if the 
leading non zero entry of rows p and q occurs in k, 
and t, respectively of M; and N;; p= 1, 2, ..., s; and 
q- 1, 2; aise: SD, with k, < ky Sa Ss k, and th <t< 


and Ro, R», wishes R,,, are 


Is, 


sats 


We can also describe an m x n row reduced echelon matrix 
interval M = (M,, N;) as follows: 


Here M = [(mj), (nij)] = [Mi, Ni]. 


Either every entry in M; and N, is 0 or there exists positive 
integers Ss}, S2; 1 < s}, Ss < mand s; and s) positive integers. ky, 
k, 5 and ti, to, ..., t,, with 1 < k, tj <nand 

(a) mj=0andn,,=0 ifi>s, and mj =0 ifj <kj and 
p > Ss and n,, = 0 ifv < ty. 


(b) m6 18155, 15j755 n 


l<v<sgp. 


= Spy, 1 Sp S89; 


(c) ki <i <4. k, andth, b= ...~<t,. 


We just give examples of 3 x 5 row reduced echelon matrix 
interval R = [Rj, Ro] 


0 1 2. 0 U3) (0: 1 —s- 0 V7 
=|}0 0 0 1 54);0 0 0 1 2 
00 0 0 O 00 0 0 =O 


The following theorem is simple and hence is left as an 
exercise to the reader. 


THEOREM 5.6: Every m xn matrix interval M = (M,, N;) is 
row equivalent to a row reduced echelon m xn matrix interval. 
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P(x) is the polynomial intervals we have introduced the 
notion of interval polynomials in [55] and have discussed their 
related properties. 


Suppose we have two sets of m-linear equations in n 
unknowns. 


say Ai X1 + Aq X2 +... + Ain Xn = Y1 
Aa X1 + Ago Xo +... + Aon Xa = V2 


AmiX1 + Am2X2 Fis oe AmnXn = Ym 


where yy2 ... Ym and Aj; 1 <i<m, 1 <j <nare given elements 
of a field F. 


Any n tuple (x;, ..., Xn) of elements of F which satisfies 
each of the equations is a solution of the systems. 


Consider another set of m-linear equations in n-unknowns. 


say Bi x + Bj >a Tzes +Bin < = Z\ 
Bo x, + By a +...+ Bon Xi =Z) 
Bini X, Bye Xo Ft Baek, = Ze 


where Z, Z2 ,..., Zm and By; 1 <i <m, 1 <j <n are given 
elements of F. 


Any n tuple (xj, ..., x!) of elements of F which satisfies 
the equations is called the solution of the systems. 


If y) = yo = ... = Ym = O and z; = Zz . = Zm = 0 Is 
homogeneous set of equations or that each of the equations is 
homogeneous. 


We now define the notion of interval system of m linear 
equations in n interval unknowns. 
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Consider 
t ' 
A,X, +. +A Xn B,.X, +--+ BX, y, Z, 
t t 
A,X, +..+A,,X%, BX, +..+B,,x, _ || ¥2 | | 22 
. . . . . | . 
t ' 
AX tet Ag Xns BuiX +--+ Ban Xn Yin | | Zim 


(5.1) 


(5.1) is a system of m linear interval equations in interval n 
. . ! Uy 
unknowns. A row matrix interval [(X1, X2, ..., Xn), (X},X55 -65 


x’ )] of elements from the field F which satisfies each of 
equations (5.1) is called the solution interval of the system. 


We will first give some examples. 


Bh = Ket Ry) SR Ry Te 
Consider Ry hoRy xX, — 4x, 


=4x, + 2%, =5XK, 7x, —2x, 


7 || 2 
=]| 3 || 0 
—1 || -8 


is an 3-linear interval equations in 3 unknowns. 


We give yet another example. 


Take 
X, —3X, +x, —4x, +x, —3x, +4x, —X, 4/) 1 
8X, +X, +X, — 5X, 2x, —-4x,+x%,+2x,|_ ||2)] 2 
Is FOR, HAR Eke . Kok, HR OM, 3/| 3 
x, + 4x, X, +X, +X; OF [2 
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this is again a 4-linear interval equations in 5 unknowns. 


Yet we consider another example. 


OX) = Ry PR ELK, x,+3x,-4x, | 3 0 
3x, +X, —2x, 2X ate -1]| 2 
5x, —X, —3x, 5 te es ee =|} 2 |) 1 
1%). PXy + 8m, Jkt RAR 4 

2X) +k, 2K) SX, —3X,;+2x, | 0) \-5 


is again a 5-linear interval equations in four unknowns. 


Now we can give these linear interval equations a matrix 
interval equations representations which is as follows. 


Ai, Ay Aun Bi, B, B., 
Ay, Ay Aon B,, B,, B,, 
Ant Ais 7 er Bu B. Be 


([(x1, seey Xn)s (eae sr Xi ny) 


yi Zy 
_ || Y2 Z, 
= oli 

Ym | 2m 


Here [A, B] [X, X"]=[Y, Z] where A = (Aj), B = (B;) 
X= (Xj, ..-, Xn) and X’ =(xj,..., X,)- 


We call [A, B] the matrix interval of coefficients of the 
system. 
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Thus this m x n matrix interval with entries from the field F 
is a function [A, B] from the set of pairs of integers (i, j); 1 <j < 
m; 1 <j <n, into the field F. The entries of the matrix interval 
[A, B] are scalar intervals [AQ, j), BG, j)] = [Ai Bij] and quite 
often it is convenient to describe this matrix interval by 
displaying it entries in a rectangular interval array having m row 
intervals and n column intervals. [X, X’] is an x 1 matrix 
interval and [Y, Z] is am x 1| matrix interval. 


Thus [A, B] [X, X’] = [Y, Z] or [AX, BX’] = [Y, Z] is 
nothing more than a shorthand notation of the system of linear 
interval equation. 


Just like in usual matrices we can define row equivalent 
m xX n matrix interval [A, B], [A’, B’] such that [AX, BX’] = 
[(0), (O)] and [A’X, BX’] = [(0), (0)] have exactly the same 
interval solutions (solution intervals). 


The proof of this result is direct and can be derived by any 
interested reader. 


We can as in case of usual simple matrices say an x n 
matrix interval M = [A, B] is an invertible n x n matrix interval 
and N = [C, D] is another n x n invertible matrix interval then 


(1) M'=[A,B]'=[A‘,B’] 
and(M')'=[A' BJ! 
=[A, B] (since (A')' =A. 
and (B')'=B). 


If both M = [A, B] and N = [C, D] are invertible n x n 
matrix interval then so is MN and (MN) '=N'M". 


Further it can be easily verified that product of n x n 
invertible matrix intervals is invertible. Further if A = [Aj, Bj] 
is invertible matrix interval; that is A row equivalent with 
I = [I,, I,] the n x n identity matrix interval. 
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If A = [A;, B;] is an n x n invertible matrix interval then 
AX = [(0), (0)]; that is [A;X;, B;X.] = [(0), (0)] has only the 
trivial solution X = [X;, X2] = [(0), (0)]. 

The system of interval equations AX = Y 

[Ai, Bi] [X1, X2] = [Y1, Yo] that is 

[A;,X), B}X2] =[Y1, Y2] has a solution 

X = [X1, X2] for each n x | matrix interval Y = [Y,, Y2]. 


We can define as in case of usual square matrices define 
these notions in case of square matrix intervals. 


Let A =[A,, B,] be an x n matrix interval over the field F. 


The principal interval minors or minor intervals of 
A = [Aj, Bi] = (Ai), (Bij)]; 1 <1, j <n are scalars 


Ay(A) = [Ax(A1), Ax(B1)] defined by 
Ax(A) = [Ax Ai, AxA2] 


where 1 <k <n. 


| eae By tas Big 


Let A = [A,, A)] be an x n matrix interval we say A is an x 
n upper triangular matrix interval if both A; and A, are upper 
triangular matrix interval. 


Consider M = [M,, Ni] 


Interval Matrices and Natural Class of Intervals | 191 


M,, M,, “ie M,, Ni Nip Ni eed Nin | 

Os “Me. ee Mee 0: NG NGS, <r ING. 
iO > MOP cts MN Oe oO” “<M, N,. 

O° 0). -22 3M. 0) 0" 0 N,, 
[0 0 M,,/( 90 0 0 Nun) | 


5 7 2 0)\f2 10 4 
. 0 1 -5 -7//0 3 1 -2 
that is M = : 
00 3 24/0 0 4 =O 
00 0 1 00 0 


is a4 x 4 upper triangular matrix interval. 


Take A = [A), B;] be an x n square matrix interval. If both 
A, and B, are n x n lower triangular matrices then we define A 
to be an x n lower triangular matrix interval. 


We will just give an example before we proceed onto define 
further structures. 


Let 
BD 0°. LTT a Oe DY) 
OF: Oo 6263.) OOk “O° SD 
PH (Pi, P= 0-0 5 DO Os D2 
0006 -1//0 00 1 4 
600 60 8)\0 00 0 7 


is a5 x 5 upper triangular matrix interval. 


Consider P’ = [P,, P2]'=[B,",P, ] 
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5.0 0) 107 SO Or “OO 30 
2° TO 0 OS. LO, 2028 
= 10 0S Or OU Ao 2. OO, 
062 6 O0}|/0 5 2 1 0 
he ie tee ee ee ee 
we see P’ is a lower triangular matrix interval. 
Consider A = [Aj, By] 
6 0 0 0]/2 0 0 0 
ah ae. Oe OO: Oh 
GAO ind: 
OFS i aie AL eS 
be a 4 x 4 lower triangular matrix interval. 
6-7 8. 0)f2 0-1 3 
‘ aa 0; Re oo Oe Sa 2. 0 
A’=[A,,B,] = ; 
0041/0 07 1 
00 0 2)\0 00 5 


is a 4 x 4 lower triangular matrix interval. 


Thus we see by taking transpose we can easily covert a 
n x n lower triangular matrix interval into a n x n upper 
triangular matrix interval and vice versa. 


We can easily verify the following statements are equivalent 
in case of an invertible n x n matrix interval A = [A,, Ao]; 


(a) There is an upper triangular matrix interval 
P= (Pi, Pa] = |(Pj),(P;)| with Py = Pa = 1 1<k 


ij 
<n such that B = AP = [Ai Pi, A>P>] = [Bi, Bo] isa 
lower triangular matrix interval. 
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(b) The principal minors of A = [Aj, A.]; that is the 
principal minors of A; and A, are all different from 0. 


We call A = [A, A] to be an x n diagonal matrix interval 
if A= [Ai, Ad] = [(Ajj), (Bii)] than Ajj = 0 and Bi =Oifiz j 
l<i,j<n. 


700 0/f8 0 0 0 
080 0llo 6 0 0 
Thus {Pi,PoI=|) 9 9 1 oflo o -1 0 
000 5/100 0 -7 


is a 4 x 4 diagonal matrix interval. 


Further it can be easily verified if A = [A;, Ao] = [(aij), (bi)] 
be an x n square upper triangular matrix interval that is if aj = 
0, bj = 0 for i >j, that is every entry below the main diagonal is 
zero. A is invertible if and only if every entry on the main 
diagonal of both A, and A, are different from zero. 


For instance 


oocmUcCOlUlUMOO 
oo co Oo 
oN MN sf 
oor N 
oO ON OO 


-7 


the upper triangular matrix interval is not invertible. 
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Prove if P = [P;, P2] 


1 1/2 1/3 1/4 1 1/2 1/3 1/4 
L2 TFS UA TS 1 Sy? Tes: 14 LS 
1/3 1/4 1/5 1/6}}1/3 1/4 1/5 1/6 
Lit 1S) TG T/T IA WS 1/6 ee 


is an invertible matrix interval and P' = [P,',P,'] has integer 
entries. 


We can define permutation matrix interval as in case of 
usual matrices. 


We will give an example or two. 


P = [Pi, Po] a 


oreo°7e 
oo or 
oor eo 
Ke COO ef 
oor fo 
oro oe 


is a permutation matrix interval. 


Now we proceed onto give example of symmetric and skew 
symmetric matrix intervals. 


Consider 
3 -1 2 44/8 4 3 1 
-1 8 4 0|/4 -1 4 2 
H = [Hi, Hp] = > 
2, 3 91/13 4 0 7 
4 0 9 2/}/1 2 7 2 


we see in this matrix interval aj = aj if14#j, 1 <1,j <4. 


Thus if A = [Aj, Bi] = [(aj), (bj)] be a square matrix 
interval. If in Ay and B,, aii = Aji and bi = bii if i #], l< i, j <n 


—_ 
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then we define A = [Aj, B,] to be a symmetric matrix interval. 
If A = [Aj, Bi] = [(aj), (b)] be a square matrix interval. 


If aij = Aji and bij = bj for 1 #] and if 1 =j then aii = bi = 0, 
1 <i<n,1<i,j <n, then we define A to be a skew symmetric 
matrix interval. 


Let P = [P,, Po] 


OB he EAS TS OREO. AW. OR. Br “a 
3° D2. ff ER OVE oO E62 Sp, Se <0 
_|J/-1 -2 0 -7 2 9//2 6 0 1 -2 3 
Nias sabe Se. Be abet oe ca. 24 5 
Oe ee ee ee ee 
28 GY AO: 2 SRO bes Se. 6, 0 


be a 6 x 6 skew symmetric matrix interval. 
We can as in case of usual square matrices write the square 
matrix interval as a sum of a square symmetric matrix interval 


and a square skew symmetric matrix interval. 


We will illustrate this situation by an example or two. 


Let 
8 4 3 5 2 -4 3 #1 
7 1 2 -2 0 7 2 
A= [Au, B,] =a 5 
10 -5 8 0O 7 O 5 
0 4 0 -7)\—2 4 0 


be a 4 x 4 square matrix interval. 
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A' = [Ai, B,] mar 
8 7 10 0O 2 q 
eee 4 1 -5 4 4 5 0 
[ Ai.B: | = . 
3 2 8 0 3 Thy oS 
5 —2 0 -7 1 2 0O 


be the transpose of the 4 x 4 square matrix interval. 
Let A= [Pi, Ps] as (Qi, Q>] 


= [Pi + Qi, P2 + Qu] 


A, +A, B, +B, 
where P; = — 1 Q,= _—., 
1 5 Q: 5 
Ae ae =B' 
P, = —— and Q, = +1 
2 5 Q 5 
8 4 3 5 8 7 10 
To BO De 2 i‘ 4 1 -5 
10 -5 8 0O 3 2 8 
0 4 0 -7 5 —2 0 
Now P,; = 
2 
16 11 13 5 


Clearly P; is a symmetric matrix of order 4 x 4. 


1 1 


. —A' 
Now consider P; = ———- 


oh OS 
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8 4 3 5) (8 7 10 0 
Bee ee 
BO! 5 RO Se De aRe -10 
Oi de OG Ap) VE: SOF Oe SF 
Ora) a’ 
W370 Hh sé 
OL SF RO SO 
5 6 0 0 


_ B, +B, 
Qi = 5 
yea Bey 2 —_ 
1 5G. SD a 0 4 
2 0 *) 4 5 
yd 8 5.4 
4 4 10 -l 
tj 10-7 2 
DOF TO. 
=f: °9! sO 4G 


It is easily verified that Q; is a symmetric matrix of order 
four. 


= B, -B; 
Q 5 
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~ 
oO 

co NN WwW 
ios) 


oN oN 
& 


is a skew symmetric matrix of order four. 


16 11 13 5 


11 2 3 2 
Now A= 


2}13 -3 16 O 2 


5 2 0 -14 


Nl 
soe | 
ae 

ar 
so 
on 
on 
nN 


= [A), Bi] 


= [5 (Pj P)), 5 (Qi + Q»)]. 


| 
+ 
| 


RO 


=3 
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The sign column matrix interval is 


a, | | D, 
a b 
X=[X, XM]=]] 27 |}? 
a, b, 
where aj, b; € {-1, 1}; 1 <i<n. 
1) /-1 
-1 1 
1 1 
Y=] 1 |.) -1 
1 }/-l 
-1]] 1 
1 1 


=[Y,, Y2] is a sign column matrix interval or sign column 
vector interval. 


The sign row matrix interval sign row vector interval 
A= [Ai, Ao] 


= [(al, ..., Am), (b1, bo, ..., bm)] where a;, b; € {-1, 1}; 1 <i<m. 


Thus 
A=([(,-1, 1, 1, 1,-1,-1, 1), Cl, 1, 1,-1,-1, 1,-1, D] 
is a sign row vector interval or (sign row matrix interval). 


These sort of matrix intervals are useful in factor analysis. 


We will give examples of orthogonal matrix interval. 
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Let 
1 -l 2 1 
X={/1 -2/)},)-1 2 
1 3 0 1 


is again a matrix interval. 


Consider 
1 -l 2 1 
; 1 1 1)\)/2 -1 0 
x'xX= : x}}1 —2],)-1 2 
-1 —2 3/\1 2 #1 
1 3 0 1 
1 1 2 


{0 dG 3h 


Thus this is an example of a orthogonal matrix interval. 
If we have X = [A, B] be a matrix interval. 
X' = [A‘, B'] be the transpose of X. 


Suppose X'X = [A‘A, B'B] = [I,, In] then we call x to be a 
orthogonal interval matrix. 
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Consider 


a 
oor So 


to be the matrix interval. 


xt = 1 00 0\)/(/0 0 0 1 
00° tO) 0° D0 6 


be the transpose of the matrix interval X. 


: 1 00 0)//0 0 0 1 
xXxX= , x 
001 0)/\0 1 0 0 


ooo KF 
o°o = © 


= [I>, I,]. 
Thus X is an orthonormal matrix interval. 
Suppose X = [A, B] be a matrix interval and X' = [A‘, B'] be 
the transpose of the matrix interval; if X'X = [A‘A, B'B] = [In, 


In] then we say X is a orthonormal matrix interval. 


Interested reader can give examples of them. 
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This result can easily be extend to square matrix interval. If 
X = [A, B] be an x n square matrix interval and X' = [A', B'] be 
the transpose of X be such that X'X = [I,, In] then we define X 
to be a orthogonal matrix interval or more specifically X is a 
orthonormal matrix interval. Suppose X'X = [D, C] where D 
and C are just diagonal matrices or symmetric matrices then we 
define X to be only a orthogonal square matrix interval. 


Thus we use matrix theory in a new way to study properties 
of interval matrices constructed using natural class of intervals. 


Chapter Six 


DSM MATRIX OF REFINED LABELS 


[46, 52] have worked with the new notion of DSm vector 
spaces and DSm super vector spaces. Lz = {..., Lj, ..., L-4, Lo, 
Li, Lo, ..., Lj, ...} = {Lj | J © Z} the set of extended labels with 
positive and negative indexes. 


Similarly Lg = Lg | q € Q} as set of labels whose indexes 
are fractions. Lg is isomorphic to Q, through the isomorphism 


fg (Lq) = 


4 


m+l1 
for any q € Q. 


On similar lines they define Lp = {Lj | j € R}; R the set of 
real numbers. 


La is isomorphic with R through the isomorphism 
r 
fr(L,) = 


m+1 
for any r € R. For more refer [46, 52]. 
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Further they have proved {Lr, +, x} 1s a field where + is the 
vector addition of labels and x is the vector multiplication of 
labels defined as DSm field of refined labels. 

If [L;, L.] is a label interval then L3,. = Lj.5 is the label in 
the middle of the label interval. 


Also L_; =—L; that occur in qualitative calculations. We just 
from [ ] recall the operations. 


Let a, b, c € R and the labels 


a= »Lo= 2 and L, = e : 
m+l m+l m+l 


Vector addition of labels 


Lz + Lp = Lath = L 


a+b ° 
m+1 


In this chapter we just indicate how we have used matrix 
theory in the construction of refined label matrices. When the 
refined labels happen to be ordered that is 


Lo <L,<...<Ln 
m ¢€ N; then we can use these matrices to built vector spaces, 
find eigen values and all properties pertaining to vector spaces 


can be derived. This has been elaborately carried out in [46, 
S21: 


So this is an unique and a new way of using matrix theory. 

Not only it is matrix of refined labels whose properties are 
studied, we have also invented to notion of DSm super vector 
spaces, which uses refined labels of super matrices. 


This study is elaborately made in [46, 52]. 


Chapter Seven 


N-MATRICES AND THEIR APPLICATIONS 


In this chapter we proceed onto recall the definition, 
properties and applications of n-matrices. This is yet another 
innovative method of using matrices. For more about these 
concepts please refer [47, 51, 54, 56]. 


DEFINITION 7.1: Let A = A; U A> where A; and A> are two 
matrices which are distinct and the entries are taken from Q or 
Ror Zor Z, or C(Z,). 


If both A; and Az are row matrices we call A to be a birow 
matrix. Similarly, if both A; and Az are column matrices we call 
A to be a bicolumn matrix, A is defined as the square bimatrix if 
both A; and Az are square matrices. If A, and A) are 
rectangular matrices we define A to be a rectangular bimatrix. 
Ifin A = A; UV A2 one of A; is a row matrix and A; is not a row 
matrix we define A to be a mixed bimatrix. That is if A; and A 
are two different matrices of order m xn ands xt we call A to 
be mixed bimatrix. 


We will illustrate this by some examples. 
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Example 7.1: Let 
A=A,U A) =(3, 1, 5, 6, 0, 2) U (3, 1 0, 0, 5), 


we call A to be a row bimatrix. 


Example 7.2: Let 


A=A,;UA)= 


NN NK NY WD 
G 


be a column bimatrix. 


Example 7.3: Let 
B=B, UB) =(9, 0, 1) U (8, 4, -4), 


B is arow bimatrix. 


Example 7.4: Let 


4 2 7 
9 0 1 2 
A=A,UA)= U!l0O 1 9 -1 8 
0 1 0 -4 
5 6 8 1 -4 


A is a rectangular bimatrix. 
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Example 7.5: Let 


0 -1 2 
3 4 4 7 0 3 
A= Ay U Ad = U 
-1 0 -1 2 0 1 
0 1 5 8 
be a square bimatrix. 
Example 7.6: Let 
PS 
-6 
M =(8, 0,-1, 4, 3, 3,7) U | 8 
1 
| —4 | 
be a mixed bimatrix. 
Example 7.7: Let 
~9 
2 
9 0 2 3 
T=T, UTh,= U 
-1 0 3 4 
—5 
6 


be a mixed bimatrix. 
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Example 7.8: Let 


9 

—2 

3 9 0 
S=S,;,US%= U 

0 -1 2 

3 

0 


be the mixed bimatrix. 


Example 7.9: Let 


3 
P=P, UP)=(-3, 0, 1,2)U 
eer ee 


be the mixed bimatrix. 
Now we can have bimatrices both mixed or otherwise. 


Suppose A = A; U Ap UA3U... U Ayj « > n= 2 and A;’s 
are distinct row matrices we define A to be n-row matrix. If 
n = 2 it becomes a row bimatrix. If n = 3 we get the row 
trimatrix and so on. 


If row matrices are replaced by column matrices we call 
them as n-column matrix. Instead of a column matrix we use 
rectangular matrices we call them as n-rectangular matrix or 
rectangular n-matrix. 


Suppose in A = A; U Ap U ... UA, the A;’s are square 
matrices, we define A to be a n-square matrix or a square n- 


matrix. 


We will illustrate all these situations by examples. 


n- Matrices and their Applications | 209 


Example 7.10: Let A= A, U Ar U A3 U Ag U As = (8, —1, 0, 2, 
1) U (0, 5, 8) UG, 2) U C7, 1, 2, 3, 4, 5, -7) U (4, 6, 8, 10, 12) 
be a 5-row matrix or row 5-matrix. 


Example 7.11: Let 


0 
—4 
1 
8 5 
, 9 2 -8 
M=M,U MUM; U My= 3 U 7 UlL3}U 
—4 
0 
5 
2 
- -7 


be the column 4-matrix. 


Example 7.12: Let V= V,UV2U V3 U V4U V5 U V6 U V7 


8 0 1 2 
1 2 3 
8 7 3 0 05 6 7 
- Ul4 5 6 U 
0 1 0 -5 007 8 
7 0 1 
00 0 -5 
1 72 3 -4 5 
9 0 0 0 
-6 7 -8 9 0O 1 00 
0 8 1 2 
-l1 2 +3 4 -5}/U/2 3 O]U 
00 5 0 
-7 8 -9 6 4 5 6 
00 0 7 
1 2 5 8 


be a 7-square matrix or square 7-matrix. 
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Example 7.13: LetT=T, UT2UT3 UT, U Ts 


123 4 5 6 

8 05 0 -8 -5 1 2 
= Ul6 5 4 2 1 

1207 5 0 -7 0 
9 12 8 3 7 
8 2 O} | = 
8 0 12 0 1 4 

3 2 -1 0 5 
5 2 0101 0 

1 5 4 2 3 
0 1:°20%1i21~=41 

VIO 1T}]U}5 6 7] U 
6 -1 05 0 8 7 

4 2 8 9 4 
5 2 10 7 0 6 

0 1 0 0 1 
| 6 2.6, 1 10) 

1 5 0 


be a rectangular 5-matrix or 5-rectangular matrix. 


Example 7.14. Let V=V,; UV2U V3U Va 


7 


34 

9 =i 

=(-9,0, ¥3,1,7+i)U | 8-i | | 
5i 

4+8i | 


8 -1 7-1 8+41 


be a 4-mixed matrix or mixed 4-matrix. 
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Example 7.15: Let A= A; U Ap U A3 U Ag U As U Ag 


9 8 
8 -i 8 0 -1 3 2 0 
=(-8,0,1,3+21)U|}4 5]/U;/-1 6 0 2 4 O 
0 6 2 4 7 0 5 -7 
g:-=2 | 
7 | ; 
. 3 1 441 
i+8 
3 8 4 2 7 
U|] 7-1 ]/U . 
0 4 5 0 8-1 
4+3 
0 —7 4 
7421 


be a 6-mixed matrix or mixed 6-matrix. 


Now we have seen just definition and examples of them. 
We proceed onto give three major applications of them. 


The first application of bimatrices (n-matrices) is in the 
construction of bilinear algebra (n-linear algebra of type I and 
type II). For more literature please refer [60-1]. 


One of the major application of these n-matrices (bimatrices 
is that they can be used in the construction of n-eigen values 
(bieigen values) and n-eigen vectors (bieigen vectors). This 
bulk operation helps in saving time as well as easy for 
calculation and comparison. Further bimatrices are used in the 
construction of bifuzzy models like BiFCMs, BiFRMs, n-FCMs 
and n-FRMs; n = 2. 


We just indicate how the bieigen values (n-eigen values) are 
calculated. We only give illustrations for theory please refer 
[60-1]. 


212 | Innovative Uses of Matrices 


Example 7.16: Let 


af 8 0 O 
M=MiUM=| ir 0 -1 2 
0 6 | 


be a bimatrix associated with some linear operator. The bieigen 
values associated with M is as follows: 


|M — (A; Ke) Ay) = |M, = Al U |M> — hl 


Sh, oO 0 
3-2, 0 
= || GO: Sle, 2 
ge ho 
0 6. 43, 


= (3 — Ay) 1 —Ay) U (8-Ay) [1-23 ) - 17] =0 U0 
= — (I-A) (3-Ay) U -(8-Az) ((1-A3 +1’) 
=0U0 


leading to A; =3 and—1 and A, =8 2A, =+ V13. 


Thus 4=3U 8,3 Ut V13,3U—-V13 


tire al Cats al ie 3 


Thus using these six bieigen values one can calculate the 
related bieigen vectors. By this way we can compare the eigen 
values and also perform the operations simultaneously. 


However the entries of the matrix M are taken from the 
field of reals. If they are taken from the rationals all the 
solution will not exist. We may have partial solution in that 
case. 
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Example 7.17: Let P=P, UP: UP3 


a OR 
r- ON OC 
oO wkW Oo SO 
ho Oo SO 


be the trimatrix with the entries from reals. To find the trieigen 
values and trieigen vectors of P. 


Consider |P — A (=A, U Ay U A3)| 


= |P — || U |P> —),| U |P3 —i)3| 


2-’, 0 
32%, +1 
= Ul 0 1-a, 0 
0 5-A, 
0 3-A, 


= (3-A1) (5-A1) U (2-Ag) (1-A2) x 
(3-Az) U (1-As) (2-As) (3-As) (4-As) 


Thus the trieigen values are (3, 2, 1), (3, 2, 2), (3, 2, 3), 
(3, 2, 4), (3, 1, 1), (G3, 1, 2), GB, 1, 3), GB, 1, 4), GB, 3, 1), GB, 3, 2), 
(3, 3, 3), (3, 3, 4), (5, 2, 1), (5, 2, 2), (5, 2, 3), (5, 2, 4), (5, 1, 2), 
(5, 1, 3), (5, 1, 4), (5, 1, 1), (5, 3, 1), (5, 3, 2), (5, 3, 3), (5, 3, 4). 


We have 24 sets of trieigen values and their corresponding 
24 set of trieigen vectors can be calculated. This is one of the 
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major innovative method of applying the bimatrix (n-matrix) 
theory in bilinear algebra or (n-linear algebra). 

The next set of applications of bimatrices (n-matrices) is in 
the fuzzy models. We just indicate how these bimatrices (n- 
matrices) used. We call a bimatrix A = A; U A; to be a fuzzy 
matrix if the entries of both A; and A; are from the unit interval 
[0, 1]. We make use of the min max or max min or min min or 
max max operations on them. 


On similar lines we define n-fuzzy matrix. A matrix S = S, 
US. U... US, (n > 2) is said to be a n-fuzzy matrix if the 
entries of each S; is from the unit interval [0, 1]; 1 <i<n. We 
use min max or max min or other such operators on S;’s. We 
also at times include —1 as a fuzzy number depending on the 
fuzzy model. The fuzzy models which make use n-fuzzy 
matrices are n- FCMs, n-FRMs and n-FAMs. 


For more about these structures one can refer the books 
[48, 54, 56]. 


Suppose we have a problem related with the school 
students. 


The school student problem is dependent on the teachers 
also dependent on the parents and on the school management. 
So the school children problem is dependent on three groups. 
Thus we can use a fuzzy trimatrix and use the Fuzzy Relational 
Maps (FRMs) model to analyse the problem. 


Let us assume S), So, ..., S¢ are the attributes associated 
with the FRM. 


Let T,, To, ..., T7 be the attributes related with the teachers. 
P,, Po, ..., P¢ are the attributes associated with the parents 


and M;, Mb, ..., Mg are the attributes associate with the 
management. 


n- Matrices and their Applications | 215 


We get using an experts opinion a trigraph and the matrix 
associated with the trigraph will be the fuzzy trimatrix which 
serves as the dynamical system of the triFRM model. If S 
denotes the triFRM models tridynamical system then S = S; U 
So U S3 


TG i PB BP. M, M,... M, 
1 S; S; 
= Sy US, US, 
Ss Ss S; 


The entries of the matrices S;, S. and S3 are from the unit 
interval [0, 1]. 


Any state vector X = X; U X) U X3 = 
Se ne Re VAIO ook UU nc Ke) 
where x € {0,1}; 1 <i<6,1<j<3. 


We find X 0 S= X; 0 8; UX. 082 U X30 83 
=Y,UY2U Y3 
=Y. 
where elements of Y; are from {0, 1} Y; will be (y, ... y;) Yo 
will be (y;,y3, ---» Y¢) and Y3 will be (y;,y3, .-.5 Ye) y! E 
{0, 1}, 1 <j<3and1<i<8. 


Thus these n-fuzzy matrices are used in constructing 
n-fuzzy models, when we have n-experts or these models can 
also be defined as multi expert models. 
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Finally we see these n-matrices when constructed over Zo ; 


p a prime and « > n > 1, can be used in the construction of n- 
codes. 

We have defined n-codes [49]. For more about the uses of 
n-matrices in algebraic coding theory refer [49]. 

The construction of n-codes using n-matrices happens to be 
one of the new ways of applying n-matrices to algebraic coding 
theory. 

Thus the concept of n-matrices and their applications are 
innovative contributions in matrix theory by the authors. 
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